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Abstract 

This article reviews the present status of the spin foam approach to the quantization of 
gravity. Special attention is payed to the pedagogical presentation of the recently introduced 
new models for four dimensional quantum gravity. The models are motivated by a suitable 
implementation of the path integral quantization of the Plebanski formulation of gravity on 
a simplicial regularization. The article also includes a self-contained treatment of the 2+1 
gravity. The simple nature of the latter provides the basis and a perspective for the analysis 
of both conceptual and technical issues that remain open in four dimensions. 
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Part I 

Introduction 
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1 Quantum Gravity: A Unifying Framework for Fundamen- 
tal Physics 

The revolution brought by Einstein's theory of gravity lays more in the discovery of the principle 
of general covariance than in the form of the dynamical equations of general relativity. General 
covariance brings the relational character of nature into our description of physics as an essential 
ingredient for the understanding of the gravitational force. In general relativity the gravitational 
field is encoded in the dynamical geometry of space-time, implying a strong form of universality 
that precludes the existence of any non-dynamical reference system - or non-dynamical background 
- on top of which things occur. This leaves no room for the old view where fields evolve on an rigid 
preestablished space-time geometry (e.g., Minkowski space-time): to understand gravity one must 
describe the dynamics of fields with respect to one another, and independently of any background 
structure. 

General relativity realizes the requirements of general covariance as a classical theory, i.e., for 
h = 0. Einstein's theory is, in this sense, incomplete as a fundamental description of nature. 
A clear indication of such incompleteness is the generic prediction of space-time singularities in 
the context of gravitational collapse. Near space-time singularities the space-time curvature and 
energy density become so large that any classical description turns inconsistent. This is reminiscent 
of the foundational examples of quantum mechanics - such as the UV catastrophe of black body 
radiation or the instability of the classical model of the hydrogen atom - where similar singularities 
appear if quantum effects are not appropriately taken into account. General relativity must be 
replaced by a more fundamental description that appropriately includes the quantum degrees of 
freedom of gravity. 

At first sight the candidate would be a suitable generalization of the formalism of quantum 
field theory (QFT). However, the standard QFT's used to describe other fundamental forces are 
not appropriate to tackle the problem of quantum gravity. Firstly, because standard QFT's are 
not generally covariant as they can only be defined if a non-dynamical space-time geometry is 
provided: the notion of particle, Fourier modes, vacuum, Poincare invariance are essential tools 
that can only be constructed on a given space-time geometry. This is a strong limitation when 
it comes to quantum gravity since the very notion of space-time geometry is most likely not 
defined in the deep quantum regime. Secondly, quantum field theory is plagued by singularities 
too (UV divergences) coming from the contribution of arbitrary high energy quantum processes. 
This limitation of standard QFT's is expected to disappear once the quantum fluctuations of the 
gravitational field, involving the dynamical treatment of spacetime geometry, are appropriately 
taken into account. But because of its intrinsically background dependent definition, standard 
QFT cannot be used to shed light on this issue. A general covariant approach to the quantization 
of gravity is needed. 

This is obviously not an easy challenge as in the construction of a general covariant QFT one 
must abandon from the starting point most of the concepts that are essential in the description 
of 'no-gravitational' physics. One has to learn to formulate a quantum theory in the absence of 
preferred reference systems or pre-existent notion of space and time. Loop quantum gravity (LQG) 
is a framework to address this task. The degrees of freedom of gravity arc quantized in accordance 
to the principles of general covariance. At the present stage there are some indications that both 
the singularity problems of classical general relativity as well as the UV problem of standard QFT's 
might vanish in the framework. 

However, these indications are not conclusive mainly because the systematic description of 
dynamical processes remain an open problem at the present stage. In this article we review the 
status of the construction of a promising approach for handling the difficult dynamical question in 
LQG: the spin foam formulation. 
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1.1 Why non-perturbative quantum gravity? 

Let us make some observations about the problems of standard perturbative quantum gravity. In 
doing so we will revisit the general discussion above, in a special situation. In standard perturbative 
approaches to quantum gravity one attempts to describe the gravitational interaction using the 
same techniques applied to the definition of the standard model. As these techniques require 
a notion of non dynamical background one (arbitrarily) separates the degrees of freedom of the 
gravitational field in terms of a background geometry r\ a \, for a.b — 1 • ■ • 4 - fixed once and for all 
- and dynamical metric fluctuations h a i,. Explicitly one writes the spacetime metric as 

9ab = T]ab + Kb- (1) 

Notice that the previous separation of degrees of freedom has no intrinsic meaning in general 
relativity. In other words, for a generic spacetime metric g a b we can write 

9ab = Vab + Kb = Vab + Kb, (2) 

where and rj a b can lead to different background light-cone structures of the underlying spacetime 
(M,g a b); equally natural choices of flat background metrics lead to different Minkowski metrics in 
this sense. This is quite dangerous if we want to give any physical meaning to the background, 
e.g., the light cone structures of the two 'natural' backgrounds will be generally different providing 
different notions of causality! Equation ([l} is used in the classical theory in very special situations 
when one considers perturbations of a given background 77^. In quantum gravity one has to deal 
with arbitrary superpositions of spacetimes; the above splitting can at best be meaningful for very 
special semi-classical states 'peaked', so to say, around the classical geometry r/ a b with small fluc- 
tuations, i.e., it might only make sense only a posteriori in the framework of a quantum gravity 
theory defined independently of any background geometry. Consequently, even when perturba- 
tions should be remain useful to approximate special situations, the splitting in Equation (TTJ) is 
inconsistent when considering general states with their arbitrary quantum excitations at all scales. 

This is specially so in gravity due to the dual role of the gravitational field that simultaneously 
describes the geometry and its own dynamical degrees of freedom. More explicitly, in the standard 
background dependent quantization the existence of a fixed background geometry is fundamental 
in the definition of the theory. In the usual treatment one imposes fields at space-like separated 
points to commute alluding to standard causality considerations. Even when this is certainly 
justified for matter fields (at least in the range of applicability of the standard model) it makes no 
sense in the case of gravity: one would be using the causal structure provided by the unphysical 
background r] a b. This difficulty has been raised several times (see, for instance, 1361]). 

Therefore, treating the gravitational field according to the splitting given in Equation (fTJ) is 
inconsistent with the fundamental nature of the gravitational degrees of freedom. Treating h a {, with 
the standard perturbative techniques is equivalent to viewing h a b as another matter field with no 
special properties. As mentioned above, gravitons would propagate respecting the causal structure 
of the un physical background to all orders in perturbation theory^ Radiative corrections do not 
affect causality. Even though such thing is quite sensible when dealing with matter fields in the 
regime where gravity is neglected it would be clearly wrong for the gravitational interaction. 

1.2 Non-renormalizability of general relativity 

If we ignore all these issues and try to setup a naive perturbative quantization of gravity we find 
that the theory is non renormalizable. This can be expected from dimensional analysis as the 

1 One can in principle avoid this by carefully devising a theory of perturbations for generally covariant systems. 
For some ideas on this see 147, 233 . 
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Figure 1: The larger cone represents the light-cone at a point according to the ad hoc background 
rjab. The smaller cones are a cartoon representation of the fluctuations of the true gravitational 
field represented by g a b- Gravitons respect the causal structure of the un physical metric to all 
orders in perturbation theory. 

quantity playing the role of the coupling constant turns out to be the Planck length £ p . The non 
renormalizability of perturbative gravity is often explained through an analogy with the (non- 
renormalizable) Fermi's four fermion effective description of the weak interaction. Fermi's four 
fermions theory is known to be an effective description of the (renormalizable) Weinberg-Salam 
theory. The non renormalizable UV behavior of Fermi's four fermion interaction is a consequence 
of neglecting the degrees of freedom of the exchanged massive gauge bosons which are otherwise 
disguised as the dimension- full coupling A Fcrmi rj at momentum transfer much lower than 

the mass of the W particle (q 2 << rrtyy). A similar view is applied to gravity to promote the 
search of a more fundamental theory which is renormalizable or finite (in the perturbative sense) 
and reduces to general relativity at low energies. From this perspective it is argued that the 
quantization of general relativity is a hopeless attempt to quantizing a theory that does not contain 
the fundamental degrees of freedom. 

These arguments, based on background dependent concepts, seem at the very least questionable 
in the case of gravity. Although one should expect the notion of a background geometry to be 
useful in certain semi-classical situations, the assumption that such structure exists all the way 
down to the Planck scale is inconsistent with what we know about gravity and quantum mechanics. 
General considerations indicate that standard notions of space and time are expected to fail near 
the Planck scale ij^. In the field of loop quantum gravity people usually tend to interpret the 
severe divergences of perturbative quantum gravity as an indication of the inconsistency of the 
separation of degrees of freedom in Equation (fTJ). According to concrete results in LQG the nature 
of spacetime is very different from the classical notion in quantum gravity. The treatment that 
uses Equation (TTJ) as the starting point is assuming a well defined notion of background geometry 
at all scales which directly contradicts these results. 

It is possible that new degrees of freedom would become important at more fundamental scales. 
It is also possible that including these degrees of freedom might be very important for the consis- 
tency of the theory of quantum gravity. However, there is a constraint that seem hardly avoidable: 

2 For instance a typical example is to use a photon to measure distance. The energy of the photon in our 
lab frame is given by E-y = hc/X. We put the photon in a cavity and produce a standing wave measuring the 
dimensions of the cavity in units of the wavelength. The best possible precision is attained when the Schwarzschild 
radius corresponding to energy of the photon is of the order of its wavelength. Beyond that the photon can collapse 
to form a black hole around some of the maxima of the standing wave. This happens for a value A c for which 
A c Ri GE-y/c 4 = hG/(\ c c 3 ). The solution is A c ss ^JhG /c 3 which is Planck length. 
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if we want to get a quantum theory that reproduces gravity in the semi-classical limit we should 
have a background independent formalism. In loop quantum gravity one stresses this viewpoint. 
The hope is that learning how to define quantum field theory in the absence of a background is a 
key ingredient in a recipe for quantum gravity. 

One of the achievements of the background independent formulation of LQG is the discovered 
intrinsic discreteness of quantum geometry. Geometric operators have discrete spectra yielding a 
physical mechanism for the cut-off of UV degrees of freedom around the Planck scale. In this way 
LQG is free of UV divergences. Is this a complete answer to the non-renormalizability problem 
of gravity mentioned above? Unfortunately, the answer to this question is in the negative at this 
stage. The reason is that the true answer to the question of renormalizabity or not renormalizability 
is not in the presence of divergences but in the degree of intrinsic ambiguity of the quantization 
recipe applied to a field theory. 

In standard background dependent quantum field theories, in order to avoid UV divergences 
one has to provide a regularization prescription (e.g., an UV cutoff, dimensional regularization, 
point splitting, etc.). Removing the regulator is a subtle task involving the tuning of certain terms 
in the Lagrangian (counter-terms) that ensure finite results when the regulator is removed. In fact 
by taking special care in the mathematical definition of the products of distributions at the same 
point, in good cases, one can provide a definition of the quantum theory which is completely free 
of UV divergences |165j (see also [3361 12181 1219] ). However, any of these regularization procedures 
is intrinsically ambiguous. The dimension of the parameter space of ambiguities depends on the 
structure of the theory. The right theory must be fixed by comparing predictions with observations 
(by the so-called renormalization conditions). In loop quantum gravity there is strong indications 
that the mathematical framework of the theory provides a regularization of divergences. It remains 
to settle the crucial issue of how to fix the associated ambiguities. 

According to the previous discussion, ambiguities associated to the UV regularization allow for 
the classification of theories as renormalizable or non nonrenormalizable quantum field theories. 
In a renormalizable theory such as QED there are finitely many ambiguities which can be fixed 
by a finite number of renormalization conditions, i.e., one selects the suitable theory by appropri- 
ate tuning of the ambiguity parameters in order to match observations. In a nonrenormalizable 
theory (e.g., perturbative quantum gravity) the situation is similar except for the fact that there 
are infinitely many parameters to be fixed by renormalization conditions. As the latter must be 
specified by observations, a non-renormalizable theory has very limited predictive power. In the 
case of gravity there is evidence indicating that the theory might be non-perturbatively renormal- 
izable coming from the investigations of non trivial fix points of the renormalization group flow in 
truncated models (see (31 1J and references therin). 

Removing UV divergences by a regularization procedure and ambiguities are two sides of the 
same coin quantum field theory. Although this can happen in different ways in particular for- 
mulations, the problem is intrinsic to the formalism of quantum field theory. In this respect, it 
is illustrative to analyze the non-perturbative treatment of gauge theories in the context of lat- 
tice gauge theory (where the true theory is studied by means of a regulated theory defined on a 
space-time discretization or lattice). It is well known that here too the regulating procedure leads 
to ambiguities; the relevance of the example resides in the fact that these ambiguities resemble 
in nature those appearing in loop quantum gravity. More precisely, let us take for concreteness 
SU(2) Yang-Mills theory which can be analyzed non-perturbatively using the standard (lattice) 
Wilson action 



In the previous equation U p £ SU(2) is the holonomy around plaquettes p, and the sum is over 
all plaquettes of a regulating (hyper-cubic) lattice. It is easy to check that the previous action 
approximates the Yang-Mills action when the lattice is shrunk to zero for a fixed smooth field 




(3) 
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configuration. This property is referred to as the naive continuum limit. Moreover, the quantum 
theory associated to the previous action is free of any UV problem due to the UV cut-off provided 
by the underlying lattice. 

Is this procedure unique? As it is well known the answer is no: one can regulate Yang-Mills 
theory equally well using the following action instead of Equation ([3]) : 



stl oci£(i- ^^^[n^^ + n^^)]) , (4) 



where H (m) (U p ) denotes to the SU(2) unitary irreducible representation matrix (of spin m) evalu- 
ated on the plaquctte holonomy U p . Or more generally one can consider suitable linear combina- 
tions 

^lym — ^ ^ Urn ^LYM- (5) 
n I 

From the view point of the classical continuum theory all these actions are equally good as they all 
satisfy the naive continuum limit. Do these theories approximate in a suitable sense the continuum 
quantum field theory as well? and are these ambiguities un-important in describing the physics of 
quantum Yang-Mills theory? The answer to both of these questions is yes and the crucial property 
that leads to this is the renormalizability of Yang-Mills theory. Different choices of actions lead 
indeed to different discrete theories. However, in the low energy effective action the differences 
appear only in local operators of dimension five or higher. This simple dimensional argument 
leads to the expectation that in the continuum limit (i.e., when the regulating lattice dependence 
is removed by shrinking it to zero) all the above theories should produce the same predictions in 
the sense that can safely ignore non-renormalizable contributions. This is indeed the case in usual 
cases, and thus the ambiguities at the 'microscopic level' do not have any effect at low energies 
where one recovers quantum Yang-Mills theory. 

Similar ambiguities are present in LQG, they appear in the regularization of the constraints in 
the canonical formulation 297 . These ambiguities are also present in the spin foam formulation 
which is the central subject of this review. The simple dimensional argument given for Yang-Mills 
theory is no longer valid in the case of gravity, and a full dynamical investigation of their effect in 
physical predictions must be investigated. The spin foam approach is best suited for getting the 
necessary insights into this important question. 



2 Why Spin Foams? 

We now provide the basic motivation for the study and definition of spin foams. We present this 
motivation from the perspective of the canonical quantization approach known as loop quantum 
gravity (see [315, 353, 31 J. In such context spin foams appear as the natural tool for studying 
the dynamics of the canonically defined quantum theory from a covariant perspective. In order to 
introduce the spin foam approach it is therefore convenient to start with a short introduction of 
the quantum canonical formulation. 



2.1 Loop quantum gravity 

Loop quantum gravity is an attempt to define a quantization of gravity paying special attention 
to the conceptual lessons of general relativity. The theory is explicitly formulated in a background 
independent, and therefore, non perturbative fashion. The theory is based on the Hamiltonian 
(or canonical) quantization of general relativity in terms of variables that are different from the 
standard metric variables. In terms of these variables general relativity is cast into the form of a 
background independent SU (2) gauge theory whose phase space structure is similar in various ways 



7 



to that of SU (2) Yang-Mills theory (the key difference being the absence of any background metric 
structure). The main prediction of loop quantum gravity (LQG) is the discreteness [337] of the 
spectrum of geometrical operators such as area and volume. The discreteness becomes important 
at the Planck scale while the spectrum of geometric operators crowds very rapidly at 'low energy 
scales' (large geometries). This property of the spectrum of geometric operators is consistent with 
the smooth spacetime picture of classical general relativity. 

Thus, from the perspective of LQG, it is not surprising that perturbative approaches would 
lead to inconsistencies. In splitting the gravitational field degrees of freedom as in Equation (fT|) 
one is assuming the existence of a background geometry which is smooth all the way down to the 
Planck scale. As we consider contributions from 'higher energies', this assumption is increasingly 
inconsistent with the fundamental structure discovered in the non perturbative treatment. How- 
ever, despite of the many achievements of LQG there remain important issues to be addressed. At 
the heart of the completion of the definition of the theory the clear-cut definition of the quantum 
dynamics remains open. The spin foam approach is one of the main avenues to exploring this 
problem. 

The background independence of general relativity implies that the canonical formulation of 
the field theory is that of a gauge theory with diffeomorphism as part of the gauge group. LQG 
is constructed by quantizing a phase space formulation of general relativity in terms of SU (2) 
connection variables. This introduces an extra SU (2) gauge symmetry group. The presence of 
gauge symmetries implies the existence of relations among phase space variables - defined on a 
spacelike initial value hypersurface - known as constraints. These constraints define the Poisson 
algebra of infinitesimal generators of gauge transformations. There are three local constraints 
G % - the Gauss constraints - generating SU(2) gauge transformations, three local constraints V a 
- the vector constraints - generating three dimensional diffeomorphisms of the initial spacelike 
hypersurface, and finally a scalar local constraint S related to the remaining gauge symmetry 
related to the four-diffeomorphism symmetry of the Lagrangian formulation. 

The canonical quantization of systems with gauge symmetries is often called the Dirac pro- 
gram. The Dirac program |137| 1217] applied to the quantization of general relativity in connection 
variables leads to the LQG approach. 

The first step in the recipe consists in finding a representation of the phase space variables of 
the theory as operators in a kinematical Hilbert space T-Lkin satisfying the standard commutation 
relations, i.e., { , } — > —i/fr[ , ]• This step has been succesfully completed in LQG. One chooses 
the polarization where the SU(2) connection is the configuration variable. Unconstrained phase 
space variables are replaced by the so-called holonomy-flux algebra which is represented by associ- 
ated operators in a kinematical Hilbert space of suitable functionals of the generalized-connection 
ip[A] which are square integrable with respect to the so-called Ashtekar-Lewandowski [28] (gauge 
invariant and diffeomorphism invariant) measure d/j. AL [A]. A key input is the use of the holonomy 
flux algebra as a starting point for quantization. Many peculiar properties of LQG follow from this 
choice this choice is motivated by having an observable algebra leading to simple (diffeomorphism 
covariant) Poisson brackets H. 

The kinematical inner product is given by 



The representation of the basic unconstrained phase space variables as suitable operators in Hkin 
used in LQG - which has the additional key property of containing a special state that is diffeo- 
morphism invariant - has been shown to be unique [239 . 

3 Interestingly, and as a side remark, the structure of the phase space of gravity with boundaries used in the 
description of concrete physical models seem to provide an extra justification of the use of the holonomy flux 
variables (see section IV E in 1 1591 for a more detailed discussion of this intriguing fact). 




(6) 
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The next step is to promote the constraints to (self-adjoint) operators in Hkin m such a way 
that the classical Poisson algebra is respected by the appropriate quantum commutator algebra (if 
this last part is achieved the quantization is said to be non- anomalous). In the case of gravity one 
must quantize the seven constraints Gi, V a , and S. Both the Gauss constraint and (finite) diffeo- 
morphism transformations have a natural (unitary) action on states on 1-Lkin- For that reason the 
quantization (and subsequent solution) is rather direct. The quantization of the scalar constraint 
presents difficulties. Concrete quantizations producing well defined operators are available: no UV 
divergences are encountered by these proposals; the fundamental discreteness of quantum geometry 
in LQG plays the role physical regulator at the Planck scale [348]. Despite of this partial success, 
problems partly related to special features of the constraint algebra (field dependent structure con- 
stant) and partly related to the non-existence of an infinitesimal generator of diffcomorphism at 
the quantum level make the issue of whether any of the proposed quantizations of the constraints 
are anomaly-free a difficult open question. 

In the third step one needs to characterize the space of solutions of the constraints and define 
the corresponding inner product that defines a notion of physical probability. This defines the 
so-called physical Hilbert space H ph . In LQG physical states are those who satisfy the quantum 
constraints which, in this sense, could be called Quantum Einstein's equations. More precisely 

Gi(A,E)\V >=0 
V a (A,E)\9 >=0, 

S{A,E)\t> >=0. (7) 

The space of solutions of the first six equations is well understood. The space of solutions of 
quantum scalar constraint remains open. For some mathematically consistent definitions of S 
the characterization of the solutions is well understood |31j . The definition of the physical inner 
product is still an open issue. The completion of this step is intimately dependent on the previously 
mentioned consistency of the quantum constraint algebra. 

Finally, physical interpretation necessitates the definition of a (complete) set of gauge invariant 
observables, i.e., operators commuting with the constraints. They represent the questions that are 
well-posed in a generally covariant quantum theory. Already in classical gravity the construction 
of gauge independent quantities is a subtle issue. At the present stage of the approach physical 
observables are explicitly known only in some special cases. Understanding the set of physical 
observables is however intimately related with the problem of characterizing the solutions of the 
scalar constraint described before. 

The spin foam approach was constructed as a means to tackle the difficult question of dynamics 
and the definition of observable quantities in LQG. It is an attempt to address the difficulties 
associated last three steps of the list given above by rethinking the problem of dynamics from a 
path integral covariant perspective. So we need first to briefly discuss the special features of the 
path integral formulation in the case of generally covariant systems. 

2.2 The path integral for generally covariant systems 

LQG is based on the canonical (Hamiltonian) quantization of general relativity whose gauge sym- 
metry is diffeomorphism invariance. In the Hamiltonian formulation the presence of gauge symme- 
tries |137j gives rise the constraints mentioned in the previous section. In the present discussion we 
will schematically represent the latter by the equation C(p,q) = for (p, q) G T. The constraints 
restrict the set of possible states of the theory by requiring them to lay on the constraint hyper- 
surface. In addition, through the Poisson bracket, the constraints generate motion associated to 
gauge transformations on the constraint surface (see Figure ([2])). The set of physical states (the so 
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called reduced phase space T re d) is isomorphic to the space of orbits, i.e., two points on the same 
gauge orbit represent the same state in T re d described in different gauges (Figure [5]) ■ 

In general relativity the absence of a preferred notion of time implies that the Hamiltonian 
of gravity is a linear combination of constraints. This means that Hamilton equations cannot be 
interpreted as time evolution and rather correspond to motion along gauge orbits of general relativ- 
ity. In generally covariant systems conventional time evolution is pure gauge: from an initial data 
satisfying the constraints one recovers a spacetime by selecting a particular one-parameter family 
of gauge-transformations (in the standard ADM context this amounts to choosing a particular 
lapse field N(x,t) and shift N a (x,t)). 




Figure 2: On the left: the geometry of phase space in gauge theories. On the right: the quantization 
path of LQG (continuous arrows). 

From this perspective the notion of spacetime becomes secondary and the dynamical inter- 
pretation of the the theory seems problematic (in the quantum theory this is refered to as the 
"problem of time"). A reason for this is the central role played by the spacetime representation 
of classical gravity solutions useful only if can be interpreted via the notion of test observers (or 
more generally test fields j^. Due to the fact that this idealization is a good approximation to 
the (classical) process of observation the notion of spacetime is a useful concept only in classical 
gravity. 

As emphasized by Einstein with his hole argument (see |315] for a modern explanation) only the 
information in relational statements (independent of any spacetime representation) have physical 
meaning. In classical gravity it remains useful to have a spacetime representation when dealing with 
idealized test observers. For instance to solve the geodesic equation and then ask dif] -invariant- 
questions such as: what is the proper time elapsed on particle 1 between two succesive crossings 
with particle 2? However, already in the classical theory the advantage of the spacetime pic- 
ture becomes, by far, less clear if the test particles are replaced by real objects coupling to the 
gravitational fielcH. 

The possibility of using the notion of test observers and test fields to construct our interpreta- 
tive framework is no longer available in quantum gravity. At the Planck scale (t p ) the quantum 

4 Most (if not all) of the applications of general relativity make use of this concept together with the knowledge 
of certain exact solutions. In special situations there are even preferred coordinate systems based on this notion 
which greatly simplify interpretation (e.g., co-moving observers in cosmology, or observers at infinity for isolated 
systems). 

5 In this case one would need first to solve the constraints of general relativity in order to find the initial data 
representing the self-gravitating objects. Then one would have essentially two choices: 1) Fix a lapse N(t) and a 
shift N a (t), evolve with the constraints, obtain a spacetime (out of the data) in a particular gauge, and finally ask 
the diff-invariant- question; or 2) try to answer the question by simply studying the data itself (without t-evolution). 
It is far from obvious whether the first option (the conventional one) is any easier than the second. 
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fluctuations of the gravitational field become so important that there is no way (not even in prin- 
ciple!) to make observations without affecting the gravitational field. The notion of test fields and 
test observer is lost and with it the usual spacetime representation. In this context there cannot 
be any, a priori, notion of time and hence no notion of spacetime is possible at the fundamental 
level. A spacetime picture would only arise in the semi-classical regime with the identification of 
some subsystems that approximate the notion of test observers. 

What is then the meaning of the path integral in such background independent context? The 
previous discussion rules out the conventional interpretation of the path integral. There is no 
meaningful notion of transition amplitude between states at different times t\ > to or equivalently 
a notion of "unitary time evolution" represented by an operator U(t\ — to). Nevertheless, a path 
integral representation of generally covariant systems exists and arises as a tool for implementing 
the constraints in the quantum theory as we argue below. 

Due to the difficulty associated with the explicit description of the reduced phase space r re( j, in 
LQG one follows Dirac's prescription. One starts by quantizing unconstrained phase space T, repre- 
senting the canonical variables as self-adjoint operators in a kinematical Hilbert space Hkin- Poisson 
brackets are replaced by commutators in the standard way, and the constraints are promoted to 
self-adjoint operators (see Figure [2]). If there are no anomalies the Poisson algebra of classical 
constraints is represented by the commutator algebra of the associated quantum constraints. In 
this way the quantum constraints become the infinitesimal generators of gauge transformations in 
Hkin- The physical Hilbert space Jiphys is defined as the kernel of the constraints, and hence to 
gauge invariant states. Assuming for simplicity that there is only one constraint we have 

i>eU phys iff exp[iNd]\^) = |V>) V JVeR, 

where U(N) — exp[iNC) is the unitary operator associated to the gauge transformation generated 
by the constraint C with parameter N. One can characterize the set of gauge invariant states, and 
hence construct H p hys, by appropriately defining a notion of 'averaging' along the orbits generated 
by the constraints in Hkin- For instance if one can make sense of the projector 

P : Uun -> H ph y S where P := J dN U(N). (8) 

It is apparent from the definition that for any tp E Hkin then Pip 6 Uphys- The path integral 
representation arises in the representation of the unitary operator U(N) as a sum over gauge- 
histories in a way which is technically analogous to standard path integral in quantum mechanics. 
The physical interpretation is however quite different. The spin foam representation arises naturally 
as the path integral representation of the field theoretical analog of P in the context of LQG. Needles 
is to say that many mathematical subtleties appear when one applies the above formal construction 
to concrete examples (see [263 and references therein). In the second part of this review we will 
illustrate in complete detail all the features of the path integral approach for generally covariant 
systems, formally described here, in the 2+1 pure gravity solvable model. In four dimension the 
spinfoam approach is a program in progress where various issues still remain open. We will describe 
the various models in the first part of this article. 

2.3 A brief history of spin foams in four dimensions 

In this section, we briefly describe the various spin foam models for quantum gravity in the liter- 
ature (for previous reviews on spin foams see 287, 35 , 296, 295, 242, 17). 

6 In order to make a Planck scale observation we need a Planck energy probe (think of a Planck energy photon) . 
It would be absurd to suppose that one can disregard the interaction of such photon with the gravitational field 
treating it as test photon. 
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2.3.1 The Reisenberger model 

According to Plebanski [301] the action of self dual Riemannian gravity can be written as a con- 
strained SU{2) BF theory 



S(B,A) = J Tr [£? A F(A)} — ipi 



B l A B J - -S ij B k A B k 



3 



(9) 



where variations with respect to the symmetric (Lagrange multiplier) tensor ipy imposes the con- 
straints 

fi« = B l A B 1 - -S ij B k A B k = 0. (10) 
3 

When B is non degenerate the constraints are satisfied if and only if B l — ±(e° A e l + Je* • fc e' A e k ) 
which reduces the previous action to that of self-dual general relativity. Reisenberger studied 
the simplicial discretization of this action in |306] as a preliminary step toward the definition of 
the corresponding spin foam model. The consistency of the simplicial action is argued by showing 
that the simplicial theory converges to the continuum formulation when the triangulation is refined: 
both the action and its variations (equations of motion) converge to those of the continuum theory. 

In reference [505 . Reisenberger constructs a spin foam model for this simplicial theory by 
imposing the constraints fT- 7 directly on the SU(2) BF amplitudes. The spin foam path integral 
for BF theory is obtained as in Section [TU The constraints are imposed by promoting the B % to 
operators x l (combinations of left/right invariant vector fields) acting on the discrete connection^]. 
The model is defined as 

Z GR = f Hdg e 6(W) Yl II A, \v Jf [<j; ...^K- (12) 

fV[A] fV[B] e ./Tr[BAF(A)] 

where Cl = J 1 A J 3 — i<5 IJ J k A Jk and we have indicated the correspondence of the different 
terms with the continuum formulation. The preceding equation is rather formal; for the rigorous 
implementation see |305] . Reisenberger uses locality so that constraints are implemented on a 
single 4-simplcx amplitude. There is however a difficulty with the this procedure: the algebra of 
operators Cl^ do not close so that imposing the constraints sharply becomes a too strong condition 
on the BF configurations!!. In order to avoid this, Reisenberger defines a one-parameter family of 
models by inserting the operator 

e-i^ 2 (13) 

instead of the delta function above. In the limit z — > oo the constraints are sharply imposed. This 
introduces an extra parameter to the model. The properties of the kernel of f2 have not been 
studied in detail. 



7 Notice that (for example) the right invariant vector field J l {XJ) = tr 1 j^U B c d / dU^ has a well defined action 
at the level of Equation H245H and acts as a B operator at the level of ( 12421 ) since 

-i X l (U) [e™^] =Tv[^UB]^ BU \,^ ~B*^ BU \ (11) 

where cr l are Pauli matrices. 

8 This difficulty also arises in the Barrett— Crane model as we shall see in Section [TOl 
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2.3.2 The Freidel Krasnov prescription 



Freidel and Krasnov |175] define a general framework to construct spin foam models corresponding 
to theories whose action has the general form 

S(B, A) = J Tr [B A F{A)] + ®{B), (14) 

where the first term is the BF action while $(-B) is a certain polynomial function of the B field. 
The formulation is constructed for compact internal groups. The definition is based on the formal 
equation 

/ V[B]V[A] e if^AF(A)]+^B) ._ e ifH£) Z [J] _ , (15) 
J J — o 

where the generating functional Z[J] is defined as 

Z[J] := [ V[B]D[A] e */Tr[SAF(A)] + Tr[BAJ] ; 



where J is an algebra valued 2-form field. They provide a rigorous definition of the generating 
functional by introducing a discretization of A4 in the same spirit of the other spin foam models 
discussed here. Their formulation can be used to describe various theories of interest such as BF 
theories with cosmological terms, Yang-Mills theories (in 2 dimensions) and Riemannian self-dual 
gravity. In the case of self dual gravity B and A are valued in sw(2), while 



$(B) = J ^ 



B i A B j - -S lj B k A B k 
3 



(17) 



according to Equation @. The model obtained in this way is very similar to Reisenberger's 
one. There are however some technical differences. One of the more obvious one is that the 
non-commutative invariant vector fields J 1 representing B l are replaced here by the commutative 
functional derivatives S/SJ 1 . 

The idea of using such generating functional techniques have regain interest in the context 
of the so called McDowell-Mansouri formulation of general relativity [253] , The interest of the 
latter formulation is that it provides an action of gravity that is given by a BF term plus a 
genuine potential term instead of constraints. The implementation of the spin foam quantization 
of such formulation is investigated in |3431 13401 1183] , see also |325j for an important remark on the 
approach. More recently similar techniques have been used in [267] , 



2.3.3 The Iwasaki model 

Iwasaki defines a spin foam model of self dual Riemannian gravitj^ by a direct lattice discretization 
of the continuous Ashtekar formulation of general relativity. The action is 



S(e,A) = J dx 



4 ^jiv\a 



O p , j k 



2 d [x A\ 



A[A r : 



Irn X 



(18) 



where A l a is an SU(2) connection. The fundamental observation of 224 is that one can write 
the discrete action in a very compact form if we encode part of the degrees of freedom of the 
tetrad in an SU{2) group element. More precisely, if we take t/ MM = e^e^j = 1 we can define 



9 Iwasaki defines another model involving multiple cellular complexes to provide a simpler representation of 
wedge products in the continuum action. A more detail presentation of this model would require the introduction 



of various technicalities at this stage so we refer the reader to 225 . 



13 



e M := e^ + idie 1 ^ £ SU(2) where at are the Pauli matrices. In this parameterization of the 'angular' 
components of the tetrad and using a hypercubic lattice the discrete action becomes 

Sa = -p £ Tr [e], e„U Xa \ , (19) 

veA 

where r M := ((3 1 ^ 2 £ p )~ 1 e^/gJ^, is the holonomy around the ^-plaquette, e the lattice constant 
and (3 is a cutoff for r p used as a regulator (r M < /3 1 / 2 ^ p e~ 1 ). The lattice path integral is defined 
by using the Haar measure both for the connection and the 'spherical' part of the tetrad e's and 
the radial part <ir M := dr^r^. The key formula to obtain an expression involving spin foams is 

e i^[ U]= J2(2 j + l)^±lM XjiU) . (20 ) 
j 

Iwasaki writes down an expression for the spin foam amplitudes in which the integration over the 
connection and the e's can be computed explicitly. Unfortunately, the integration over the radial 
variables r involves products of Bessel functions and its behavior is not analyzed in detail. In 3 
dimensions the radial integration can be done and the corresponding amplitudes coincide with the 
results of Section 114.51 

2.3.4 The Barrett Crane model 

The appealing feature of the previous models is the clear connection to loop quantum gravity, since 
they are defined directly using the self dual formulation of gravity (boundary states are SU(2)- 
spin networks). The drawback is the lack of closed simple expressions for the amplitudes which 
complicates their analysis. There is however a simple model that can be obtained as a systematic 
quantization of simplicial 5*0(4) Plebanski's action. This model was introduced by Barrett and 
Crane in |53] and further motivated by Baez in |34j . The basic idea behind the definition was that 
of the quantum tetrahedron introduced by Barbieri in [52 and generalized to 4d in [36]. The beauty 
of the model resides in its remarkable simplicity. This has stimulated a great deal of explorations 
and produced many interesting results. 

2.3.5 Markopoulou Smolin causal spin networks 

Using the kinematical setting of LQG with the assumption of the existence of a micro-local (in the 
sense of Planck scale) causal structure Markopoulou and Smolin define a general class of (causal) 
spin foam models for gravity [261, 258] (see also |262j ). The elementary transition amplitude 
A Sl ^s I+1 from an initial spin network sj to another spin network sj+i is defined by a set of simple 
combinatorial rules based on a definition of causal propagation of the information at nodes. The 
rules and amplitudes have to satisfy certain causal restrictions (motivated by the standard concepts 
in classical Lorentzian physics). These rules generate surface-like excitations of the same kind we 
encounter in the more standard spin foam model but endow the foam with a notion of causality. 
Spin foams J-^.^. Sf are labeled by the number of times these elementary transitions take place. 
Transition amplitudes are defined as 

JV-l 

{si,s f )=^2l[A SI ^ SI+1 . (21) 

N 1=0 

The models are not related to any continuum action. The only guiding principles are the restrictions 
imposed by causality, simplicity and the requirement of the existence of a non-trivial critical 
behavior that would reproduce general relativity at large scales. Some indirect evidence of a 
possible non trivial continuum limit has been obtained in some versions of the model in 1+1 
dimensions. 
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2.3.6 Gambini Pullin model 

Gambini and Pullin [186] introduced a very simple model obtained by modification of the BF 
theory skein relations. As we argued in Section [14] skein relations defining the physical Hilbert 
space of BF theory are implied by the spin foam transition amplitudes. These relations reduce the 
large kinematical Hilbert space of BF theory (analogous to that of quantum gravity) to a physical 
Hilbert space corresponding to the quantization of a finite number of degrees of freedom. Gambini 
and Pullin define a model by modifying these amplitudes so that some of the skein relations are now 
forbidden. This simple modification frees local excitations of a field theory. A remarkable feature 
is that the corresponding physical states are (in a certain sense) solutions to various regularizations 
of the scalar constraint for (Riemannian) LQG. The fact that physical states of BF theory solve the 
scalar constraint is well known [347] . since roughly F(A) = implies EEF(A) — 0. The situation 
here is of a similar nature, and - as the authors argue - one should interpret this result as an 
indication that some 'degenerate' sector of quantum gravity might be represented by this model. 
The definition of this spin foam model is not explicit since the theory is directly defined by the 
physical skein relations. 

2.3.7 Capovilla Dell Jacobson theory on the lattice 

The main technical difficulty that we gain in going from 3-dimensional general relativity to the 
4-dimensional one is that the integration over the e's becomes intricate. In the Capovilla-Dell- 
Jacobson [1041 1103] formulation of general relativity this 'integration' is partially performed at the 
continuum level. The action is 



where e • F A F := e Fa^Fd. Integration over r] can be formally performed in the path integral 
and we obtain 



One serious problem of this formulation is that it corresponds to a sector of gravity where the 
Weyl tensor satisfy certain algebraic requirements. In particular flat geometries are not contained 
in this sector. 

2.3.8 The Engle Pereira-Rovelli Livine (EPRL) 

A modification of the Barrett-Crane model was recently introduced in [162, 163 and extended for 
arbitrary Immirzi parameter in [158] . The basic idea was to relax the imposition of the Plebanski 
constraints that reduce BF theory to general relativity in the quantum theory. The anomalous 
commutation relations of the £rn field in the quantum theory imply that the commutation of the 
Plebanski constraints does not define a closed algebra. Imposing the constraints strongly as in 
the Barrett-Crane model implies the imposition of additional conditions that are not present in 
the classical theory. There is a natural way to relax the constraints and this leads to a simple 
model that has more clear relationship with the canonical picture provided by LQG. The detailed 
description of this model will be the main subject of the following discussion of spin foam models 
for general relativity in 4 dimensions. 

10 The six components B IJ are associated to invariant vector fields in the Lorentz group as a direct consequence of 
the discretization procedure in spin foams. This makes them non commutative and render the simplicity constraints 
(functionals of the B field alone) non commutative. The origin of the non commutativity of B fields is similar to 
the non commutativity of fluxed in canonical LQG |26j . 




(22) 




X 



(23) 
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2.3.9 Freidel Krasnov (FK) 



A very similar set of models were independently introduced by Freidel and Krasnov in [176] , 
Indeed these models are arguably the same as the EPRL model for a suitable range of the Immirzi 
parameter (the vertex amplitudes coincide for 7 < 1 as we will discuss in the sequel). However, 
the logic used in the derivation is different. The idea is to express the BF path integral in terms 
of the coherent state representation in order to impose the Plebanski constraints semiclassically 
in terms of expectation values. The coherent intertwiner basis relevant for the definition of spin 
foams and the 4d quantum gravity amplitude was also derived by Livine and Speziale in 246 . 
Freidel and Krasnov introduced the linear version of Plebanski constraints that where later used in 
order to provide a simpler derivation and generalization of the EPRL model to arbitrary Immirzi 
parameter. 

The last two models are going to be the subject of most of our discussions in what follows. 
They are the most developed candidates in four dimensions and we will review their properties in 
detail in the following part. 
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Part II 

Preliminaries: LQG and the canonical 
quantization of 4d gravity 

In this part we briefly review the basic results of the loop quantum gravity approach to the 
canonical quantization of gravity. This part is relevant for the interpretation of the new spin foams 
models presented in the next part. The reader interested in the canonical formulation of general 
relativity in terms of connection variables is referred to the text books [3151 13531 [20] and the review 
article [31]. For a pedagogical introduction see [296]. 

In the following section we review the canonical analysis of general relativity formulated as 
a constrained BF theory, i.e., in the so-called Plebanski formulation [301] . The study is done in 
terms of the variables that, on the one hand, allow for the closest comparison with the discrete 
variables used in the spin foam approach, while, on the other hand, allow for the introduction the 
basic elements of the canonical quantization program of loop quantum gravity. The latter being, 
in my view, the clearest possible setting in which spin foams are to be interpreted. 

In Section H] we present the basic elements of loop quantum gravity (LQG). In Section [5] we 
give a short description of the spin foam representation of the path integral of gravity. Concrete 
realizations of this will be reviewed in 4d in Part IIIII and in 3d in Part IIVI 
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3 Classical General Relativity in Connection Variables 



The Hamiltonian analysis of general relativity is the basic starting point for canonical quantization. 
Loop quantum gravity and spin foams are based on the possibility of formulating Hamiltonian 
general relativity in terms of Yang-Mills-like connection variables. The primitive ancestor of 
these formulations is Ashtekar's self-dual complex connection formulation of general relativity [19j . 
Modern, LQG and spin foams are based on a certain relative of these variables that is often 
referred to as Ashtekar-Barbero variables [50 ] 122 1J . These variables can be introduced directly 
by means of a canonical transformation of the phase space of the Palatini formulation of general 
relativity [350] , or can be derived directly from the canonical analysis of the so called Hoist action 
for general relativity [220 . More simply the new variables also follow from the addition of the 
Nieh-Yan topological invariant to the Palatini action |125j . The Ashtekar-Barbero connection 
parametrization of the phase space of general relativity also arises naturally from the consideration 
of the most general diffcomorphism invariant action principle that can be written for the field 
content of the Palatini first order formulation |312] (i.e. a Lorentz connection oj^ j = — w„ and a 
co-tetrad with a spacetime indeces and I, J = 0, 1, 2, 3 internal Lorentz indeces). 

The new spin foam models are based on Plebanski's realization that one can obtain classical 
general relativity by suitable constraining the variations of the B field in the simple kind of topo- 
logical theory called BF theory |301] . For that reason, the best suited action principle of gravity 
for the study of the new spin foam models corresponds to one of the Plebanski type. However, the 
details of such treatment are only partially presented in the literature (see for instance [198] . [362] ). 
The material of the following subsection is meant to complete this void, the study is largely inspired 
in the perspective exploited in references |270) 1359] for the canonical formulation of constrained 
BF theory. As we will show below, the hamiltonian analysis of such formulation of gravity is at 
the heart of the replacement of the Plebanski quadratic constraints by linear ones that has been 
so useful in the definition of the new vertex amplitudes [176] . 



3.1 Gravity as constrained BF theory 

In order to keep the presentation simple (and not to bother with the ± signs appearing in dealing 
with raising and lowering of Lorentz indices) we present here the Hamiltonian formulation of 
Riemannian gravity. This is enough to illustrate the algebraic structures that are necessary to 
introduce in the present contexlFI. 

Our starting point is Plebanski's Riemannian action for general relativity which can be thought 
of as Spin(4) BF theory plus constraints on the B-field. It depends on an so(4) connection w, a Lie- 
algebra-valued 2-form B and Lagrange multiplier fields A and fi. Writing explicitly the Lie-algebra 
indexes, the action is given by 



S[B,lj,X] 



(*B IJ + -B IJ ) A Fjj(u) + Xijkl B IJ A B KL 

1 



(24) 



where /i is a 4-form and Xijkl = —^jikl = —^ijlk = ^klij is a tensor in the internal space 
satisfying £ IJKL Xjjkl — 0, and 7 is the Immirzi parameter. The Lagrange multiplier tensor 
A/jifi has then 20 independent components. The previous action is closely related to the one 
introduced in [1051 13591 1271] (there is a simple analog of the previous action in 3d [90]). Variation 
with respect to A imposes 20 algebraic equations on the 36 components of B. They are 

e ^ B jJ B KL = e e IJKL ( 25 ) 

where e = ^opqrB^ B^ e^ vp(T [126] , The solutions to these equations are 

B = ±*(eAe), and B = ±eAe, (26) 



I am grateful to M. Montesinos and M. Velazquez for discussions and contributions to the present section. 
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in terms of the 16 remaining degrees of freedom of the tetrad field e„. If one substitutes the first 
solution into the original action one obtains Hoist's formulation of general relativity [220] 



S[e,u] = - I Tr 



(*e A e + — e A e) A F(A) 

7 



(27) 



This property is the key to the definition of the spin foam model for gravity of the next section. 
3.2 Canonical analysis 

Now one performs the usual 3+1 decomposition by assuming the existence of a global time function 
t whose level hypersurfaces S t define a foliation of the spacetime manifold The previous action 
takes the following form 

S\II,H,U),X] = J dt J Tr[II Aw + w D u AU + HAF + X 1 (7* if — if ) A ( 7 *n — II)], (28) 

where 

tt/J _ *p>IJ , 1 nU 
11 ab - U ab + ~ B ab 

1 

Hl J = *Btf + ±Bg, (29) 

for a, b — 1,2,3 are St-cotangent-space abstract indices, A' = 7 2 / (1 — 7 2 ) 2 A, and the trace is the 
obvious contraction of internal indices. By choosing an internal direction n 1 (and in the gauge 
where n 1 = (1,0,0,0)) we can now introduce a convenient reparametrization of the unconstrained 
phase space by writing the 18 canonical pairs (11^, w^) by (P(7; ±)* h , ± Aj) where 

nr,±y ab = \e* Jk ni k b ±±-n%, ^ = ^^±7^- (30) 

For notational convenience we will also introduce 



H(7;±)l = ^tff±i-if°\ 



We can rewrite the previous action as 

S[P( 7 ; ±), ± A, A] = J dt J P( 7 ; +)i A + A 1 + P( 7 ; -)i A ± l + N t G* + m JB ! + 
+H A F + A' ( 7 *if - H) A ( 7 *II - IL)], (31) 
where TV 1 , rf, and A are Lagrange multipliers imposing the constraints: 

B l = (£> (+A)+( - A) A P( 7 ; +Y) « 

2 

(E* = _L e ^(+A - - A), A P( 7 ; +) k « 0, 
2 7 

(¥ l 3 k HJk - h 01 ) A (i ei jfe nifc - n °0 + (i«i)«o 
iVjk Hjk ~ H ° l ) A (ie'" i n 0s - n' m )+ 

~y Jm TjOk Tjlm\ a n^i TTjfc ttO£\ r 



/j-z^ pa 



+ (ieTif ufc - if lm ) A (§e J , fc IF* - If") 
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In the previous equation we explicitly recall that the last three equations are a rewriting of the 
simplicity constraints (|25|) . We can rewrite them in terms of F(7; — ) and P(-; — ) as 

B* = £> ( - A)+(+A) A P( 7 ; +Y ps 0, 

2 

<S l = -U^(+A - "A)j A P( 7 ; +)* « 0, 
27 

P J ' =H( 7 ;-) (3 AP( 7 ;-) 4) «0, 

H« =E(-;-) (l AP(-;-) 3) w0, 

7 7 

Iir j = H( 7 ; A ¥(-; -V + H(-; -V A P( 7 ; -V 

7 7 

-i[H( 7 ; -) fc A P(-; -) fc + H(-; -) k A P( 7 ; -) fe ]£« « 0. (32) 
3 7 7 

The 20 constraints F J , IP 7 , and IIP-' are referred to as the simplicity constraints. The general 
solution of the simplicity constraints (|26[) requires, in the 3+1 setting, the introduction of four new 
fields: the lapse N and the shift vectors N a . The solution is 

P(7;-) l a6 = o, 

1 (1-r 2 ) 

P(-;-)a 6 = ^ Zi P(7;+)l h , 
7 2 



H( -.:-),', = V iL-j!l e ^(p( 7;+ )), 



1 fl-7 2 ) 

H(-;-); = ^-^V c P( 7; +)* a , 
7 2 

with 

■ir:..D) ^^ M '^ + >^ + >?. (33) 

8 ^det(P(i;+)) 

This solution breaks the internal Lorentz gauge as the first condition above does not commute 
with the boost constraint IB 1 : as it will become clear below it amounts to choosing the so-called 
time gauge. The above equations parametrize the solutions of the Plebanski constraints in terms 
of the 9 components of P(7; +)a& pl us 4 extra parameters given by a scalar TV and a space tangent 
vector N a £ T(Y,) and the 3 parameters in the choice of an internal direction n 1 (here implicitly 
taken as n 1 = (1,0,0,0)). These are exactly the 16 parameters in the co-tetrad in (|26|) . 
The components of H that are necessary to write the term H A F in the action (|40|) are 

H Q a l = 1 F( r ,+y ca N c + -e* a 
7 

\e%H^ =I>(T,+) k C aN c + Nel 

The conservation in time of the set of constraints above lead to 6 secondary constraints. As shown 
in [312J these follow from 



<C ij = e^MPfx-f 1 oc Ne u A 



d e i ) + ^L( + A+-A) i Ae" 



(34) 



Notice now that the set of 9 constraints B 1 « and (D*- 13 * 1 « are indeed equivalent to the following 
simple condition 

IV l a = i( + A + -A)l-rj [ (e)«0, (35) 
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where is the torsion free spin connection compatible with the triad e l a , i.e. the unique solution 
of Cartan's structure equations 

^4]+4 r [ a 4]=0' (36) 

Remark: The primary constraint P(7;— ) ks (first line in Equation h33\l ) will have central 
importance in the definition of the new spin foam models. Explicitly, from i30\) , one has 

This is exactly the linear simplicity constraints that we will rediscover in the simplicial setting of 
Sections [?] and and which will be imposed in a suitable sense at all tetrahedra of the cellular 
decomposition in order to satisfy the Plebanski constraints and produce a state sum model for 
gravity out of that of BF theory. The consistency requirement are expected to hold in the 
discrete setting from the fact that one is imposing the linear simplicity constraint for all tetrahedra, 
thus P( 7 ; — ) ~ is valid for all times. 



Now from the fact that 



1 



{lVl(x), P( 7 ; -y bc (y)} = -e abc 5VS(x, y) (38) 

we conclude that the given pair are second class. We can substitute their solution in the action. 
In particular, from eqs. (|30[) and (|35|) . + Jk\ becomes the Ashtekar-Barbero connection 

+K = K + iK (39) 

where we have defined K l a = (+A^ - "A^)/2. With all this 
S[P( 7 ;+), + A,A^,Ay = 
= ^JdtJ P( 7 ;+), A+A* + 

+N Z D +A A P( 7 ; +) 1 + NJP{ T , +) fe A ( 7 F 0i + F' k e ijk ) + Ne l A (— + F^ k e ijk ). (40) 

7 

The previous action corresponds to the standard Hamiltonian formulation of general relativity in 
terms of SU (2) connection variables. This becomes more transparent if we now change to a more 
standard notation by introducing the socalled densitized triad 

E? =~ f e abc V( r ,+) i bc (41) 
in terms of which the last line of Equation (|33l) becomes 



j and ■■ (42) 

From now on we denote the Ashtekar-Barbero [51] connection A l a given by 

Al = + Al^r a+1 Kl (43) 
The Poisson brackets of the new variables are 

{E«(x),Al(y)} = K1 S a b S;S(x,y), {E° (x) , E<? (y)} = {A{{x), A{{y)} = 0. (44) 
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All the previous equations follow explicitly from (|40j) except for (A J a (x), A\{y)} = which follows 
from the special property of T l a in three dimensions that 



8W[E 
~SEf 



(45) 



where W[E] is the generating functional for the spin connection^!. The action (|40|) becomes 




(46) 



where the constraints are explicitly given by: 



V b (E*,Al) = E?F ab 



(l + 1 2 )K a G t 



(47) 




for the 18 phase space variables (A a ,Ej). In addition to imposing conditions among the canonical 
variables, first class constraints are generating functionals of (infinitesimal) gauge transformations. 
From the 18-dimensional phase space of general relativity we end up with 11 fields necessary 
to coordinatize the constraint surface on which the above seven conditions hold. On that 11- 
dimensional constraint surface, the above constraint generate a seven-parameter-family of gauge 
transformations. The reduce phase space is four dimensional and therefore the resulting number 
of physical degrees of freedom is two, as expected. 

The constraint (14U)) coincides with the standard Gauss law of Yang-Mills theory (e.g. V • E = 
in electromagnetism). In fact if we ignore (l47|) and (|48|) the phase space variables (A a ,Ej) 
together with the Gauss law (|49|) characterize the physical phase space of an SU(2j3 Yang-Mills 
(YM) theory. The gauge field is given by the connection A a and its conjugate momentum is the 
electric field Ej. Yang-Mills theory is a theory defined on a background spacetime geometry. 
Dynamics in such a theory is described by a non vanishing Hamiltonian - the Hamiltonian density 
of YM theory being H = E a Ef + B a Bf. General relativity is a generally covariant theory and 
coordinate time plays no physical role. The Hamiltonian is a linear combination of constraints!^! 
Dynamics is encoded in the constraint equations (14"T|) . (|48I) . and (|4"9"|) . In this sense we can regard 
general relativity in the new variables as a background independent relative of SU(2) Yang-Mills 
theory. We will see in the sequel that the close similarity between these theories will allow for the 
implementation of techniques that are very natural in the context of YM theory. 

12 In fact equation l|36|l implies that rj,[_E] = r^[Ai?] for a constant A. This homogeneity property plus the 
existence of W^-E] imply that W[£] = f E^F' a as it can be verified by direct calculation. 

13 The constraint structure does not distinguish 50(3) from 5(7(2) as both groups have the same Lie algebra. 
From now on we choose to work with the more fundamental (universal covering) group SU(2). In fact this choice 
is physically motivated as 51/(2) is the gauge group if we want to include fermionic matter [272] . 

14 In the physics of the standard model we are used to identifying the coordinate t with the physical time of a 
suitable family of observers. In the general covariant context of gravitational physics the coordinate time t plays the 
role of a label with no physical relevance. One can arbitrarily change the way we coordinatize spacetime without 
affecting the physics. This redundancy in the description of the physics (gauge symmetry) induces the appearance 
of constraints in the canonical formulation. The constraints in turn are the generating functions of these gauge 
symmetries. The Hamiltonian generates evolution in coordinate time t but because redefinition of t is pure gauge, 
the Hamiltonian is a constraint itself, i.e. H = on -lull l.'i",' 217| . More on this in the next section. 
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To conclude this section let us point out that the real connection formulation of general relativity 
presented here is a peculiar property of four dimensions due to the special property (|4"5"]) . Never- 
theless, there are means to obtaining real connection formulations for gravity and supergravity in 
higher dimensions as recently shown by Thiemann and collaborators [75] [76] [77] [78] [79 81, 82] [80]. 

3.2.1 Constraints algebra 

Here we simply present the structure of the constraint algebra of general relativity in the new 
variables. 

{G(a),G(j3)} = G([a,j3]), (50) 
where a — a l Ti £ su(2), f3 = f3 l Ti € su(2) and [a, /3] is the commutator in su(2). 

{G(a),V(N a )} = -G(J? N a). (51) 

{G(a),S(N)} = 0. (52) 

{V{N a ), V(M a )} = V([N, M] a ), (53) 

where [N, M] a = N b dbM a — M b dbN a is the vector field commutator, and Jzf/v denotes the Lie 
derivative along the vector field iV a . The previous constraints define the subalgebra od spacial dif- 
feomorphisms and SU(2) internal gauge transformations. This property allows to implement them 
in the quantum theory separately from the scalar constraint. If we include the scalar constraint 
the remaining Poisson brackets are 

{S(N), V(N a )} = -S(J? N N), (54) 

and finally 

{S(N),S(M)} = V(S a ) + terms proportional to the Gauss constraint, (55) 

where for simplicity we are ignoring the terms proportional to the Gauss law (the complete ex- 
pression can be found in |31j ) and 

j^a j^b £ij 

S a = 'J (Nd b M ~ Md b N). (56) 



Notice that instead of structure constants, the r.h.s. of (|55[) is written in terms of field dependent 
structure functions. For this reason it is said that the constraint algebra does not close in the BRS 
sense. 



3.3 Geometric interpretation of the new variables 

The geometric interpretation of the connection A l a , defined in (I43[) , is standard. The connection 
provides a definition of parallel transport of SU(2) spinors on the space manifold E. The natural 
object is the SU (2) element defining parallel transport along a path e C S also called holonomy 
denoted h e [A\ , or more explicitly 



h e [A] = Pexp- 




e 



where P denotes a path-order-exponential. 

The densitized triad - or electric field - Ef also has a simple geometrical meaning. As it follows 
from (|42[) . Ef encodes the full background independent Riemannian geometry of S. Therefore, any 
geometrical quantity in space can be written as a functional of Ef. One of the simplest is the area 
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As[Ef] of a surface S C E whose expression we derive in what follows. Given a two dimensional 
surface in S C £ - with normal 

dx b dx c 

Ttabc (58) 



where er 1 and cr 2 are local coordinates on 5 - its area is given by 



Ml ab } = [ Vhd^da 2 , 
Js 



(59) 



where h = det(h ab ) is the determinant of the metric h ab — q ab — n 2 n a rib induced on S by q ab . 
From equation gSJ) it follows that qq ab = EfE bl so that det(q ab ) = det(£f). Contracting the 
previous equality with n a rib, namely 

qq ab n a n b = E^E b 5^n a n b . (60) 

Now observe that q nn = q ab n a n b is the nn-matrix element of the inverse of q ab . Through the well 
known formula for components of the inverse matrix we have that 

„„ = det(g afc ~ n- 2 n a n b ) = h 
dct(q a6 ) q' 

But q a i, — n n a rib is precisely the induced metric h a b- Replacing q nn back into (|B(7)) we conclude 
that 

h = E'; E':d' >n„n,.. (62) 
Finally we can write the area of S as an explicit functional of Ef : 

MEt\ = J ^EfE b 5^n a n b da 1 da 2 . (63) 

This simple expression for the area of a surface will be very important in the quantum theory. 



4 Loop Quantum Gravity and Quantum Geometry in a Nut- 
shell 

Loop quantum gravity is a proposal for the implementation of the quantization program established 
in the 1960s by Dirac, Wheeler, and DeWitt, among others (for recent reviews see [3T1 13501 1317] ). 
The technical difficulties of Wheeler's 'gcomctrodynamics' are circumvent by the use of connection 
variables instead of metrics [20l HU [50] . At the kinematical level, the formulation is similar to 
that of standard gauge theories. The fundamental difference is however the absence of any non- 
dynamical background field in the theory. 

The configuration variable is an S'C/(2)-connection A l a on a 3-manifold S representing space. 
The canonical momenta are given by the densitized triad The latter encode the (fully dynam- 
ical) Riemannian geometry of E and are the analog of the 'electric fields' of Yang-Mills theory. 

In addition to diffeomorphisms there is the local SU(2) gauge freedom that rotates the triad 
and transforms the connection in the usual way. According to Dirac, gauge freedoms result in 
constraints among the phase space variables which conversely are the generating functionals of 
infinitesimal gauge transformations. In terms of connection variables the constraints are 

g t =V a Ef = 0, C a = E b F b k a = 0, S .■ a ir;i: b l-„,,. .... 0. (64) 

where T> a is the covariant derivative and F ba is the curvature of A l a . Qi is the familiar Gauss 
constraint - analogous to the Gauss law of electromagnetism - generating infinitesimal SU(2) 
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gauge transformations, C a is the vector constraint generating spacc-diffeoniorphism, and S is the 
scalar constraint generating 'time' reparameterization (there is an additional term that we have 
omitted here for simplicity - see Equations (|47l) to (|49l) , and Section 13.2.11 for the precise form of 
the constraints and their relationship with the gauge symmetry groups). 

Loop quantum gravity is defined using Dirac quantization. One first represents (|64|) as operators 
in an auxiliary Hilbert space % and then solves the constraint equations 

G^ = 0, C„tf = 0, S* = 0. (65) 

The Hilbert space of solutions is the so-called physical Hilbert space H p h ys ■ In a generally covariant 
system quantum dynamics is fully governed by constraint equations. In the case of loop quantum 
gravity they represent quantum Einstein's equations. 

States in the auxiliary Hilbert space are represented by wave functionals of the connection 
^f(A) which are square integrable with respect to a natural diffcomorphism invariant measure, the 
Ashtekar-Lewandowski measure [2"51 (we denote it Jz? 2 [A] where A is the space of (generalized) 
connections). Roughly speaking l 15 l. this space can be decomposed into a direct sum of orthogonal 
subspaces % = @ 7 "H 7 labeled by a graph 7 in S. The fundamental excitations are given by the 
holonomy ht [ A) £ SU(2) along a path £ in E: 

ht(A) =V exp / A. (66) 
Jt 

Elements of % 7 are given by functions 

* f , 7 (A) = f(h il (A),he 2 (A),...,h en (A)), (67) 

where hi is the holonomy along the links £ € 7 and / : SU(2) n — > C is (Haar measure) square 
integrable. They are called cylindrical functions and represent a dense set in H denoted Cyl. 
The "momenta" conjugate to the holonomies are given by the so-called flux operators 



E(S, a) = [ a l Ei (68) 
Js 



acros a 2-suface s C £ and labelled by a smearing field a € su(2). It has been shown that the 
associated (Poisson) holonomy-flux algebra admits a unique quantization in a Hilbert space with a 
diffeomorphism invariant states 1170) 1239] . The Hilbert space H mentioned above is precisely that 
unique representation. 

Gauge transformations generated by the Gauss constraint act non-trivially at the cndpoints 
of the holonomy, i.e., at nodes of graphs. The Gauss constraint (in (15^]) ) is solved by looking at 
SU(2) gauge invariant functionals of the connection (^f 2 [A]/G). The fundamental gauge invariant 
quantity is given by the holonomy around closed loops. An orthonormal basis of the kernel of 
the Gauss constraint is defined by the so called spin network states 1 ^ 1 ,{j e }.{, n }(A) [304, 323, 33 . 
Spin-network^^ are defined by a graph 7 in S, a collection of spins {je} - unitary irreducible 
representations of SU(2) - associated with links I G 7 and a collection of SU(2) intertwiners 
{(.„} associated to nodes n € 7 (see Figure [3]). The spin- network gauge invariant wave functional 
i & J ,{3,},{i. n }(A) is constructed by first associating an SU(2) matrix in the ^-representation to the 



15 The construction is rather a projective limit |28| . This leads to non-trivial requirements (cylindrical consistency) 
on the structure (operators) of the theory. This point is very important for spin foams specially in considering the 
refining limit of amplitudes. We will revisit this issue in Section 17.51 

16 Spin-networks were introduced by Penrose [291 , 292 290 , 293,, in a attempt to define 3-dimensional Euclidean 
quantum geometry from the combinatorics of angular momentum in QM. Independently they have been used in 
lattice gauge theory as a natural basis for gauge invariant functions on the lattice. For an account of their applications 
in various contexts see 13381. 
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holonomies he(A) corresponding to the link £, and then contracting the representation matrices at 
nodes with the corresponding intertwiners i n , namely 

%.W.i*y(A) = II II 3iMA% (69) 

where ji[h e (A)] denotes the corresponding ^-representation matrix evaluated at corresponding link 
holonomy and the matrix index contraction is left implicit. 



1 

2 




Figure 3: Spin-network state: At 3-valent nodes the intertwiner is uniquely specified by the cor- 
responding spins. At 4 or higher valent nodes an intertwiner has to be specified. Choosing an 
intertwiner corresponds to decompose the n-valent node in terms of 3-valent ones adding new vir- 
tual links (dashed lines) and their corresponding spins. This is illustrated explicitly in the figure 
for the two 4-valent nodes. 

The solution of the vector constraint is more subtle [35] . One uses group averaging techniques 
together with the diffeomorphism invariance of the kinematical inner product in H. The diffeo- 
morphism constraint does not exist in the quantum theory. Only finite diffeomorphisms can be 
defined. As a result solutions (diffeomorphism invariant states) correspond to generalized states. 
These are not in H but are elements of the topological dual Cyi*0- However, the intuitive idea 
is quite simple: solutions to the vector constraint are given by equivalence classes of spin-network 
states up to diffeomorphism. Two spin-network states are considered equivalent if their underlying 
graphs can be deformed into each other by the action of a diffeomorphism. 

This can be regarded as an indication that the smooth spin-network category could be replaced 
by something which is more combinatorial in nature so that diffeomorphism invariance becomes 
a derived property of the classical limit. LQG has been modified along these lines by replacing 
the smooth manifold structure of the standard theory by the weaker concept of piecewise linear 
manifold [364 . In this context, graphs defining spin- network states can be completely character- 
ized using the combinatorics of cellular decompositions of space. Only a discrete analog of the 
diffeomorphism symmetry survives which can be dealt with in a fully combinatorial manner. We 
will take this point of view when we introduce the notion of spin foam in the following section. 

4.1 Quantum geometry 

The generalized states described above solve all of the constraints (f6"4"|) but the scalar constraint. 
They are regarded as quantum states of the Riemannian geometry on S. They define the kine- 
matical sector of the theory known as quantum geometry. 

Geometric operators acting on spin network states can be defined in terms of the fundamental 

17 According to the triple Cyl CUC Cyl*. 
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triad operators Ef . The simplest of such operators is the area of a surface S classically given by 



A S {E) = [ dx 2 ^Tr[n a n b E»Et>] (70) 
«/ s 

where n is a co normal. The geometric operator A${E) can be rigorously defined by its action on 
spin network states |337l 13221 [29] . The area operator gives a clear geometrical interpretation to 
spin-network states: the fundamental 1-dimensional excitations defining a spin-network state can 
be thought of as quantized 'flux lines' of area. More precisely, if the surface S C £ is punctured 
by a spin- network link carrying a spin j, this state is an eigenstate of A,s{E) with eigenvalue 
proportional to i 2 P ^ ' j(j + 1). In the generic sector - where no node lies on the surface - the 
spectrum takes the simple form 

as({j}) = 8*iP P £ Vj^h+T), (71) 

i 

where i labels punctures and i is the Imirzi parameter |22lf^l . as{{j}) is the sum of single puncture 
contributions. The general form of the spectrum including the cases where nodes lie on S has been 
computed in closed form [29] , 

The spectrum of the volume operator is also discrete [3371 13221 12521 150] . If we define the 
volume operator V a (E) of a 3-dimensional region a C X then non vanishing eigenstates are given 
by spin-networks containing n-valent nodes in a for n > 3. Volume is concentrated in nodes. For 
new results on the volume see |661 165] . Other geometric quantities have been considered in the 
quantization; in the particular in studies of coupling LQG to matter [348 the introduction of a 
metric operator is necessary. For another proposal of length operator see |64] . 



4.2 Quantum dynamics 

In contrast to the Gauss and vector constraints, the scalar constraint does not have a simple 
geometrical meaning. This makes its quantization more involved. Regularization choices have to 
be made and the result is not unique. After Thiemann's first rigorous quantization [346] other well 
defined possibilities have been found |195) I129| 1130] . This ambiguity affects dynamics governed by 

5* = 0. (72) 

The difficulty in dealing with the scalar constraint is not surprising. The vector constraint - 
generating space diffeomorphisms - and the scalar constraint - generating time reparameterizations 
- arise from the underlying 4-diffeomorphism invariance of gravity. In the canonical formulation the 
3+1 splitting breaks the manifest 4-dimensional symmetry. The price paid is the complexity of the 
time re-parameterization constraint S. The situation is somewhat reminiscent of that in standard 
quantum field theory where manifest Lorentz invariance is lost in the Hamiltonian formulation^!. 

From this perspective, there has been growing interest in approaching the problem of dynamics 
by defining a covariant formulation of quantum gravity. The idea is that (as in the QFT case) 
one can keep manifest 4-dimensional covariance in the path integral formulation. The spin foam 
approach is an attempt to define the path integral quantization of gravity using what we have 
learned from LQG. 

18 The Imirzi parameter i is a free parameter in the theory. This ambiguity is purely quantum mechanical (it 
disappears in the classical limit). It has to be fixed in terms of physical predictions. The computation of black hole 
entropy in LQG fixes the value of t (see 1221 ). 

19 There is however an additional complication here: the canonical constraint algebra does not reproduce the 
4-diffcomorphism Lie algebra. This complicates the geometrical meaning of S. 
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In standard quantum mechanics path integrals provide the solution of dynamics as a device 
to compute the time evolution operator. Similarly, in the generally covariant context it provides 
a tool to find solutions to the constraint equations (this has been emphasized formally in various 
places: in the case of gravity see for example [206] . for a detailed discussion of this in the context 
of quantum mechanics see [310 J. Recall discussion of Section [2T2T 

Let us finish by stating some properties of S that do not depend on the ambiguities mentioned 
above. One is the discovery that smooth loop states naturally solve the scalar constraint opera- 
tor [228, 339 . This set of states is clearly to small to represent the physical Hilbert space (e.g., 
they span a zero volume sector). However, this implies that S acts only on spin network nodes. 
Its action modifies spin networks at nodes by creating new links according to Figure |jp°l . This is 
crucial in the construction of the spin foam approach of the next section. 

k k 




i i 



Figure 4: A typical transition generated by the action of the scalar constraint 



5 Spin Foams and the Path Integral for Gravity in a Nut- 
shell 

The possibility of defining quantum gravity using Feynman's path integral approach has been 
considered since Misner [268] , Given a 4- manifold Ai with boundaries Si and £2, and denoting 
by G the space of metrics on M, the transition amplitude between |[<7 a f>]) on Si and ![(?„(,]) on £2 
is formally 

([q ab ]\Wa b }) = Jv[g] e iS ^\ (73) 
[9] 

where the integration on the right is performed over all space-time metrics up to 4-diffcomorphisms 
[g] £ Gj 'Diff(M) with fixed boundary values up to 3-diffeomorphisms [q a b], W a b\i respectively. 

There are various difficulties associated with ([73]) . Technically there is the problem of defining 
the functional integration over [g] on the RHS. This is partially because of the difficulties in defining 
infinite dimensional functional integration beyond the perturbative framework. In addition, there 
is the issue of having to deal with the space G/oiff(M), i.e., how to characterize the diffeomorphism 
invariant information in the metric. This gauge problem (3-diffeomorphisms) is also present in the 
definition of the boundary data. There is no well defined notion of kinematical state |[<?ab]) as the 
notion of kinematical Hilbert space in standard metric variables has never been defined. 

The situation is different in the framework of loop quantum gravity. The notion of quantum 
state of 3-geometry is rigorously defined in terms of spin- network states. They carry the diff- 
invariant information of the Riemannian structure of £. In addition, and very importantly, these 
states are intrinsically discrete (colored graphs on E) suggesting a possible solution to the functional 
measure problem, i.e., the possibility of constructing a notion of Feynman 'path integral' in a 

20 This is not the case in all the available definitions of the scalar constraints as for example the one defined 
in (129111301 . 
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combinatorial manner involving sums over spin network world sheets amplitudes. Heuristically, '4- 
geometries' are to be represented by 'histories' of quantum states of 3-geometries or spin network 
states. These 'histories' involve a series of transitions between spin network states (Figure [5]), 
and define a foam-like structure (a '2-graph' or 2-complex) whose components inherit the spin 
representations from the underlying spin networks. These spin network world sheets are the so- 
called spin foams. 

The precise definition of spin foams was introduced by Baez in |34j emphasizing their role as 
morphisms in the category defined by spin network^!. A spin foam T : s — > s' , representing 
a transition from the spin-network s = (7, {je}, {in}) into s' = (7', {jc}, {>«'})) is defined by a 
2-complex J bordered by the graphs of 7 and 7' respectively, a collection of spins {j/} associated 
with faces / G J and a collection of intertwiners {t e } associated to edges e € 3 ■ Both spins and 
intertwiners of exterior faces and edges match the boundary values defined by the spin networks s 
and s' respectively. Spin foams J- : s — > s' and T' : s' s" can be composed into J-J-' : s — > s" 
by gluing together the two corresponding 2-complexes at s' . A spin foam model is an assignment 
of amplitudes A [J-] which is consistent with this composition rule in the sense that 

A[FF'] = A[F]A[F']. (74) 

Transition amplitudes between spin network states are defined by 

(*>*'W = E A ^> ( 75 ) 

J-;s— >s' 

where the notation anticipates the interpretation of such amplitudes as defining the physical scalar 
product. The domain of the previous sum is left unspecified at this stage. We shall discuss this 
question further in Section [Vj This last equation is the spin foam counterpart of Equation (1731) . 
This definition remains formal until we specify what the set of allowed spin foams in the sum are 
and define the corresponding amplitudes. 




Figure 5: A typical path in a path integral version of loop quantum gravity is given by a series of 
transitions through different spin-network states representing a state of 3-geometries. Nodes and 
links in the spin network evolve into 1-dimensional edges and faces. New links are created and 
spins are reassigned at vertexes (emphasized on the right). The 'topological' structure is provided 
by the underlying 2-complex while the geometric degrees of freedom are encoded in the labeling of 
its elements with irreducible representations and intertwiners. 

The background-independent character of spin foams is manifest. The 2-complex can be 
thought of as representing 'space-time' while the boundary graphs as representing 'space'. They do 

21 The role of category theory for quantum gravity had been emphasized by Crane in I122lll2llll24 |. 



29 



not carry any geometrical information in contrast with the standard concept of a lattice. Geometry 
is encoded in the spin labelings which represent the degrees of freedom of the gravitational field. 

In standard quantum mechanics the path integral is used to compute the matrix elements of the 
evolution operator U(t). It provides in this way the solution for dynamics since for any kinemat- 
ical state \& the state U(t)^> is a solution to Schrodinger's equation. Analogously, in a generally 
covariant theory the path integral provides a device for constructing solutions to the quantum 
constraints. Transition amplitudes represent the matrix elements of the so-called generalized 'pro- 
jection' operator P (i.e., (s,s') h — (sP,s') recall the general discussion of Sections l2.2p such 
that P^ is a physical state for any kinematical state ^. As in the case of the vector constraint 
the solutions of the scalar constraint correspond to distributional states (zero is in the continuum 
part of its spectrum). Therefore, Hphys is not a proper subspace of % and the operator P is not a 
projector (P 2 is ill definedf^l. In Section [14] we give an explicit example of this construction. 



In the notation of the previous section states in Hphys are elements of Cyl* 
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Part III 

The new spin foam models for four 
dimensional gravity 

In this part we study the non-perturbative quantization of four dimensional general relativity from 
the spin foam perspective. Together with Part [HI his part of the article form a self contained body, 
thus can be studied completely independently from the rest of the content of this review. 

The new spin foam models for four dimensional quantum gravity are introduced from a perspec- 
tive that, in some aspects, is independent from the one taken in the original works related to the so 
called EPRL model [1581 11631 1162] as well as the one used in the definition of the FK models [176] . 
Our starting point is the quantization of BF theory in the path integral formulation, which leads to 
a well defined unambiguous state sum model or topological field theory (see Section [5]). The space 
of histories of BF theory path integral will be constrained to satisfy the so-called linear simplicity 
constraints that reduce BF theory to gravity. 

In Section [5] we review the quantization of BF theory and its spin foam representation. In 
Section [7] we present the Riemannian EPRL model together with various mathematical tools that 
will be useful in the following sections. In Section [5] the Lorentzian version of the EPRL model 
is reviewed. In Section [9] we present the FK model. In Section [10] we recall the definition of the 
Barrett-Crane model. 

Finally, there are various review works in the literature [320, 319 where the new models are 
introduced from a rather minimalistic perspective by simply postulating the defining amplitudes 
and deriving their relation to gravity a posteriori. Such choice is indeed quite advantageous if 
the novel simplicity of the new models is to be emphasised. In this article we have chosen a 
complementary more constructive view. This strategy allows for a systematic presentation of the 
ingredients that go into the construction of the new models from the perspective of continuum 
general relativity and BF theory. A possible drawback of this choice is that most of the simplicial 
geometry intuitions used in other derivations are almost completely avoided. We hope that this 
lost will be compensated by the potential advantages of an alternative viewpoint. 
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6 Spinfoam Quantization of BF Theory 



Here we follow the perspective of [35]. Let G be a compact group whose a Lie algebra g has an 
invariant inner product here denoted (), and M. a d-dimensional manifold. Classical BF theory is 
defined by the action 

S\B,u] = J (BAFH), (76) 
M 

where B is a g valued (d — 2)-form, w is a connection on a G principal bundle over M.. The 
theory has no local excitations: all solutions of the equations of motion are locally related by 
gauge transformations. More precisely, the gauge symmetries of the action are the local G gauge 
transformations 

SB = [B, a], Su = d w a, (77) 

where a is a g-valued 0-form, and the 'topological' gauge transformation 

SB = d u r], Slu = 0, (78) 

where d^ denotes the covariant exterior derivative and rj is a g- valued 0-form. The first invariance 
is manifest from the form of the action, while the second is a consequence of the Bianchi identity, 
d w F(w) = 0. The gauge symmetries are so vast that all the solutions to the equations of motion 
are locally pure gauge. The theory has only global or topological degrees of freedom. 

Remark: In the special case G = SU(2) and d = 3 BF theory is (Riemannian) general relativ- 
ity where the field B l a is given by the cotretrad frames e l a of general relativity in the first order 
formalism. This simple example will be studied in more detail in Section Another case of 
interest is G = Spin(4) and d — 4 as it will provide the basis for the construction of the spinfoam 
models for four dimensional quantum (Riemannian) general relativity studied in Sections^ Q3 
and \1(A The relationship with general relativity stems from the fact that constraining the field 
B 1 ^ = t IJ KL e a: e b i n the action ([76]) - where e\ L is interpreted as the tetrad co-frame - produces 
the action of general relativity in four dimensions. In the physically relevant cases of the above 
examples one needs to deal with non-compact groups - G = SX(2,R) and G = SL(2, C) respectiv- 
ely. The non-compacteness of the gauge group leads to certain infrared divergencies of transition 
amplitudes (infinite volume factors). We avoid such complications at this stage and concentrate 
on the compact G case. The infinite volume divergences will be solved in the particular case of 
interest which is the spinfoam models for four dimensional Lorentzian gravity whose construction 
is reviewed in Section 

For the moment we assume Ai to be a compact and orientable. The partition function, Z, is 
formally given by 

Z = { V[B]V[uo] exp(i I (B A F{u))). (79) 
Formally integrating over the B field in (|79p we obtain 

Z = J V[u] 8 (F(u>)) . (80) 

The partition function Z corresponds to the 'volume' of the space of flat connections on A4. 

In order to give a meaning to the formal expressions above, we replace the d-dimensional 
manifold M. with an arbitrary cellular decomposition A. We also need the notion of the associated 
dual 2-complex of A denoted by A*. The dual 2-complex A* is a combinatorial object defined by 
a set of vertices v € A* (dual to d-cells in A) edges e £ A* (dual to (d— l)-cells in A) and faces 
/ £ A* (dual to (d— 2)-cells in A). In the case where A is a simplicial decomposition of M the 
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Figure 6: Left: a triangulation and its dual in two dimensions. Right: the dual two complex; faces 
(shaded polygone) are dual to O-simpleces in 2d. 




Figure 7: Left: a triangulation and its dual in three dimensions. Right: the dual two complex; 
faces (shaded wedge) are dual to 1-simpleces in 3d. 




Figure 8: Left: a triangulation and its dual in four dimensions. Right: the dual two complex; 
faces (shaded wedge) are dual to triangles in 4d. The shaded triangle dual to the shaded face is 
exhibited. 
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structure of both A and A* is illustrated in Figures El [3 and [8] in two, three, and four dimensions 
respectively. 

For simplicity we concentrate in the case when A is a triangulation. The field B is associated 
with Lie algebra elements Bf assigned to faces / E A*. We can think of it as the integral of the 
(d— 2)-form B on the (d— 2)-cell dual to the face / € A*, namely 



B f = J B. (81) 

(d-2)-ccll 

In other words Bf can be interpreted as the 'smearing' of the continuous (d— 2)-form B on the 
(d— 2)-cells in A. We use the one-to-one correspondence between faces / € A* and (d— 2)-cells in 
A to label the discretization of the B field B / . The connection u> is discretized by the assignment 
of group elements g e € G to edges e € A*. One can think of the group elements g e as the holonomy 
of ui along e G A*, namely 



<]<■ 



Pexp(- [ u), (82) 

J e 



where the symbol "P exp " denotes the path-order-exponential that reminds us of the relationship 
of the holonomy with the connection along the path e € A*. 

With all this the discretized version of the path integral (|79l) is 



Z(A) = f U dg e J] dB f e iB f u ' = f J[ dg e J[ S(g ei ■ ■ ■ g e J, (83) 

J eGA* /GA* J eGA* /GA* 

where Uf — g ei ■ ■ ■ g en denotes the holonomy around faces, and the second equation is the result of 
the B integration: it can be thus regarded as the analog of (|80]l 1^1 . The integration measure dBf 
is the standard Lebesgue measure while the integration in the group variables is done in terms of 
the invariant measure in G (which is the unique Haar measure when G is compact). For given 
h € G and test function F(g) the invariance property reads as follows 

J dgF(g) = J dgFig- 1 ) = J dgF(gh) = J dgF(hg) (84) 

The Peter- Weyl's theorem provides a useful formula or the Dirac delta distribution appearing in 
(l83|) . namely 

S(g) = J2 d p^[p(9% (85) 
p 

where p are irreducible unitary representations of G. From the previous expression one obtains 

Z ( A )= E / Il^e H d Pf Tr[p f (gl...g:)]. (86) 

C:{p}^{f} J eeA* / £A* 

Integration over the connection can be performed as follows. In a triangulation A, the edges 
e <G A* bound precisely d different faces; therefore, the g e 's in ([SB")) appear in d different traces. 
The relevant formula is 

P inv(Pir-- ,Pd) ■= I dg e Pl (g e )® p 2 {g e )®---® p d {g e ). (87) 



23 It is important to point out that the integration over the algebra valued B field does not exactly give the group 
delta function. For instance, in the simple case where G = SU(2) with B g su(2) integration leads to the SO (3) 
delta distribution (which only contains integer spin representations in the mode expansion H85II ). Generally, one 
ignores this fact and uses the G-delta distribution in the models found in the literature. 
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For compact G it is easy to prove using the invariance (and normalization) of the the integration 
measure that Pf nv = {Pf nv ) 2 i s the projector onto lnv[pi <g) p 2 <g> • • • ® Pd] • I n this way the 
spin foam amplitudes of 5*0(4) BF theory reduce to 

z BF (A)= yi n d «n ? ™^'--^ ( 8§ ) 

Cf.{f}^ P f /£A* eeA* 

In other words, the BF amplitude associated to a two complex A* is simply given by the sum over 
all possible assignments of irreducible representations of G to faces of the number obtained by the 
natural contraction of the network of projectors Pf nv according to the pattern provided defined by 
the two-complex A*. 

There is a nice graphical representation of the partition function of BF theory that will be 
very useful for some calculations. On the one hand, using this graphical notation one can easily 
prove the discretization independence of the BF amplitudes. On the other hand this graphical 
notation will simplify the presentation of the new spin foam models of quantum gravity that will 
be considered in the following sections. This useful notation was introduced by Oeckl 280, 285 
and used in [199] to give a general prove of the discretization independence of the BF partition 
function and the Turacv-Viro invariants for their definition on general cellular decompositions. 

Let us try to present this notation in more detail: The idea is to represent each representation 
matrix appearing in (|86|) by a line (called a wire) labeled by an irreducible representation, and 
integrations on the group by a box (called a cable). The traces in equation (|86j) imply that there 
is a wire, labelled by the representation pf, winding around each face / e A*. In addition, there 
is a cable (integration on the group) associated with each edge eeA*. As in ([951 . there is a 
projector P? nv is the projector into Inv[pi ® p 2 ® • • • <8> Pd] associated to each edge. This will be 
represented by a cable with d wires as shown in (|89|) . Such graphical representation allows for a 
simple diagrammatic expression of the BF quantum amplitudes. 



Pi 92 P3 Pd 



P inv(Pl:P2,P3, ■■■ ,Pd) = 




We can now express the product of projection operators defining the BF quantum amplitudes 
in equation (I95[) in an entirely graphical way. Let us show this in 2, 3 and 4 dimensions. The 
generalization to arbitrary dimension is obvious. 

According to Figure ([6]) in two dimensions edges e € A* are shared by two faces. This means 
that there are two representation matrices associated to the group element g e and hence two wires 
labelled by two in principle different pj. The endpoint of the open wires are connected to the 
neighbouring vertices and form close loops around each face (due to the trace in (I86[) V There is a 
cable with two wires on each e € A*. The BF amplitude is obtained by plugging together all the 
cables in the way dictated by the dual 2-complex A*. In two dimensions the result is 



z BF (A) = £ n d P 

Cf.{f}^p f /eA* 



where it is understood that the open ends of the wires are connected to neighbouring vertices 
according to the combinatorial structure provided by the cellular 2-complex A*. In other words, 
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the graph on the right hand side of the previous equation is just the basic building block of a closed 
cable- wire diagram involving the whole complex A. 

The construction is now obvious in three dimensions. The only difference is that now three 
faces share each edge and hence cables have three wires labelled by three different representations. 
The wires are connected according to the structure provided by the dual two complex A* (see now 
Figure [7]) so that a close wire-loop is obtained for each face / £ A* . The result is 



I I I 



z BF (A) = e n d p 

Cf.{f}^ Pf feA* 



pj 

V Pi 


\pi 


pK 


/ 3 3 



TXT 



At each vertex we have six wires corresponding to the six faces / £ A* sharing a vertex v £ A*. 
As in the 2d case, the cable-diagram diagram on the right is only a piece of a global diagram 
embracing the whole complex A. 

Finally, in four dimensions the construction follows the same lines. Now edges are shared by 
four faces; each cable has now four wires. The cable wire diagram giving the BF amplitude is 
disctated by the combinatorics of the dual two complex A* . From Figure [8] we get 




The 10 wires corresponding to the 10 faces / £ A* sharing a vertex v £ A* are connected to 
the neighbouring vertices through the 5 cables (representing the projectors in (|95p and Figure [Ml) 
associated to the 5 edges e £ A* sharing the vertex v £ A*. 

6.1 Special Cases 

We end this section with some simple examples of BF theory that will be useful in following 
applications. 

6.1.1 SU(2) BF theory in 2d: the simplest topological model 

The amplitude of BF theory in two dimensions can be entirely worked out due to the simplicity 
of the projectors Pf nv in this case. For concreteness here we take the special case G — SU(2) as 
this structure group will play an important case in the gravity context, while the techniques used 
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here will provide a nice warming-up exercise for what follows. Irreducible representations of SU(2) 
are labelled by spins (half- integers): from now on we set pf = jf £ N/2. In two dimensions, the 
relevant projectors are Pf nv : jf ® jf — > Inv[j/ <g> jf] which in turn are non trivial if and only if 
jf = j%. In our graphical notation the projector is 

j j' j 

±L a U 

KM') = Q =^ff T) n> 03) 

where the lines on the right hand side denote the identity representation matrices in the repre- 
sentation j. Replacing this in the expression of the BF amplitude given in equation ([9H|) one gets 




ZlUA) = £ I] + 1) II 79TTTT j ) > Q > ( 94 ) 

j f£A* eGA* 

where the sum over jf in (f9T)|) has collapsed to a single sum as according to Spin foam 

amplitudes vanish unless all the spins assigned to faces are equal, i.e. jf — j for all / G A*. In 
addition to a factor dj = (2j + 1) for each face / £ A* one gets a factor l/(2j + 1) for each edge 
e £ A* from the projectors fsee !93[) . Finally, for each vertex v £ A* one obtains a factor (2j + 1) 
coming from the remaining closed loop that is left over once (j§3")l is used in all the edges associated 
with v (see Finally, if we denote Nf, N e , and N v the number of faces, the number of 

edges, and the number of vertices in A* respectively; then it follows that the BF partition function 
becomes 

Z| d F (A) = £(2j + 1)"/-".+". = ^(2j + 1)*, (95) 

3 j 

where \ is th e Euler character of the two dimensional manifold M. . In this simple two dimensional 
example one discovers explicitly two important properties of the partition function (and quantum 
amplitudes) of BF theory: first, quantum amplitudes are independent of the regulating triangu- 
lation used in their definition, and second their discretization independent value depend only on 
global properties of the spacetime defining topological invariants of the manifold M.. One can 
genuinly write 

Zf F {M) = Zf F {A). (96) 



Extra remarks on two dimensional BF theory 

Two dimensional BF theory has been used as the basic theory in an attempt to define a manifold 
independent model of QFT in 245 . It is also related to gravity in two dimensions in two ways: 
on the one hand it is equivalent to the so-called Jackiw-Tcitclboim model [226, 345 , on the other 
hand it is related to usual 2d gravity via constraints in a way similar to the one exploited in four 
dimensions (see next section). The first relationship has been used in the canonical quantization 
of the Jackiw-Teitelboim model in [120] . The second relationship has been explored in [289J 



6.1.2 SU{2) BF theory and 3d Riemannian gravity 

In three dimensions edges e £ A* are shared by three faces and thus the relevant projector is 
represented by a cable with three wires. In the case G — SU(2) the singlet component Invjji ® 
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Therefore, we can write the projector as 
j m k 

(97) 

A\ 

where the three valent open graphs on the right hand side represent the unique normalized invariant 
vector in j ® m ® fc. Now we can play the same game as in the previous example. Replacing the 
expression of the projectors in equation (|91j) we get 



Z BF {A) = 




Hence, the amplitude can be written as a sum over spins jf associated to vertices of a product of 
face and vertex amplitudes, namely 



Z BF (A) = 




The vertex amplitude is a 6j-symbol and depends on six spins labelling the corresponding faces 
associated with a vertex. 

Extra remarks on three dimensional BF theory 

Three dimensional BF theory and the spin foam quantization presented above is intimately related 
to classical and quantum gravity in three dimensions (for a classic reference see |106j ). We will 
discuss this relationship in detail in Part IIVI The state sum as presented above matches the 
quantum amplitudes first proposed by Ponzano and Regge in the 1960s based on their discovery of 
the asymptotic expressions of the 6j symbols |303j and is often referred to as the Ponzano-Regge 
model. Divergences in the above formal expression require regularization. We will see in Part IIVI 
that natural regularizations are available and that the model is well defined [STJ 12751 1180] . For a 
detailed study of the divergence structure of the model see [94, 95, 96, 97]. The quantum deformed 
version of the above amplitudes lead to the so called Turaev-Viro model [358] which is expected to 
correspond to the quantization of three dimensional Riemannian gravity in the presence of a non 
vanishing positive cosmological constant. For the definition of observables in the latter context as 
well as in the analog four dimensional analog see [50] . 

The topological character of BF theory can be preserved by the coupling of the theory with 
topological defects that play the role of point particles. In the spin foam literature this has 



ji ® js] is one dimensional when it is non-trivial. 

j m k 
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been considered form the canonical perspective in [2731 1274] and from the covariant perspective 
extensively by Freidel and Louapre 182 . These theories have been shown by Freidel and Livine to 
be dual, in a suitable sense, to certain non-commutative helds theories in three dimensions [179, 



Concerning coupling BF theory with non topological matter see |166[ 1149] for the case of 
fermionic matter, and [341 for gauge helds. A more radical perspective for the definition of matter 
in 3d gravity is taken in |167j . For three dimensional supersymmetric BF theory models see 2 13 , 32 

Recursion relations for the 6j vertex amplitudes have been investigated in [S3] I151J . They 
provide a tool for studying dynamics in spin foams of 3d gravity and might be useful in higher 
dimensions [93] . 

6.1.3 SU{2) BF theory: The Oouguri model 

In four dimensions the structure is essentially the same. Now edges are shared by four repre- 
sentations. In the case of G = SU(2) the vector space Inv[ji (g) j 2 ® ja ® 3a] has a dimension 
that is generically greater than one. Consequently, if one wants to write the BF quantum parti- 
tion function as a sum over representations of products of faces and vertex amplitudes (the spin 
foam representation) one has to express the projectors in terms of a basis of invariant vectors in 
ji <8> 32 <£> 33 ® 3 a- Graphically, this can be done as follows 



where the four-valent diagrams on the right hand side represent elements of an orthonormal basis 
of invariant vectors in j\ <g> j<z ® j'3 ® ji and its dual, which are labelled by a half integer l. As in 
the previous lower dimensional cases we replace the previous equation in the expression (|92j) to 
obtain 



Eg. 



Ji 32 33 3a 



3i 32 33 3a 




(100) 




(101) 



which can be re-written in the spinfoam representation as 



39 



Zbf (a) = n ^ e n li 

Cf.{f) :!■ C e :{e}-H e 




J4 (102) 



The vertex amplitude in the previous expression is often called a 15j-symbol. 

6.1.4 SU(2) x SU(2) BF theory: a starting point for 4d Riemannian gravity 

For completeness we now present the BF quantum amplitudes in the case G — SU(2) x SU(2). 
This special case is of fundamental importance in the construction of the gravity models presented 
in the following sections. The construction mimics that of the previous section in all details. In fact 
the product form of the structure group implies the simple relationship Zbf(SU(2) x SU(2)) = 
Zbf(SU(2)) 2 . Nevertheless, it is important for us to present this example in explicit form as it will 
provide the graphical notation that is needed to introduce the gravity models in a simple manner. 
As in the previous cases the spin foam representation of the BF partition function follows from 
expressing the projectors in (|92[) in the an orthonormal basis of intertwiners. From the previous 
example and the product form of the structure group we have 



/jl P3 p4 



1 1 1 1 



■ — .— ■— .— .-j- ,_)_ 

Ji 32 33 3d Jj 32 3r 3a 



3l 32 33 3i Jl J2 33 Ji 

E " " < 103 ) 

"* /TK /TK 



where pj = jj ® j~f , and jf and 1^ are half integers labelling left and right representations of 
SU(2) that defined the irreducible unitary representations of G = SU{2) x SU(2). Accordingly, 
when replacing the previous expression in (|92p one gets 



Zbf(A) 



E 

:{/}- 



/-'/ 




(104) 



and equivalently 
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C/:{/}^P/ /6A* C e :{e}-> tc 




(105) 



from which we obtain finally the spin foam representation of the SU(2) x SU(2) partition function 
as a product of two SU(2) amplitudes, namely 

Cf-{f}^Pf /eA* C e :{e}-> te «GA* 




J 4 'I 




it (106) 



Extra remarks on four dimensional BF theory 



The state sum (|86j) is generically divergent (due to the gauge freedom analogous to (|78[)). A 
regularized version defined in terms of SU q (2) x SU q (2) was introduced by Crane and Yetter |124[ 
123 . As in three dimensions, if an appropriate regularization of bubble divergences is provided, 
(|86p is topologically invariant and the spin foam path integral is discretization independent. 

As in the three dimensional case BF theory in any dimensions can be coupled to topological 
defects [37]. In the four dimensional case defects are string-like |169j and can carry extra degrees 
of freedom such as topological Yang-Mills fields [269 . The possibility that quantum gravity could 
be defined directly form these simple kind of topological theories has also been considered outside 
spin foams |344j (for which the UV problem described in the introduction is absent) is attractive 
and should, in my view, be considered further. 

It is also possible to introduce one dimensional particles in four dimensional BF theory and 
gravity as shown in |174j . 



6.2 The coherent states representation 

In this section we introduce the coherent state representation of the SU(2) and Spin(4) path 
integral of BF theory. This will be particularly important for the definition of the models defined 
by Freidel and Krasnov in [176] that we will address in Section [9] as well as in the semiclassical 
analysis of the new models reported in Section 1131 The relevance of such representation for spin 
foams was first emphasized by Livine and Speziale in [246J . 
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6.2.1 Coherent states 



Coherent states associated to the representation theory of a compact group have been studied by 
Thiemann and collaborators [351] [M9] [3351 EH HI] 1313 [356] |3521 [46] ST] [171] see also [70]. Their 
importance for the new spin foam models was put forward by Livine and Speziale in [246] where 
the emphasis is put on coherent states of intertwiners or the so-called quantum tetrahedron (see 
also |117j ). Here we follow the presentation of |176j . 

In order to built coherent states for Spin(4) we start by introducing them in the case of SU(2). 
Starting from the representation space J#j of dimension dj = 2j + 1 one can write the resolution 
of the identity in terns of the canonical orthonormal basis \j, m) as 

lj =$>\m)0\m|, (107) 

m 

where —j<m< j. There exists an over complete basis \j,g) £ J#j labelled by g £ SU{2) such 
that 

li = d j J dg\j,g){j,g\, (108) 

SU(2) 

The states \j, g) £ are SU (2) coherent states defined by the action of the group on maximum 
weight states \j,j) (themselves coherent), namely 

\j,g) = g\j,j) = ]T \j,m)Di nj (g), (109) 

m 

where ^ (g) are the matrix elements of the unitary representations in the \j, m) (Wigner matri- 
ces). Equation (|108l) follows from the orthonormality of unitary representation matrix elements, 
namely 



dj / dg \j,g)(j,g\,= dj £ \j,m)(j,m'\ / dgD^Di (g) = ]T \j,m){j,m\, (110) 

where in the last equality we have used the orthonormality of the matrix elements graphically 
represented in (1§3"]) . The decomposition of the identity (|108p can be expressed as en integral on 
the two-sphere of directions S 2 — SU(2)/U(1) by noticing that D 3 m -{g) and D 3 m -{gh) differ only 
by a phase for any group element h from a suitable U(l) C SU(2). Thus one has 



lj = dj / dn\j,n){j,n\, (111) 
Js 2 

where n £ S 2 is integrated with the invariant measure of the sphere. The states \j, n) form (an 
over-complete) basis in Jffj. SU(2) coherent states have the usual semiclassical properties. Indeed 
if one considers the operators J % generators of su{2) one has 

(j,n\J%n)=jn i , (112) 

where n l is the corresponding three dimensional unit vector for n £ S 2 . The fluctuations of J 2 are 
also minimal with A J 2 = h 2 j, where we have restored h for clarity. The fluctuations go to zero 
in the limit h — > and j — > oo while hj is kept constant. This kind of limit will be used often as 
a notion of semiclassical limit in spin foams. The state \j, n) is a semiclassical state describing a 
vector in M 3 of length j and of direction n. It will convenient to introduce the following graphical 
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notation for equation (| 111 |) 



dn 



si 



(113) 



Finally, an important property of SU (2) coherent states stemming from the fact that 



1 3l j) I 2 ? 2 ) I 2 5 2 ) ' ' ' I 2 ? 2 ) — I 2 1 2 } 



is that 

|j,n) = |i,< 2j - (114) 

The above property will be of key importance in constructing effective discrete actions for spin 
foam models. In particular, it will play a central role in the study of the semiclassical limit of the 
EPRLand FK modesl studied in Sections [JJ El and O In the following subsection we provide an 
example for Spin(A) BF theory 



6.2.2 Spin(A) BF theory: amplitudes in the coherent state basis 

Here we study the coherent states representation of the path integral for Spin(4) BF theory. The 
construction presented here can be extended to more general cases. The present case is however of 
particular importance for the study of gravity models presented in Sections [JJ [SJ and With the 
introductions of coherent states one achieved the most difficult part of the work. In order to express 
the Spin{4) BF amplitude in the coherent state representation one simply inserts a resolution of 
the identity in the form (jllll) on each and every wire connecting neighbouring vertices in the 
expression (| 1041) for the BF amplitudes. The result is 




(115) 



where we have explicitly written the n± £ S 2 integration variables only on a single cable. One 
observes that there is one n± € S* 2 per each wire coming out at an edge e € A*; as wires are in 
one-to-one correspondence with faces / e A* the integration variables n^y € S 2 are labelled by 
an edge and face subindex. In order to get an expression of the BF path integral in terms of an 
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affective action we restore at this stage the explicit group integrations represented by the boxes in 
the previous equation. One gets, 

Zbf(A) = £ J] d.-d j; / n d,- d jtf dn; f dni f 

C/:{/}^P/ /£A* J eG A* 

II II ^^(K,!^ (116) 

where we have used the coherent states property (| 1 14[) . and |n ± ) is a simplified notation for ||, n^). 
Also we denote by gf v the group elements associated to the five boxes for each given vertex (see 
eq. (|115[) ). The previous equation can be finally written as 

Z BF (A)= £ II d iJ d ^ / II d j; f d jt f dn ef d <f d 9evd9+ exp(^ ±)n± [ 9 ± ]), (117) 

Cf-{f}^Pf /eA* J eG A* 

where the discrete action 

^±n ± [5 ± ]= E ^.nj^] (H8) 

with 

5 

-S'J.nls] = E ^bln^abl^ 1 ^! n 6a ), (119) 

a<b=l 

and the indices a, 6 label the five edges of a given vertex. The previous expression is exactly equal 
to the form (I55|) of the BF amplitude. In the case of the gravity models studied in what follows, the 
coherent state path integral representation will be the basic tool for the study of the semiclassical 
limit of the models and the relationship with Regge discrete formulation of general relativity. 

7 The Riemannian EPRL model 

In this section we introduce the Riemannian version of the Engle-Pereira-Rovelli-Livine (EPRL) 
|162| 1158] , The section is organized as follows. The relevant Spin(A) representation theory is 
introduced in Section 17.11 In Section 17.21 we present and discuss the linear simplicity constraints 
- classically equivalent to the Plcbanski constraints - and discuss their implementation in the 
quantum theory. In Section [7T3l we introduce the EPRL model of Riemannian gravity. In Section [7T4l 
we prove the validity of the quadratic Plebanski constraints - reducing BF theory to general 
relativity - directly in the spin foam representation. In Section T7.5l we discuss a certain modification 
of the EPRL model. In Section 171)1 we present the coherent state representation of the EPRL model 
which is essential for the semiclassical analysis of Section [T3] The material of this section will also 
allow us to describe the construction of the closely related (although derived from a different 
logic) Riemannian FK constructed in |176j . The idea that linear simplicity constraints are more 
convenient for dealing with the constraints that reduce BF theory to gravity was pointed out by 
Freidel and Krasnov in this last reference. We have shown in Part Ql] that they arise directly from 
the canonical analysis of the Plebanski formulation of gravity. 

7.1 Representation theory of Spin(A) and the canonical basis 

In order to introduce the definition of the EPRL model we need to briefly review the representation 
theory of Spin(A). The group Spin(4) = SU(2) x SU(2) with Lie algebra spin(A) = su(2) © su(2) 
can explicitly be described in terms of generators Jj_ such that 

[4,4]=eV±- (120) 
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Unitary irreducible representations ^+ j- of Spin(A) are given by the product of unitary irre- 
ducible represetantions of SU(2), i.e., they are labelled by two half-integers j . A standard basis 
in J^j+ j- is given by vectors \j + , j~ , m + , m~) of eigenstates of the Casimirs C\ = J? + J 2 and 
C2 = J\ — </ 2 and the components Jl and , explicitly 

Ci\j + ,j~ ~,m + ,m~) = + 1) + 1)) \j + ,™ + 

C2\j + ,j~,m + ,m-) = + 1) + 1)) \j + , j~,m + ,m~) 

J±\j + ,j~, rn + , m~) = rtv^ m + , mT). (121) 

The definition of the EPRL model requires the introduction of an (arbitrary) subgroup SU (2) C 
Spin(A). This subgroup will be shown, a posteriori , to be the internal gauge group of the gravita- 
tional phase space in connection variables (appearing in the classical canonical study of Section [321 
and leading to the quantization of Section |4]). In the quantum theory, the representation theory of 
this SU(2) subgroup will be hence important. This importance will soon emerge as apparent from 
the imposition of the constraints that define the EPRL as well as the FK model of Section [9] 

The link between the unitary representations of spin(A) and those of su(2) comes from the fact 
that the former can be expressed as a direct sum of the latter according to 

j + +r 

j=\j + -j~\ 

As the unitary irreducible representations of the subgroup SU(2) € Spin(A) are essential in un- 
derstanding the link of the EPRL model and the operator canonical formulation of LQG it will be 
convenient to express the action of the generators of the Lie algebra spin (4) in a basis adapted to 
the above equation. In order to do this we first notice that the Lie algebra (1120)) can be equiva- 
lcntly characterized in terms of the generators of a rotation subgroup L l and the remaining "6oos£" 
generators K % as follows 

[L 3 ,L ± ]=±L ± [I + ,I-]=2I 3 
[L + ,K + ] = [L_,K_} = [L 3 ,K 3 ] = 

[K 3 ,L ± ] =±K± [L±,K T ] = ±2 K 3 [L 3 ,K ± ] = ± K ± 

[K 3 ,K±] =±L ± [K+,KJ\ = 2 L 3 , (123) 

where K± = K 1 ± iK 2 and L± = L Y ±iL 2 respectively. The action of the previous generators in 
the basis \j + ,j~,j, to) can be shown to be 

L 3 \j + ,j~,j,m) = m\j + ,j-,j,m), 

L+ \j + ,j~,j,m) = V U + m + 1)0' - m ) \i + ,3~, m + 1), 
L~\j + ,j~,j,m) = V (j + m ) U - m + 1 ) \j + 1 j ~ i J, m - !)) 
K 3 \j + ,r,j,m) = 

a jVf - m2 \j + >j~>j - l- m ) + 7j"i|j + . j~,3,™>) - a j+1 y/(j + l) 2 - TO 2 |j+, j~, j + l,m), 

K + \j + ,j~,j,m) = a ^(j - m)(j - to - - l,m+ 1) 

+7j V U -m)(j + m + 1) | j + , j~ , j, to + 1) 

+a j+ i v /(.?+m + l)(j + TO + 2)|j+ , j- , j + 1, m + 1) , 

K~\j + ,j~,j,m) = -a, \/ {j + m)(j+m- 1) | j+, - l,m - 1) 

+lj V(j + m )(j - 171 + l )\j + ,3~, 3, m-1) 

-a j+ i y/(J - TO + 1) (j - TO + 2) | j+, j- , j + 1, TO - 1) , (124) 
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where 



j + (j + + i)-rcr + i) l (P - (j + + j- + i) 2 )(.f - (j + -rn mn 

7j = TcJTij a ^V ' (125) 

The previous equations will be important in what follows: they will allow for the characterisation 
of the solutions of the quantum simplicity constraints in a direct manner. This concludes the 
review of the representation theory that is necessary for the definition of the EPRL model. 



7.2 The linear simplicity constraints 

As first shown in |176] , the quadratic simplicity constraints (|25p - and more precisely in their dual 
version presented below (|154l) - are equivalent in the discrete setting to the linear constraint on 
each face of a given tetrahedron 

D) = L) - l -K) « 0, (126) 

where the label / makes reference to a face / € A*, and where (very importantly) the subgroup 
SU(2) C Spin(4) that is necessary for the definition of the above constraints is chosen arbitrarily 
at each tetrahedron, equivalent on each edge e £ A*. Such choice of the rotation subgroup is 
the precise analog of the time gauge in the canonical analysis of Section [31 Moreover, the linear 
simplicity constraint above is the simplicial version of the canonical primary constraint (|3T[) . The 
EPRL model is defined by imposing the previous constraints as operator equations on the Hilbert 
spaces defined by the unitary irreducible representations of Spin(A) that take part in the state-sum 
of BF theory. We will show in Section 17.41 that the models constructed on the requirement of a 
suitable imposition of the linear constraints (|126l) satisfy the usual quadratic Plebanski constraints 
in the path integral formulation (up to quantum corrections which are neglected in the usual 
semiclassical limit). 

It will be useful in what follows to write the above constraint in its equivalent form 

D) = (1 - 7 ) Jp - (1 + l)Jj l « 0, (127) 

From the commutation relations p20[) of previous section we can easily compute the commu- 
tator of the previous tetrahedron constraints and conclude that in fact it does not close, namely 



The previous commutation relations imply that the constraint algebra is not closed and cannot 
therefore be imposed as operator equations on the states summed over in the BF partition function 
in general. There are two interesting exceptions to the previous statement: 

1. The first one is to take 7 = ±1 in which case the constraint (|187|) can be imposed strongly. 
More precisely, for the plus sign it reduces to the condition J| + = (or equivalently for the 
other sign Jj_ = 0) which amounts to setting all the left quantum numbers in the BF state- 
sum for the group Spin(4) to zero. It is immediate to see that if we impose such restriction 
on the Spin(A) BF amplitude (1106)) it simply reduces to the amplitude (|102p defining SU(2) 
BF theory. Thus, this possibility does not lead to an acceptable model for quantum gravity 
in 4d. 
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2. The second possibility is to work in the sector where — 0. We will show later that this 
corresponds to the famous Barret-Crane model [53] with the various limitations discussed in 
Section [TDJ 

The EPRL model is obtained by restricting the representations appearing in the expression of 
the BF partition function so that at each tetrahedron the linear constraints (|18T[) the strongest 
possible way that is compatible with the uncertainties relations stemming from (|188[) . In addition 
one would add the requirement that the state-space of tetrahedra is compatible with the state- 
space of the analogous excitation in the canonical context of LQG so that arbitrary states in the 
kinematical state of LQG have non trivial amplitudes in the model. 



On the weak imposition of the linear simplicity constraints 

We now discuss the weak imposition of the linear simplicity constraints in the quantum framework. 
There are essentially three views in the literature: two of them, discussed below, concerns directly 
the way the EPRL model has been usually introduced. The third possibility is the semiclassical 
view based on the coherent state representation leading to the FK model (see Section [§]) . 

7.2.1 Riemannian model: The Gupta Bleuler criterion 

Due to the fact that the constraints do not form a closed (first class) algebra in the generic case 
one needs to devise a weaker sense in which they are to be imposed. One possibility is to consider 
the Gupta-Bleuler criterion consisting in selecting a suitable class of states for which the matrix 
elements on vanish. In this respect one notices that if we chose the subspace Jtfj C J#j+ j- one 
has 

(j + ,r,j,q\D 3 f \j + ,j-,j,m) =S 9 , m m(l - 



(j + ,r,J,q\Df\j + ,r,j,m) = S q±hmy /(J ±m + T m)(l - % 

One immediately observes that matrix elements of the linear constraints vanish in this subclass if 
one can chose 

There are two cases: 

1. Case 7 < 1: Following |136j . in this case one restricts the Spin(4) representations to 

^ = ^±7)1, (130) 

which amounts to choosing the maximum weight component j = y + j~ in the expansion 
(|183p . Simple algebra shows that condition (|189[) is met. There are indeed other solutions 
[134] of the Gupta-Bleuler criterion in this case. 

2. Case 7 > 1: In this case [13] one restricts the Spin(4) representations to 

] ± = { 1 ±l)L + 1 —±, (131) 



which amounts to choosing the minimum weight component j = j + — j in the expansion 
(|183|) . This choice of j is the solution to condition (|189l) . 
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7.2.2 Riemannian model: The Master constraint criterion 



Another criterion for weak imposition can be developed by studying the spectrum of the Master 
constraint Mf = DfDf. Strong imposition of the constraints would amount to looking for the 
kernel of the master constraint Mf. However, generically the positive operator associated with the 
master constraint does not contain the zero eigenvalue in the spectrum due to the open nature of 
the constraint algebra (|188|) . The proposal of [160] . that we follow here, is to look for the minimum 
eigenvalue among spaces € Explicitly 

Mf\i/j>=m j ±Jip> ) (132) 

where 

m J±J - (1 - lfj+{3+ + !) + (! + 7 2 )j-(i- + 1) - (1 - l 2 )\j(j + 1) - i+(i+ + 1) - + 1)] 

There are two cases: 

1. Case 7 < 1: In this case the minimum eigenvalue is obtained for 

3 ± = (1 ± 7 )i/2, (133) 

i.e., one has 

m j ± ,j > m j 

where rrij = h 2 (l — j 2 )j where we have restored the explicit dependence on H so that it 
is apparent that the selected eigenvalue vanishes in the the semiclassical limit h — > and 
k — > oo with hk =constant. 

2. Case 7 > 1: In this case the minimum eigenvalue is obtained for 

, 7 ± = ( 7 ±1)j72, (134) 

for which rrij = h 2 ^ 2 — l)j and and thus vanishes in the the semiclassical limit h — > and 
k — > 00 with hk =constant. 

The master constraint criterion will be used here as the basic justification for the definition of 
the EPRL model. Notice however, that the master constraint is only invariant under the subgroup 
SU(2) C Spin(A) chosen at each tetrahedron. Its lack of Lorentz invariance is clear from the 
fact that its validity is justified from the canonical analysis in the time gauge. Nevertheless, as 
shown in [325], one can also write down a criterion (equivalent to the master constraint one) in 
a manifestly Spin{A) invariant way. This reflects the arbitrariness in the choice of the internal 
subgroup SU(2) C Spin{4) and implies the consistency of the time gauge. In order to get a gauge 
invariant constraint one starts from the master constraint and uses the D l j = classically to rewrite 
it entirely in terms of the Spin(4) Casimirs, namely 

M f = (l- 7 ) 2 4 + (l+j) 2 J* -2(1-7 2 )J+- J_ 

* (l+7 2 )C 2 -2C l7 , (135) 

where C\ and C2 are the two spin{A) Casimirs given in equation (I121[) . The Lorentz invariant 
version of the master constraint is hence 

My = (1 + j 2 )C 2 - 2Ci7 = (136) 

One can verify for the cases 7 < 1 and 7 > 1 the Lorentz invariant criterion is just equivalent to 
the master constraint analysis above. 
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7.2.3 Riemannian model: Restrictions on the Immirzi parameter 

We have just shown how the states in J^+ j- can be restricted in order to satisfy the linear 
simplicity constraints according to different criteria. From all these the EPRL model is defined 
most simply using the Master constraint criterion. In all cases the restriction takes the form of the 
selection of a subspace M'j C given by unitary irreducible representations of a subgroup 

SU(2) £ Spin(A). As we will show in what follows this implies that the boundary states of the 
spin foam models defined by the EPRL construction are in correspondence to SU(2) spin network 
states: the states of the kinematical Hilbert space of LQG reviewed in Part|TIJ This is very nice as 
it suggests that the spin foam amplitudes can be interpreted as dynamical transition amplitudes 
of the canonical quantum theory. Now in order for this to be so, one would like arbitrary spins 
j £ N/2 to be admissible in the state sum. We shall see now that this imposes restrictions on the 
possible values of the Immirzi parameter. However, these restrictions seem so far of little physical 
importance as they are a special feature of the Riemannian model that is absent in its Lorentzian 
relative. As above there are two distinct cases: 

1. Case 7 < 1: If one wants the spin j £ N/2 to be arbitrary then the only possibilities are 
7 = or 7 = 1. Rational 7 values would restrict the spins j to a subclass in N/2. This 
restriction is not natural from the viewpoint of LQG. ITs relevance if any remains misterious 
at this stage. 

2. Case 7 > 1: In this case 7 £ N would allow any j £ N/2. Rational choices of 7 have the 
same effect as in the previous case. 

7.2.4 Riemannian model: Overview of the solutions of the linear simplicity con- 
straints 

Wc have reviewed different criteria given in the literature for the definition of the subset of quantum 
states that are selected in the construction of the EPRL model. From all the criteria the Master 
constraint one is the most direct and clear-cut. For the Riemannian model the result is: 

1. Case 7 < 1: 

j ± = (l± 1 )j/2, (137) 

2. Case 7 > 1: 

i ± = ( 7 ±l)j/2, (138) 

where j £ N/2 is arbitrary - except from the requirement that j^ 1 £ N/2 which depends of the 
value of the Immirzi parameter according to the previous section. From now on we denote the 
subset of admissible representation of Spin(A) as 

c Irrep(Spm(4)). (139) 

It follows from the discussion in the previous sections that for states ^> £ Jfcy satisfy the constraints 
(|187|) in the following semiclassical sense: 

(K) - 1 L))^> = e sci (140) 

where the symbol G sc (order semiclassical) denotes a quantity that vanishes in limit h — > 0, j — > 00 
with hj =constant. Equivalently one has 

[(i- 7 )j;-(i + 7 )jl]* = ^ c . (i4i) 
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This last equation can be written graphically as: 




The previous equation will be of great importance in the graphical calculus that will alow us to 
show that the linear constraint imposed here at the level of states imply the vanishing of the 
quadratic Plebanski constraints (|25[) and their fluctuations, computed in the path integral sense, 
in the appropriate large spin semiclassical limit. 



7.3 Presentation of the EPRL amplitude 

Here we complete the definition of the EPRL models by imposing the linear constraints on the BF 
amplitudes constructed in Section [6j We will also show that the path-integral expectation value 
of the Plebanski constraints (l25j) . as well as their fluctuations, vanish in a suitable semiclassical 
sense. This shows that the EPRL model can be considered as a lattice definition of the a quantum 
gravity theory. 

We start with the Riemannian model for which a straightforward graphical notation is available. 
The first step is the translation of Equation (|122p in terms of the graphical notation introduced in 
Section[6] Concretely, for 7 < 1 one has j ± = (1 ± j)j/2 <G J£y — as defined in (|139|) — becomes 



(1-7)! 



(1 + 7)1 



e 



(1-7); 



(1+7)! 



a 



(143) 



(l-7f! 



V (l+7)| 



For 7 > 1 we have 



(7- 



J(l + 7)| 



(7- 



(1+7): 



e 



(7- 



T) 1 



(144) 



U + 7)5 



The implementation of the linear constraints of Section (|7.2j) consist in restricting the represen- 
tations pf of Spin(A) appearing in the state sum amplitudes of BF theory as written in Equation 
(|104[) to the subclass pf € J^y C Irrep(»!spm(4)), defined above, while projecting to the highest 
weight term in (| 143[) for 7 < 1 . For 7 > 1 one must take the minimum weight term in (|144p . The 
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action of this projection will be denoted Wj : ^(1+7)5/2, |(i-7)h72 — > graphically 




(145) 



Explicitly, one takes the expression of the BF partition function (|95[) and modifies it by replacing 
the projector Pf nv (pi, • • • , p±) with p\, ■ ■ ■ p^ € by a new object 

*%r! (71, • • • . k) = Pinv (Pi ■ ■ ■ Pi) ® • • • ® »$4 (PI — P4) (146) 

with jx, • • • J4 € N/2 implementing the linear constraints described in the previous section. Graphi- 
cally the modification of BF theory that produces the EPRL model corresponds to the replacement 




(147) 



on the expression (|104p , where we have dropped the representation labels from the figure for 
simplicity. We have done the operation (|145|) on each an every of the four pairs of representations. 
The Spin(4) integrations represented by the two boxes at the top and bottom of the previous 
graphical expression restore the full Spin(A) invariance as the projection (|145p breaks this latter 
symmetry for being based on the selection of a special subgroup SU(2) C Spin(4) in its definition 
(see section [7j5] for an important implication). One should simply keep in mind that green wires in 
the previous two equations and in what follows are labeled by arbitrary spins j (which are being 
summed over in the expression of the amplitude ([148[) ), while red and blue wires are labelled by 
j + = (1 + r f)j/2 and j~ = |1 — r y\j/2 respectively. With this (I104|) is modified to 

zg rl (A) = J2 II d n-7ii d (i+7)iII p ^0-i.--- ,u) = 

pfEJf /6A' e 



- e rid 

Pf ejf /eA* 




(148) 



The previous expression is defines the EPRL model amplitude. 



7.3.1 The spin foam representation of the EPRL amplitude 

Now we will work out the spin foam representation of the EPRL amplitude which at this stage 
will take no much more effort than the derivation of the spin foam representation for Spin(A) BF 
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theory as we went from equation (|104|) to (|106|) in Section |6l The first step is given in the following 
equation 




which follows basically from the invariance of the Haar measure (I84[) (in the last line we have used 
(jlOOp ). More presicely, the integration of the subgroup SU{2) G Spin(4), represented by the green 
box on the right, can be absorbed by suitable redefinition of the integration on the right and left 
copies of SU(2), represented by the red and blue boxes respectively. With this we can already 
write the spin foam representation of the EPRL model, namely 




(150) 



where the vertex amplitude (graphically represented) depends on the 10 spins j associated to the 
face- wires and the 5 intertwiners associated to the five edges (tetrahedra) . As in previous equations 
we have left the spin labels of wires implicit for notational simplicity. We can write the previous 
spin foam amplitude in another form by integrating out all the projectors (boxes) explicitly. Using, 
(fM)]) we get 
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thus replacing this in (j!48[) we get 



^h(a)=e n vniWE n ( i52 ) 

i/ /gA* t e veA* 




form already in (|15ip . more explicitly 




(153) 



The previous is the form of the EPRL model as derived in [158] . 
7.4 Proof of validity of the Plebanski constraints 

In this section we prove that the quadratic constraints are satisfied in the sense that their path 
integral expectation value and fluctuation vanish in the appropriate semiclassical limit. 
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7.4.1 The dual version of the constraints 

In this section we rewrite the Plebanski constraints (|25p in an equivalent version where spacetimc 
indices are traded with internal Lorentz indices, namely 

t lJKL B\iB% L - e e^ pa « 0. (154) 

The constraints in this form are more suitable for the translation into the discrete formulation. 
More precisely, according to (f8Tj) . the smooth fields is now associated with the discrete quan- 
tities -B t / r f angl( , s , or equivalently Bj J as, we recall, faces / € A* are in one-to-one correspondence 
to triangles in four dimensions. The constraints (|154[) are local constraints valid at every space- 
time point. In the discrete setting, spacetime points are represented by four-simplexes or (more 
addapted to our discussion) vertices v £ A*. With all this the constraints (|154j) are discretized as 
follows: 

Triangle (or diagonal) constraints: £ijklBI j Bf L = 0, (155) 

for all / € v, i.e., for each and every face of the 10 possible faces touching the vertex v. 

Tetrahedron constraints: cijklB 1 / Bf, L = 0, (156) 

for all /, /' € v such that they are dual to triangls sharing a one-simplex, i.e., belonging to the 
same tetrahedron out of the five possible ones. 

4-simplex constraints: ^uklB 1 ^ B¥ L = e v , (157) 

for any pair of faces /, / G v that are dual to triangles sharing a single point. The last constraint 
will require a more detailed discussion. At this point let us point out that the constraint (|15T|) is 
interpreted as a definition of the four volume e v of the four-simplex. The constraint requires that 
such definition be consistent, i.e., the true condition is 

cuklBj ' Bf L — eijKhB 1 ^ Bp L = euKLBj'/, Bp^ = ••• = &„ (158) 

for all five different possible pairs of / and / in a four simplex, and where we assume the pairs /-/ 
are ordered in agreement with the orientation of the complex A*. 

7.4.2 The path integral expectation value of the Plebanski constraints 

Here we prove that the Plebanski constraint are satisfied by the EPRL amplitudes in the path 
integral expectation value sense. 

The triangle constraints: 

We start from the simplest case: the triangle (or diagonal) constraints (|155|) . We choose a face 
/ € v (dual to a triangle) in the cable- wire-diagram of Equation (|148|) . This amounts to choosing a 
pair of wires (right and left representations) connecting two nodes in the vertex cable wire diagram. 
The two nodes are dual to the two tetrahedra - in the four simplex dual to the vertex - sharing 
the chosen triangle. From (|29[) one can show that 

e IJKL ByBf L oc (1 + l) 2 JJ ■ Jj - (1 - 7) 2 4 ' Jf, (159) 
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where J* denotes the self-dual and anti-self-dual parts of II £ . The path integral expectation value 
of the triangle constraint is then 



((1 + 7 ) 2 JJ ■ JJ - (1 - ifJt ■ Jj) cx (160) 




where the double graspings on the anti-self-dual (blue) wire and the self-dual (red) wire represent 
the action of the Casimirs JJ ■ JJ and J^ ■ J^ on the cable-wire diagram of the corresponding 
vertex. Direct evaluation shows that the previous diagram is proportional to h jf which vanishes in 
the semiclassical limit h — > 0, j — > oo with hj ^constant. We use the notation already adopted in 
(|140p and call such quantity ff sc . This concludes the proof that the triangle Plebanski constraints 
are satisfied in the semiclassical sense. 

The tetrahedra constraints: 

The proof of the validity of the tetrahedra constraints (|156p . In this case we also have 




= ff sc . (161) 



where we have chosen an arbitrary pair of faces. In order to prove this let us develop the term on 
the right. The result follows from 
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(162) 



where in the first line we have used the fact that the double grasping can be shifted through the 
group integration (due to gauge invariance (|84|l ). and in the first and second terms on the second 
line we have used Equation (I142[) to move the graspings on self-dual wires to the corresponding 
anti-self-dual wires. Equation (|161|) follows immediately from the pervious one; the argument 
works in the same way for any other pair of faces. Notice that the first equality in Equation 
(|162| implies that we can view the Plebanski constraint as applied in the frame of the tetrahedron 
as well as in a Lorentz invariant framework (the double grasping defines an intertwiner operator 
commuting with the projection Pf nv represented by the box). An analogous statement also holds 
for the triangle constraints f|14.1[) . 

The 4-simplex constraints 

Now we show the validity of the four simplex constraints in their form (|158| . As we show below, 
this last set of constraints follow from the Spin(A) gauge invariance of the EPRL node (i.e., the 
validity of the Gauss law) plus the validity of the tetrahedra constraints (|156|) . Gauge invariance 
of the node takes the following form in graphical notation 




(163) 



where the above equation represents the gauge invariance under infinitesimal left SU (2) rotations. 
An analogous equation with insertions on the right is also valid. The validity of the previous 
equation can again be related to the invariance of the Haar measure used in the integration on the 
gauge group that defines the boxes (|85|). 

Now we chose an arbitrary pair / and / (where, recall, / is one of the three possible faces whose 
dual triangle only shares a point with the corresponding to /) and will show how the four volumen 
e v defined by it equals the one defined by any other admissible pair. The first step is to show that 
we get the same result using the pair /-/ and /-/, where / is another of the three admissible faces 
opposite to /. The full result follows from applying the same procedure iteratively to reach any 
admissible pair. It will be obvious from the treatment given below that this is possible. Thus, for 
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a given pair of admissible faces we have 




where going from the first line to the second and third lines we have simply used (|163|) on the 
bottom graspings on the right and left wires. The last line results from the validity of (|156p : notice 
that the second terms in the second and third lines add up to SC as well as the third terms in the 
second and third line. There is an overall minus sign which amounts for an orientation factor. It 
should be clear that we can apply the same procedure to arrive at any admissible pair. 
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7.5 Modifications of the EPRL model 



7.5.1 P ep ri is not a projector 

Let us study in a bit more detail the object P^ prl ( ji , • • • ,ji)- We see that it is made of two 
ingredients. The first one is the projection to the maximum weight subspace Jifj for 7 > 1 in the 
decomposition of j- for = (1 ± j)j/2 (j^ — (7 ± l)j/2 for 7 > 1) in terms of irreducible 
representations of an arbitrarily chosen SU(2) subgroup of Spin{4). The second ingredient is to 
eliminate the dependence on the choice of subgroup by group averaging with respect to the full 
gauge group Spin(A). This is diagramaticaly represented in ()147|) . However PjL r ;(ji, ■ ■ • , ]a) is not 
a projector, namely 

K v h (ii> • • • , kf + P^ri (h, ■ • • j k)- (165) 

Technically this follows from (|146l) and the fact that 

&(Pi---P4),(^ 1 ®---®^- 4 )]^0 (166) 

i.e., the projection imposing the linear constraints (defined on the frame of a tetrahedrom or edge) 
and the Spin(A) (or Lorentz) group averaging - rendering the result gauge invariant - do not 
commute. The fact the P^ pr i(ji, • • • , J4) is not a projection operator has important consequences 
in the mathematical structure of the model: 

1. From (|148[) one can immediately obtain the following expression for the EPRL amplitude 

Z eprl (A)= J2 II d |i-7li d (i+7)in P 4ri(iir--.i4). (167) 

Ps e.x feA* ' e 

This expression has the formal structure of expression (|95|) for BF theory. The formal simi- 
larity however is broken by the fact that -P e e pr ;(ji, ■ • • , ji) is not a projection operator. From 
the formal perspective is the possibility that the amplitudes be defined in term of a network 
of projectors (as in BF theory) might provide an interesting structure that might be of rele- 
vance in the definition of a discretization independent model (see discussion in Part [Vj) . On 
the contrary, the failure of -P e e pr ;(ji, • • • , ji) to be a projector may lead, in my opinion, to 
difficulties in the limit where the complex A is refined: the increasing of the number of edges 
might produce either trivial or divergent amplitudes 0. 

2. Another difficulty associated with Pe prl (ji, ■ • • , J4) 2 ^ P^prliki ' ' ' >k) is the failure of the 
amplitudes of the EPRL model, as defined here, to be consistent with the abstract notion of 
spin foams as defined in [34]. This is a point of crucial importance under current discussion 
in the community. The point is that the cellular decomposition A has no physical meaning 
and is to be interpreted as a subsidiary regulating structure to be removed when computing 
physical quantities. Spin foams (as defined in Section [S]) can fit in different ways on a given 
A, yet any of these different embeddings represent the same physical process (like the same 
gravitational field in different coordinates). Consistency requires the spin foam amplitudes to 
be independent of the embedding, i.e., well defined on the equivalence classes of spin foams as 
defined by Baez in |34j (the importance of these consistency requirements was emphasized in 
[54]). The amplitude (|167p fails this requirement due to -P e e pr ;(ji, ■ • • , ji) 2 ^ Pt w i{j\i ••' , J4 ) ■ 

24 This is obviously not clear from the form of II167I I. We are extrapolating the properties of (Pf rl ) N for large N 
to those of the amplitude (1 1 6 7ft in the large number of edges limit implied by the continuum limit. 
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7.5.2 The Warsaw proposal 



The above difficulties have a simple solution in the Riemannian case. As proposed in 45, 229 one 
can obtain a consistent modification of the EPRL model by replacing P^ pr i in f|16T[) by a genuine 
projector P^, graphically 



Kpriih ■ ■ ■ k) 




a/3 





(168) 




where 

9l = = Averse 

It is easy to check that by construction 

The variant of the EPRL model proposed in 1229] takes then the form 

z eprl (A) = n d u-7ii d (i+7)*n p -0'i'--- .j*) 

is /£A* e 



(169) 



(170) 



ee n v^v^ n ^ 



eeA* 



n 

i>GA* 




Thus in the modified EPRL model edges e G A* are assigned pairs of intertwiner quantum numbers 
l% and if; and an edge amplitude given by the matrix elements g% , c (where v s and Vt stand for 
the source and target vertices of the given oriented edge). The fact that edges are not assigned a 
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single quantum number is not really significative; one could go to a basis of normalized eigenstates 
of and rewrite the modified model above as a spin foam model where edges are assigned a 
single (basis element) quantum number. As the nature of such basis and the quantum geometric 
interpretation of its elements is not clear at this stage, it seems simpler to represent the amplitudes 
of the modified model in the above form. 

The advantages of the modified model are important,; however, a generalization of the above 
modification of the EPRL model in the Lorentzian case is still lacking. Notice that this modification 
does not interfere with the results on the semiclassical limit (to leading order) as reviewed in Section 
If 31 The reason is that the matrix elements g e a a — > 5 a p in that limit PQ . 

7.6 The coherent states representation 

We have written the amplitude defining the EPRL model by constraining the state sum of BF 
theory. For semiclassical studies that we will review in Section [T3] it is convenient to express the 
EPRL amplitude in terms of the coherent states basis. The importance of coherent states in spin 
foam models was put forward in |246j and explicitly used to re-derive the EPRL model in [247] . 
The coherent state technology was used by Freidel and Krasnov in [176] to introduce a new kind 
of spin foam models for gravity: the FK models. In some cases the FK model is equivalent to the 
EPRL model; we will review this in detail in Section [9] 

The coherent state representation of the EPRL model is obtained by replacing (|113|) in each of 
the intermediate SU(2) (green) wires in the expression (1148)) of the EPRL amplitudes, namely 




(171) 



The case 7 < 1 

In this case the coherent state property (|114p implies 
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(172) 

where we used in the last line the fact that for 7 < 1 the representations j of the subgroup 
SU(2) 6 Spin(4) are maximum weight, i.e., j = j + + j~. Doing this at each edge we get 



il /6A* 




(173) 



where we have explicitly written the n S S 2 integration variables on a single cable. The expression 
above is very similar to the coherent states representation of Spin(A) BF theory given in Equation 
(1115|) . In fact one would get the above expression if one would start form the expression (fTT5l) and 
would set n^j = n~j — n e f while dropping for example all the sphere integrations corresponding 
to the ri^f (or equivalently n~j). Moreover, by construction the coherent states participating in 
the previous amplitude satisfy the linear constraints (|127|) in expectation values, namely 

(j,n ef \D}\j,n ef ) = 

= (j,n e/ |(l- 7 )J+' 1 + (l + 7 )J7 J b-,n e/ ) -0. (174) 

Thus the coherent states participating in the above representation of the EPRL amplitudes solve 
the linear simplicity constraints in the usual semiclassical sense. The same manipulations leading 
to (|176|) in Section [5] lead to a discrete effective action for the EPRL model, namely 

KTn^Y. II d (1 - 7 )££ d ( 1+7) £X / II d jBf dn ef dg-dg+ exp (SJ?Jg±]), (175) 
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where the discrete action 



S y ±' to*] = V (S v „ \g-]+S v „ [g + ]) (176) 

u£A* 

with 



S ln[9]= E 2jab\n-(n a b\g a 1 gb\nba), (177) 

a<b=l 

and the indices a, b label the five edges of a given vertex. The previous expression is exactly equal 
to the form of the BF amplitude. In the case of the gravity models presented here, the coherent 
state path integral representation (analogous to (JTTTJ)) will be the basic tool for the study of the 
semiclassical limit of the models and the relationship with Regge discrete formulation of general 
relativity. 



The case 7 > 1 

The case 7 > 1 is more complicated [55]. The reason is that the step (I172j) directly leading to 
the discrete action in the previous case is no longer valid as the representations of the subgroup 
SU{2) G Spin(A) are now minimum instead of maximum weight. However, the representations 
j + = j~ + j are maximum weight. We can therefore insert coherent states resolution of the 
identity on the right representations and get: 




(178) 



where we are representing the relevant part of the diagram appearing in equation (|171[) . In the 
last line we have used that j + = j + j~ (i.e. maximum weight), and the graphical notation 
mo-on = (m\n) as it follows from our previous conventions. With all this, one gets 

ZlZrt = (179) 
II d ( i_ 7 )i/ d (i +7 )i£ / II d Je/d (1+7)£££ dn e/ dm e /d 5 -d ff + exp (S^Jff^), 

Of /£A* 2 2 J eG A* 2 

where the discrete action 

ugA* 
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with 

S ']±,n,m[.9 ± ] = 3ab(l + 7) l og{(mab\gt b \mba)) + jab{l - l)log((m a b\g~ b \mba)) + 

l<a<b<5 

+2j ab {\og((n ab \m ab )) + \og((m ba \n ba ))) . (181) 

7.7 Some additional remarks 

It is important to point out that the commutation relations of basic fields — reflecting the simple 
algebraic structure of spin (4) — used here is the one induced by the canonical analysis of BF 
theory presented previously. The presence of constraints generally modifies canonical commutation 
relations in particular in the presence of second class constraints. For some investigation of the 
issue in the context of the EPRL and FK models see [13] ■ In [TT] it is pointed out that the 
presence of secondary constraints in the canonical analysis of Plebanski action should translated in 
additional constraints in the holonomies of the spin foam models here considered (see also [12]). A 
possible view is that the simplicity constraints are here imposed for all times and thus secondary 
constraints should be imposed automatically. 

There are alternative derivations of the models presented in the previous sections. In particular 
one can derive them from a strict Lagrangean approach of Plebanski's action. Such viewpoint is 
taken in [86] [85] [91]. The path integral formulation of Plebansky theory using commuting in- 
fields was studied in [212] , where it is shown that only in the appropriate semiclassical limit the 
amplitudes coincide with the ones presented in the previous sections (this is just another indication 
that the construction of the models have a certain semiclassical input; see below). The spin foam 
quantization of the Hoist formulation of gravity via cubulations was investigated in |48) . The 
simplicity constraints can also be studied from the perspective of the U(N) formulation of quantum 
geometry [153[ I155[ 1154] which also has been explored in the Lorentzian sector |150j . Such U(N) 
treatment is related to previous work which has been extended to a completely new perspective 
on quantum geometry with possible advantageous features [98] 1249] . For additional discussion on 
the simplicity constraints see [143] . 

8 The Lorentzian EPRL Model: Representation Theory 
and Simplicity Constraints 

In this section we introduce some elements of SL(2, C) representation theory and show how the so 
called linear simplicity constraints are solved in the Lorentzian EPRL model [2941 11611 1158] . 

8.1 Representation theory of SL(2, C) and the canonical basis 

Unitary irreducible representations J4? p .k of SL(2, C) are infinite dimensional and are labelled by 
a positive real number p £ K + and a half-integer k £ N/2. The two Casimirs are C\ = iJijJ = 
L 2 — K 2 and C2 = \*JijJ ij — K ■ L where U are the generators of an arbitrary rotation subgroup 
and K l are the generators of the corresponding boosts. The Casimirs act on \p, k) £ J^ Pi u as 
follows 

d\p,k) = (k 2 -p 2 ~l)\p,k) 

C 2 \p,k) = 2pk\p,k). (182) 

For detail on the representation theory of SL(2,C) see [3311119611197] , As in the Riemannian case, 
the definition of the EPRL model requires the introduction of an (arbitrary) subgroup SU(2) C 
SL(2, C). This subgroup correspond to the internal gauge group of the gravitational phase space 
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in connection variables in the time gauge. In the quantum theory, the representation theory of this 
SU (2) subgroup will be hence important. This importance will soon emerge as apparent from the 
imposition of the constraints that define the EPRL. The link between the unitary representations 
of SL(2, C) and those of SU(2) is given by the decomposition 

oo 

= ©«#?•• ( 183 ) 

As the unitary irreducible representations of the subgroup SU(2) £ SL(2, C) are essential in 
understanding the link of the EPRL model and the operator canonical formulation of LQG it 
will be convenient to express the action of the generators of the Lie algebra s/(2,C) in a basis 
adapted to the above equation. In order to do this we first notice that the Lie algebra (I120j) can be 
equivalently characterized in terms of the generators of the rotation subgroup If and the remaining 
boost generators K % as follows 

[L 3 ,L±]=±L± [L +1 L_]=2L 3 
[L+,K+] = [£_,#_] = [L 3 ,K 3 }=0 

[K 3 , L±] =±K ± [L ± ,K T ] = ±2 K 3 [L 3 , K±] =±K± 

[K 3 ,K ± ] = T L ± [K + ,K_}=-2 L 3 , (184) 

where K± = K 1 ± iK 2 and L± = L x ±iL 2 respectively. The action of the previous generators in 
the basis \p, k, j, m) can be shown to be [197J 

L 3 \p,k,j,m) = m\p,k,j,m), 

L + \p, k, j, m) = y/(J + m+l)(j - m)\p, k,j,m + 1), 
L~\p,k,j,m) = \J (j + m)(j -m + l)\p, k, j, m - 1), 



K 3 \p,k,j,m) = QijVi 2 -m 2 \p,k, j - l,m) +^ j m\p,k,j,m) - a j+1 y/(j + l) 2 - m 2 |p, k,j + l,m), 

K + \p,k,j,m) = a 3 \J {j - m) (j - m - 1) \p, k, j - l,m+l) 

+7jV(j -m)(j + m + l)\p, k,j, m + 1) 

+a j+ly / (j + m + + m + 2)\p, k,j + l,m+ 1), 

K~\p,k,j,m) = -a jy /(j + m)(j + m - l)\p,k,j -l,m- 1) 

+7j V U + m)(j -m + l)\p, k, j, m - 1) 

-a j+iy /(j - m + 1)0" - m + 2)b. k,j + l,m-l), (185) 
where 

kp . (j 2 - k 2 )(j 2 +p 2 ) 

"^WTT) a * = i iW-D (186) 

This concludes the review of the representation theory that is necessary for the definition of the 
EPRL model. 

8.2 Lorentzian model: the linear simplicity constraints 

As in the Riemannian case it can be shown as in Section (17.41) that the quadratic simplicity 
constraints (j!54[) are implied, in the discrete setting, by the linear constraint on each face 

D) = L) - -K) w 0, (187) 
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where the label / makes reference to a face / € A*, and where (very importantly) the subgroup 
SU(2) C SL(2, C) that is necessary for the definition of the above constraints is chosen arbitrarily 
at each tetrahedron, equivalent on each edge e £ A*. The commutation relations of the previous 
tetrahedron constraint are 



[D},D},} = 5 f f> 



(1 - - 



2,, ; , e \D k -S eel ^^-e%L k f . il.ssi 



The previous commutation relation implies that the constraint algebra is not closed and cannot 
therefore be imposed as operator equations of the states summed over in the BF partition function 
in general. There are two interesting exceptions: 

1. The first one is to take 7 = ±z in which case the constraint (|187p could in principle be 
imposed strongly. The self-dual sector is however not present in J4? p ,k- 

2. The second possibility is to work in the sector where L\ = 0. Which, as in the Riemannian 
case, it corresponds to the Barret-Crane model [53] with the various limitations discussed in 
Section ]W\ 

The EPRL model is obtained by restricting the representations appearing in the expression of 
the BF partition function so that at each tetrahedron the linear constraints (j!87[) the strongest 
possible way that is compatible with the uncertainties relations stemming from (|188p . In addition 
one would add the requirement that the state-space of tetrahedra is compatible with the state- 
space of the analogous excitation in the canonical context of LQG so that arbitrary states in the 
kinematical state of LQG have non trivial amplitudes in the model. 



On the weak imposition of the linear simplicity constraints 

We now discuss the weak imposition of the linear simplicity constraints in the quantum framework. 
There are essentially three views in the literature: two of them, discussed below, concerns directly 
the way the EPRL model has been usually introduced. The third possibility is the semiclassical 
view based on the coherent state representation leading to the FK model (see section 



8.2.1 Lorentzian model: The Gupta-Bleuler criterion 

As in the Riemannian case, the fact that the constraints D l j do not form a closed (first class) algebra 
in the generic case requires a weaker sense in which the constraints are imposed. One possibility is 
to consider the Gupta-Bleuler criterion consisting in selecting a suitable class of states for which the 
matrix elements on vanish. In this respect one notices that if we chose the subspace Jifj C J%? Pt k 
one has 

(p,k,j,q\D f \p,k,j,m) = 5 q , m m(l -) 



(p,k,j,q\Df\p,k,j,m) = S q±hm ^/(j ±m+ =pm)(l - — ). 

One immediately observes that matrix elements of the linear constraints vanish in this subclass if 
one can choose 

*-s£iT 7 (189) 

Now the situation simplifies with respect to the Riemannian case as there is only one single case. 
The linear simplicity constraint is satisfied as long as p — 7j(j + !)/&• Therefore, one has an 
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integer worth of possibilities. In particular, given j one has k = j — r for all possible half-integers 
such that k > [134] . As we will see in what follows the best choice from the point of view of the 
Master constraint criterion is k — j from which it follows 

p = 1 { ] + l). (190) 

8.2.2 Lorentzian model: The Master constraint criterion 

Another criterion for weak imposition can be developed by studying the spectrum of the Master 
constraint 

M f = DfDf = L 2 — — 2L ■ K + —r;K 2 

= ^[{l+l 2 )L 2 -2 1 C 2 -C 1 }. (191) 

7 

Strong imposition of the constraints Dy would amount to looking for the kernel of the master 
constraint Mf. However, generically the positive operator associated with the master constraint 
does not have the zero eigenvalue in the spectrum due to the open nature of the constraint algebra 
(|188p . The proposal of [160] is to look for the minimum eigenvalue among spaces J&j € J%>,k- 
Explicitly 

M f \ip >=m p<k ,i\i> >, (192) 

where 

m p ,k,j = ^b 2 + 1 - k 2 + (1 + 7 2 )i(j + 1) - 2-ykp]. 

For given j the minimum eigenvalue is obtained when k = j and p = jj in which case one gets 
m P,kj > mj for rrij = h 2 (l + r ) 2 )j + 1. Again we see that the important condition 

M f \-(j,j,j,rn) = ^ c (193) 

holds (recall the notion G sc from Section f7.2.4|) . Interestingly, the states that satisfy the previous 
semiclassical criterion satisfy the Gupta-Bleuler criterion not exactly and only in the large spin 
sense. However, the advantage of the master constraint criterion is clear from the point of view of 
cylindrical consistency as the choice p = 7(5 + 1) obtained in the previous subsection would imply 
that ,7=0 spin network links would be have a non trivial amplitude in the path integral defined in 
Section [7.31 Such thing would be in sharp conflict with diffeomorphism invariance in the model. 

Therefore, as in the Riemannian case we will use the master constraint criterion as the basic 
ingredient in the definition of the EPRL model. One can write a Lorentz invariant constraint [326 
from the master constraint and using the = classically. This leads to the following Lorentz 
invariant condition in terms of the SL(2,C) Casimirs 

M} 1 = (1 - 7 2 )C 2 - 2C l7 , (194) 

where C\ and C2 are the two si (2, C) Casimirs given in equation (|182[) . Surprisingly, the conditions 
p = jj and k — j minimizing the master constraint solve its Lorentz invariant version strongly! 
This resolves the cylindrical consistency issue pointed out above. More recently, it has been shown 
that the constraint solution p — jj and k = j also follows naturally from a spinor formulation of 
the simplicity constraints 363 , 150, 248 . 

8.2.3 Lorentzian model: Restrictions on the Immirzi parameter 

Unlike the Riemannian model, all values of spin j are allowed, and no restrictions on the Immirzi 
parameter arise. 



66 



8.2.4 Lorentzian model: Overview of the solutions of the linear simplicity constraints 

In the Lorentzian model one restricts the SL(2, C) representations of BF theory to those satisfying 

P = 13 k=j (195) 
for j £ N/2. From now on we denote the subset of admissible representation of SL(2, C) as 

Jt-y C Irrep(SL(2,C)) (196) 

The admissible quantum states ^ are elements of the subspace J#j C J^yjj (i-e., minimum weight 
states) satisfy the constraints (|187l) in the following emiclassical sense: 

(K) - 1 L))^ = e sci (197) 

where as in the Riemannian case the symbol G sc (order semiclassical) denotes a quantity that 
vanishes in limit h — > 0, j — > oo with hj ^constant. 

8.3 Presentation of the EPRL Lorentzian model 

As briefly discussed in Section HI unitary irreducible representations of SL(2,C) are infinite di- 
mensional and labelled by a positive real number p £ R + and a half-integer k £ N/2. These 
representation are the ones that intervene in the harmonic analysis of square integrable functions 
of SL(2,C) |196j . Consequently, one has an explicit expression of the delta function distribution 
(defined on such test function), namely 

6(g) = W dp (p 2 + fc 2 ) ]T DV&M (198) 

k jR+ j,m 

where -Dj^'m' (flO w ith j > k and j > m > —j (similarly for the primed indices) are the matrix 
elements of the unitary representations p — k in the so-called canonical basis [331] . One can use 
the previous expression the Lorentzian version of Equation (|86[) in order to introduce a formal 
definition of the BF amplitudes, which now would involve integration of the continuous labels pf 
in addition of sums over discrete quantum numbers such as k, j and m. The Lorentzian version of 
the EPRL model can be obtained from the imposition of the linear simplicity constraints, discussed 
in Section \8?2\ to this formal expression. As the continuum labels pf are restricted to pf = yjf 
the Lorentzian EPRL model becomes a state-sum model as its Riemannian relative. Using the 
following graphical notation 

D P j £ j , m ,{9)= >> m ■ I ■ (199) 

the amplitude is 



J eprl 
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where the boxes now represent SL{2, C) integrations with the invariant measure. The previous 
amplitude is equivalent to the its spin foam representation 




The vertex amplitude is divergent due to the presence of a redundant integration over SX(2, C), it 
becomes finite by dropping an arbitrary integration, i.e. removing any of the 5 boxes in the vertex 
expression [161J . 



8.4 The coherent state representation 

It is immediate to obtain the coherent states representation of the Lorentzian models. As in the 
Riemannian case, one simply inserts resolution of the identities (1108)) on the intermediate SU (2) 
(green) wires in (|200l) from where it results 




(200) 



9 The Freidel— Krasnov Model 

Shortly after the appearance of the paper 162 , Freidel and Krasnov [176] introduced a set of 
new spin foam models for four dimensional gravity using the coherent state basis of the quantum 
tetrahedron of Livine and Speziale [246] . The idea is to impose the linearized simplicity constraints 
(|126p directly as a semiclassical condition on the coherent state basis. As we have seen above, 
coherent states are quantum states of the right and left tetrahedra in BF theory which have a 
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clear-cut semiclassical interpretation through their property (|112j) . We have also seen that the 
imposition of the linear constraints (|126|) a la EPRL is in essence semiclassical as they are strictly 
valid only in the large spin limit. In the FK approach one simply accept from the starting point 
that, due to their property of not defining set that is closed under commutation relations, the 
Plebansky are to be imposed semiclassically. One defines new models by restricting the set of 
coherent states entering in the coherent state representation of Spin(4) BF theory (|115j) to those 
that satisfy condition (|126j) in expectation values. They also emphasize how the model [162 
corresponds indeed to the sector 7 = 00 which has been shown to be topological |240] . 

The case 7 < 1 

For 7 < 1 the vertex amplitude is identical to the EPRL model. This is apparent in the coherent 
state expression of the EPRL model (|175j) . Thus we have 

^;r(A) = E II d u-^ d u + 7)i 

is feA* 



n 



d,-,, wd, 
(1+7)5 (7- 



dn 



ef 




|| (201) 



From the previous expression we conclude that the vertex amplitudes of the FK and EPRL model 
coincide for 7 < 1 

^JT/fc = ^eprf (202) 

Notice however that different weights are assigned to edges in the FK model. This is due to 
the fact that one is restricting the Spin(A) resolution of identity in the coherent basis in the 
previous expression, while in the EPRL model the coherent state resolution of the identity is 
used for SU (2) representations. This difference is important and has to do with the still un- 
settled discussion concerning the measure in the path integral representation (see Section [1 2. II and 
references therein). 

The case 7 > 1 

For the case 7 > 1 the FK amplitude is given by 
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The study of the coherent state representation of the FK model for 7 > 1, and comparison with 
equation (|178|l for the EPRL model, clearly shows the difference between the two models in this 
regime. 

Z h = H II d (1 - 7) ^ d ( 1+7) ^ / II d ll _ 7l¥ d {1+ ^dn ef d 9 -dgt v exp^i;^]), (204) 

j f feA* J eG A* 

where the discrete action 

S' s l VJ?] = E K ^ + ^ r F 5 + ]), (205) 

where 5(7) = sign(l — 7) and 

5 

Slntff] = E 2 Jafc ln ( n afclfa 1 5fc|«ha). ( 206 ) 
a<6=l 

with the indices a, & labelling the five edges of a given vertex. 

10 The Barrett-Crane Model 

The Barrett-Crane models [531 H22 can be easily obtained in the present framework as a suitable 
limiting case 7 — > 00 of the previous constructions. One way to get it is to take the limit 7 — > 00 
of the linear simplicity constraint (I126[) at the classical level. In this way the constraint becomes 

D i =L}«0 (207) 

with which the constraint algebra closes. One can impose the constraints strongly in this case. 
The solution is obtained by setting jf — in all faces at the level of the amplitude (|148|) . An 
immediate consequence of this is that jf = jj and all tetrahedra are assigned a unique intertwiner: 
the so-called Barrett-Crane intertwiner. Explicitly, setting jt = in (jl48[) we get 




(208) 
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where since jf — jj we are arbitrarily choosing jf as the spins summed over in the state sum. 
Eliminating redundant group integrations we get the previous amplitude in the spin foam repre- 
sentation, namely 




(209) 



Observe that the vertex amplitude depend only on the 10 corresponding jj, t his is the Riemannian 
Barrett-Crane vertex. It is interesting to write the analog of equation (I147[) . The Barrett-Crane 
model is obtained by the modification of the BF projectors in (l^5|) so that 

Z BC (A) = II d Ml P ^c(Ji + , • • • Jt) (210) 

3+ /£A' 

with 




P e BC {h---h) = \BC){BC\= = EE l' ) ( 211 ) 




where \BC) is the Barrett-Crane intertwiner. The tetrahedra state space is one-dimensional once 
the four quantum numbers jf has been fixed. Indeed, it has been shown that the state space of 
strong solutions of the the quadratic quantum Plebanski constraints is one dimensional and it is 
given by the BC intertwiner \BC) [307]. 



10.1 The coherent states representation of the BC amplitude 

One can insert coherent states resolution of the identity (jlllj) in the 10 wires obtaining the following 
coherent state representation of the Barrett-Crane amplitude 

Z bc = J2 II d "+ / II *i.fdn ef dg+ exp(^ n [ 5 +]), (212) 
where the discrete action 

Sf°n\9 + ]='E E 2j+l n <r la6 |( 5 +)- 1 9 +|n ab ). (213) 

»£A l a<b=l 

Notice that there are only 10 normals n a b as in this case n a \, = nt, a . The large spin asymptotics of 
the Barrett-Crane vertex has been investigated in [53] • 



71 



11 Boundary Data for the New Models and Relationship 
with the Canonical Theory 

So far we have considered cellular complexes with no boundary. Transition amplitudes are expected 
to be related to the definition of the physical scalar product as discussed in Section 12.21 In order 
to define them one needs to consider complexes with boundaries. Boundary states are defined on 
the boundary of the dual two-complex A* that we denote OA* . The object dA* is a one-complex 
(a graph). According to the construction of the model (Section [7] and [HJ boundary states are in 
one-to-one correspondence with SU{2) spin networks. This comes simply from the fact that links 
(one-cells) I € dA* inherit the spins labels (unitary irreducible representations of the subgroup 
SU(2)) of the boundary faces while nodes (zero-cells) n g dA* inherit the intertwiner levels of 
boundary edges. 

At this stage one can associate the boundary data with elements of a Hilbert space. Being 
in one-to-one correspondence with SU(2) spin networks, a natural possibility is to associate to 
them an element of the kinematical Hilbert space of LQG. More precisely, with a given coloured 
boundary graph 7 with links labelled by spins je and nodes labelled by interwiners i n we associate 
a cylindrical function \& j,{j e },{i. n } € Jf 2 (SU(2) Ne ), where here Ni denotes number of links in 
the graph 7 (see Equation 169]) . In this way, the boundary Hilbert space associated with dA* is 
isomorphic (if one used the natural AL measure) with the Hilbert space of LQG truncated to 
that fixed graph. Moreover, geometric operators such as volume and area defined in the covariant 
context are shown to coincide with the corresponding operators defined in the canonical formulation 
|135[ 1136] . Now, if cellular complexes are dual to triangulations then the boundary spin networks 
can have at most four valent nodes. This limitation can be easily overcome: as in BF theory 
the EPRL amplitudes can be generalized to arbitrary complexes with boundaries given by graphs 
with nodes of arbitrary valence. The extension of the model to arbitrary complexes has been first 
studied in [2311 1230] , it has also been revisited in |134] . 

Alternatively, one can associate the boundary states with elements of Jz? 2 (Spin(4) Ne ) (in the 
Riemannian models) - or carefully define the analog of spin network states as distributions in the 
Lorentzian casd^l. In this case one gets special kind of spin network states that are a subclass of the 
so-called projected spin networks introduced in [TH] 1241] in order to define an heuristic quantization 
of the (non-commutative and very complicated) Dirac algebra of a Lorentz connection formulation 
of the phase space of gravity (T4] [10] [16] [9] [8] [7] [15]. The fact that these special subclass of 
projected spin networks appear naturally as boundary states of the new spin foams is shown in 

Due to their similarity for 7 < 1 the same relationship between boundary data and elements of 
the kinematical Hilbert space hold for the FK model. However, the such simple relationship does 
not hold for the model in the case 7 > 1. 

The role of knotting in the Hilbert space of the EPRL model is studied in [38] . 

25 The definition of a polymer-like gauge invariant states as elements of a Hilbert space of a gauge theory with 
non compact group is an open issue. The naive attempts fail, the basic problem is that gauge invariant states are 
not square intcgrable functions with respect to the obvious generalization of the AL measure for fixed graph: this 
is due to the infinite volume of the gauge group (the invariant measure is not normalizable). For that reason the 
discussion of the Lorentzian sector with G = SL(2, C) is usually formal. See |177l for some insights on the problem 
of defining a gauge invariant Hilbert space of graphs for non compact gauge groups. 
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12 Further developments and related models 



12.1 The measure 

The spin foam amplitudes discussed in the previous sections have been introduced by constraining 
the BF histories through the simplicity constraints. However, in the path integral formulation, the 
presence of constraints has the additional effect of modifying the weights with which those histories 
are to be summed: second class constraints modify the path integral measure (in the spin foam 
context this issue was raised in [84 ). As pointed out before, this question has not been completely 
settled yet in the spin foam community. The explicit modification of the formal measure in terms 
of continuous variables for the Plebansky action was presented in |100j . A systematic investigation 
of the measure in the spin foam context was attempted in [157] and |208j . As pointed out in [84] , 
there are restrictions in the manifold of possibilities coming from the requirement of background 
independence. The simple BF measure chosen in the presentation of the amplitudes in the previous 
sections satisfy these requirements. There are other consistent possibilities; see for instance |72j for 
a modified measure which remains extremely simple and is suggested from the structure of LQG. 

12.2 Relation with the canonical formulation: the scalar constraint 

An important question is the relationship between the spin foam amplitudes and the canonical 
operator formulation. The question of wether one can reconstruct the Hamiltonian constraints 
out of spin foam amplitudes has been analysed in detail in three dimensions. For the study of 
quantum three dimensional gravity from the BF perspective see [275] . we will in fact present this 
perspective in detail in the three dimensional part of this article. For the relationship with the 
canonical theory using variables that are natural from the Regge gravity perspective see [55J [HI] ■ 
In four dimensions the question has been investigated in [3] in the context of the new spin foam 
models. In the context of group field theories this issue is explored in |244j . Finally, spin foams can 
in principle be obtained directly from the implementation of the Dirac program using path integral 
methods this avenue has been explored in [2131 1214] from which a discrete path integral formulation 
followed [209 . The question of the relationship between covariant and canonical formulations in 
the discrete setting has been analyzed also in |142] , 

12.3 Spin foam models with timelike faces 

By construction all tetrahedra in the FK and EPRL models are embedded in a spacelike hypersur- 
face and hence have only spacelike triangles. It seem natural to ask the question of whether a more 
general construction allowing for timelike faces is possible. The models described in previous sec- 
tions have been generalized in order to include timelike faces in the work of Conrady |119IU14llll8] . 
An earlier attempt to define such models in the context of the Barrett-Crane model can be found 
in [299] . 

12.4 Coupling to matter 

The issue of the coupling of the new spin foam models to matter remains to a large extend 
un-explored territory. Nevertheless some results can be found in the literature. The coupling 
of the Barrett-Crane model (the 7 — > 00 limit of the EPRL model) to Yang-Mills fields was 
studied in |288j . More recently the coupling of the EPRL model to fermions has been investigated 
in [2111 [67] . A novel possibility of unification of the gravitational and gauge fields was recently 
proposed in [6]. 
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12.5 Cosmological constant 

The introduction of a cosmological constant in the construction of four dimensional spin foam 
models has a long history. Barrett and Crane introduced a vertex amplitude [54] in terms of the 
Crane and Yetter model [124] for BF theory with cosmological constant. The Lorentzian quantum 
deformed version of the previous model was studied in |277j . For the new models the coupling 
with a cosmological constant is explored in terms of the quantum deformation of the internal gauge 
symmetry in [1331 12071 1210] as well as (independently) in }168j . 

12.6 Spin foam cosmology 

The spin foam approach applied to quantum cosmology has been explored in 68, 360, 216, 73, 
328 , 327 . The spin foam formulation can also be obtained from the canonical picture provided by 
loop quantum cosmology (see [83] and references therein). This has been explored systematically 
in [25] El [23] [102]. Very recent results are [250] 

12.7 Constraints are imposed semiclassically 

As we have discussed in the introduction of the new models, Heisenberg uncertainty principle 
precludes the strong imposition of the Plebanski constraints that reduce BF theory to general 
relativity. The results on the semiclassical limit of these models seem to indicate that metric 
gravity should be recovered in the low energy limit. However, its seems likely that the semiclassical 
limit could be related to certain modifications of Plebanski's formulation of gravity [234 , 235 , 1237[ 
12361 [238 . A simple interpretation of the new models in the context of the bi-gravity paradigm 
proposed in |342] could be of interest. 

12.8 Group field theories associated to the new spin foams 

As already pointed out in [34] spin foams can be interpreted in close analogy to Feynman diagrams. 
Standard Feynman graphs are generalized to 2-complexes and the labeling of propagators by mo- 
menta to the assignment of spins to faces. Finally, momentum conservation at vertices in standard 
feynmanology is now represented by spin-conservation at edges, ensured by the assignment of the 
corresponding intertwiners. In spin foam models the non-trivial content of amplitudes is contained 
in the vertex amplitude which in the language of Feynman diagrams can be interpreted as an 
interaction. This analogy is indeed realized in the formulation of spin foam models in terms of a 
group field theory (GFT) [30911308] . 

The GFT formulation resolves by definition the two fundamental conceptual problems of the 
spin foam approach: diffcomorphism gauge symmetry and discretization dependence. The difficul- 
ties are shifted to the question of the physical role of A and the convergence of the corresponding 
perturbative series. 

In three dimensions this idea has been studied in more detail. In [256] scaling properties of 
the modification of the Boulatov group field theory introduced in [181] was studied in detail. 
In a further modification of the previous model (known as coloured tensor models [202]) new 
techniques based on a suitable I/TV expansion imply that amplitudes are dominated by spherical 
topology [201] ; moreover, it seem possible that the continuum limit might be critical as in certain 
matrix models |204[ 189] I203| I205| 1332] . However, it is not yet clear if there is a sense in which 
these models correspond to a physical theory. The naive interpretation of the models is that they 
correspond to a formulation of 3d quantum gravity including a dynamical topology. 
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12.9 Results on the semiclassical limit of EPRL-FK models 



13 The semiclassical limit 

Having introduced the relevant spin foam models in the previous sections we now present the results 
on the large spin asymptotics of the spin foam amplitudes suggesting that on a fixed discretization 
the semiclassical limit of the EPRL-FK models is given by Regge's discrete formulation of general 
relativity [55 1 155 ] . 

The semiclassical limit of spin foams is based on the study of the the large spin limit asymptotic 
behaviour of coherent state spin foam amplitudes. The notion of large spin can be defined by the 
rescaling of quantum numbers and Planck constant according to j — > Xj and h — > h/X and taking 
A >> 1. In this limit the quantum geometry approximates the classical one when tested with 
suitable states (e.g. coherent states). However, the geometry remains discrete during this limiting 
process as the limit is taken on a fixed regulating cellular structure. That is why one usually 
makes a clear distinction between semiclassical limit and the continuum limit. In the semiclassical 
analysis presented here one can only hope to make contact with discrete formulations of classical 
gravity; hence the importance of Regge calculus in the discussion of this section. 

The key technical ingredient in this analysis is the representation of spin foam amplitudes in 
terms of the coherent state basis introduced in Section 16.21 Here we follow [55J [55] [55] [5j|] [57] . 
The idea of using coherent states and discrete effective actions for the study of the large spin 
asymptotics of spin foam amplitudes was put forward in [1151 1116] . The study of the large spin 
asymptotics has a long tradition in the context of quantum gravity dating back to the studied of 
Ponzano-Regge [303] . More directly related to our discussion here are the early works [53] 155] . 
The key idea is to use asymptotic stationary phase methods for the amplitudes written in terms 
of the discrete actions presented in the previous section. 

13.1 Vertex amplitudes asymptotics 

In this section we review the results of the analysis of the large spin asymptotics of the EPRL 
vertex amplitude for both the Riemannian and Lorentztian models. We follow the notation and 
terminology of [55] and related papers. 

13.1.1 SU(2) 15j-symbol asymptotics 

As SU (2) BF theory is quite relevant for the construction of the EPRL-FK models, the study of 
the large spin asymptotics of the SU (2) vertex amplitude is a key ingredient in the analysis of [55] . 
The coherent state vertex amplitude is 



which depends on 10 spins j a b and 20 normals n a i, ^ rib a . The previous amplitude can be expressed 
as 




(214) 




o=l 
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(215) 




(216) 



a<b=l 
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and the indices a, b label the five edges of a given vertex. The previous expression is exactly equal 
to the form (|86p of the BF amplitude. In the case of the EPRL model studied in Sections [7] and 
[51 the coherent state representation - see equations 11751 11791 and 12001 - is the basic tool for the 
study of the semiclassical limit of the models and the relationship with Regge discrete formulation 
of general relativity. 

In order to study the asymptotics of (|215j) one needs to use extended stationary phase methods 
due to the fact the the action (|216|) is complex (see |115| 1116] ). The basic idea is that in addition 
to stationarity one requires real part of the action to be maximal. Points satisfying these two 
conditions are called critical points. As the real part of the action is negative definite, the action 
at critical points is purely imaginary. 

Notice that the action (|216|) depends parametrically on 10 spins j and 20 normals n. These 
parameters define the so-called boundary data for the four simplex v G A*. Thus, there is an 
action principle for every given boundary data. The number of critical points and their properties 
depend on these boundary data, hence the asymptotics of the vertex amplitude is a function of 
the boundary data. Different cases are studied in detail in |55j , here we present their results in the 
special case where the boundary data describe a non-degenerate Regge geometry for the boundary 
of a four simplex, these data are referred to as Regge-like, and satisfy the gluing constraints. For 
such boundary data the action (|216[) has exactly two critical points leading to the asymptotic 
formula 

' 7V + exp(z^Aj Qb ef b )+iV_exp(- J ^Aj Qb ef ; 6 ) , (217) 

a<b a<ib 



15j(Aj,n) ^ 



where ab the appropriate diahedral angles defined by the four simplex geometry; finally the N± 
are constants that do not scale with A. 



13.1.2 The Riemannian EPRL vertex asymptotics 

The previous result together with the fact that the EPRL amplitude for 7 < 1 is a product of 
SU (2) amplitudes with the same n in the coherent state representation (|175l) implies the asymptotic 
formula for the vertex amplitude to be given by the unbalanced square of the above formula [56] . 
namely 
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One can write the previous expression as 
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where 



a<b 



(219) 



is the Regge like action for Xjj a b — A a b the ten triangle areas (according to the LQG area spectrum 
l|7ip). Remarkably, the above asymptotic formula is also valid for the case 7 > 1 [55]. The first 
term in the vertex asymptotics is in essence the expected one: it is the analog of the 6j symbol 
asymptotics in three dimensional spin foams. Due to their explicit dependence on the Immirzi 
parameter, the last two terms are somewhat strange from the point of view of the continuum field 
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theoretical view point. However, this seems to be a peculiarity of the Riemannian theory alone as 
the results of [55] for the Lorentzian models show. Non geometric configurations are exponentially 
surpressed. Finally, recently is has been shown that it is possible to restrict the state sum in order 
to obtain complex amplitudes |156j and [2651 1266] . 



13.1.3 Lorentzian EPRL model 

To each solution one can associate a second solution corresponding to a parity related 4-simplex 
and, consequently, the asymptotic formula has two terms. It is given, up to a global sign, by the 
expression 
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where N± are constants that do not scale. Non geometric configurations are exponentially sur- 
pressed. By taking into account discrete symmetries it might be possible to further improve the 
asymtotics of the model [329] , 



13.2 Full spin foam amplitudes asymptotics 

In [116] Freidel and Conrady gave a detailed description of the coherent state representation of 
the various spin foam models described so far. In particular they provided the definition of the 
effective discrete actions associated to each case which we presented in (|204[) . This provides the 
basic elements for setting up the asymptotic analysis presented in [115] (the first results on the 
semiclassical limit of the new spin foam models) which is similar to the studies of the asymptotic of 
the vertex amplitude reviewed above but more general in the sense that the semiclassical limit of 
a full spin foam configuration (involving many vertices) is studied. The result is technically more 
complex as one studies now critical points of the action associated to a coloured complex which in 
addition of depending on group variables g it depends on the coherent state parameters n. The 
authors of [115] write Equation (|204l) in the following way 

Z h=H II (221) 

is /6A* 

where 

WZ*Vf)= I II d |i_ 7 |^ d (1+7) ^ dn e/^ e 7^+ exp(^;[ 5 ± ]). (222) 

They show that those solutions of the equations of motion of the effective discrete action that are 
non geometric (i.e. the contrary of Regge like) are not critical and hence exponentially suppressed 
in the scaling jf — > Xjt with A >> 1. If configurations are geometric (i.e. Regge like) one has 
two kind of contributions to the amplitude assymptotics: those coming from degenerate and non- 
degenerate configurations. If one (by hand) restricts to the non-degenerate configurations then one 
has 

W l^f) ~ xmn Jn^n f) eM^SL 3 JA\ Jf )), (223) 

where n e , n v , and n/ denote the number of edges, vertices, and faces in the two complex A* 
respectively. For more recent results on asymptotics of spin foam amplitudes of arbitrary simplicial 
complexes see [215] . 
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13.3 Fluctuations: two point functions. 

The problem of computing the two point function and higher correlation functions in the context 
of spin foam has received lots of attention recently. The framework for the definition of the 
correlation functions in the background independent setting has been generally discussed by Rovelli 
in [316] and correspods to a special application of a more general proposal investigated by Oeckl 
|2S3|2Sl||2S3|n3[miliniiniBI2HI]- It was then applied to the Barrett-Crane model in HIE], 
where it was discovered that certain components of the two point function could not yield the 
expected result compatible with Regge gravity in the semiclassical limit. This was used as the 
main motivation of the weakening of the imposition of the Plebanski constraints leading to the 
new models. Soon thereafter it was argued that the difficulties of the Barrett-Crane model where 
indeed absent in the EPRL model [2] . The two point function for the EPRL model was calculated 
in [55] and it was shown to produce a result in agreement with that of Regge calculus [75 1254] in 
the limit 7 — > 0. 

The fact that, for the new model, the double scaling limit 7 — > and j — > 00 with 7j=constant 
defines the appropriate regime where the fluctuation behave as in Regge gravity (in the leading 
order) has been further clarified in |255j . This indicates that the quantum fluctuations in the new 
models are more general than simply metric fluctuations. The fact the the new models are not 
metric at all scales should not be surprising as we know that the Plebanski constraints that produce 
metric general relativity out of BF theory has been implemented only semiclassically (in the large 
spin limit). At the deep Planckian regime fluctuations are more general than metric. However, it 
not clear at this stage why this is controlled by the Immirzi parameter. 

All the previous calculations involve a complex with a single four-simplex. The first computation 
involving more than one simplex was performed in [2571171] for the case of the Barrett-Crane model. 
Certain peculiar properties were found and it is not clear at this stage whether these issues remain 
in the EPRL model. Higher order correlation functions have been computed in [330] . the results 
are in agreement with Regge gravity in the 7 — > limit. 



78 



Part IV 

Three dimensional gravity 

In this part we shall review the quantization of three dimensional gravity introducing in this simpler 
context what in this case can be called the spin foam representation of the quantum dynamics. We 
present this case in great detail as it is completely solvable. More precisely, it allows for an explicit 
realization of the path integral quantization for a generally covariant system as generally discussed 
in Section [2.2l Moreover, the path integral takes the form of a sum of spin foam amplitudes which 
share all the kinematical quantum geometric properties with their four dimensional relatives. 

The simplicity of this theory allows for the clear-cut illustration of some central conceptual 
difficulties associated with the quantization of a generally covariant system. We will use the 
insights provided by 3d gravity to support the discussion of certain difficult conceptual issues for 
the spin foam approach in Section [V] 
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14 Spin Foams for 3-Dimensional Gravity 



Three dimensional gravity is an example of BF theory for which the spin foam approach can be 
implemented in a rather simple way. Despite its simplicity the theory allows for the study of many 
of the conceptual issues to be addressed in four dimensions. In addition spin foams for BF theory 
are the basic building block of 4-dimensional gravity models. For a beautiful presentation of BF 
theory and its relation to spin foams see [35] . For simplicity we study the Riemannian theory; the 
Lorentzian generalization of the results of this section have been studied in [172] . 

14.1 Spin foams in 3d quantum gravity 

Here we derive the spin foam representation of LQG in a simple solvable example: 2+1 gravity. 
For the definition other approaches to 3d quantum gravity see Carlip's book [106] , 

14.2 The classical theory 

Riemannian gravity in 3 dimensions is a theory with no local degrees of freedom, i.e., a topological 
theory. Its action (in the first order formalism) is given by 

S[e,w] = J Tr(e A F(u})), (224) 

M 

where M = E X K (for E an arbitrary Riemann surface), cj is an SC/(2)-connection and the 
triad e is an su(2)-valued 1-form. The gauge symmetries of the action are the local SU(2) gauge 
transformations 

8e=[e,a], Suj — d^a, (225) 

where a is a sw(2)-valued 0-form, and the 'topological' gauge transformation 

be = d„i], 5lo = 0, (226) 

where d„ denotes the covariant exterior derivative and r) is a su(2)-valued 0-form. The first 
invariance is manifest from the form of the action, while the second is a consequence of the Bianchi 
identity, d^F(uj) = 0. The gauge symmetries are so large that all the solutions to the equations of 
motion are locally pure gauge. The theory has only global or topological degrees of freedom. 

Upon the standard 2+1 decomposition, the phase space in these variables is parametrized by the 
pull back to E of to and e. In local coordinates one can express them in terms of the 2-dimensional 
connection A l a and the triad field E!j = e bc e^Sjk where a = 1,2 are space coordinate indices and 
i,j = 1, 2, 3 are su(2) indices. The Poisson bracket is given by 

{Ai(x),E b j (y)}=5^5 i j S^(x,y). (227) 

Local symmetries of the theory are generated by the first class constraints 

D bEj - 0, F l ab (A) = 0, (228) 

which are referred to as the Gauss law and the curvature constraint respectively. This simple 
theory has been quantized in various ways in the literature [106] . here we will use it to introduce 
the spin foam representation. 
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14.3 Spin foams from the Hamiltonian formulation 



The physical Hilbert space, H p h ys , is denned by those 'states in %un that are annihilated by 
the constraints. As discussed in |315l 1353] . spin network states solve the Gauss constraint - 

D a Ef\s) = - as they are manifestly SU(2) gauge invariant. To complete the quantization one 
needs to characterize the space of solutions of the quantum curvature constraints F* b , and to 
provide it with the physical inner product. As discussed in Section l2.2i we can achieve this if we 
can make sense of the following formal expression for the generalized projection operator P: 

P = I D[N] exp(i I Tr[NF{A)]) = ]J S[F{A)], (229) 

J £ ICS 

where N(x) G su(2). Notice that this is just the field theoretical analog of equation ([5]). P will 
be defined below by its action on a dense subset of test-states called the cylindrical functions 
Cyl C Hkin ■ If P exists then we have 

(sPU[N], s') = (sP, s') V s,s' G Cyl, N(x) G su(2) (230) 

where U[N] = exp(i J Tr[iNF(A)]). P can be viewed as a map P : Cyl — > Kp C Cyl* (the space 
of linear functionals of Cyl) where Kp denotes the kernel of the curvature constraint. The physical 
inner product is defined as 

(s',s) p :=(s'P,s), (231) 
where (, ) is the inner product in T-Lkim and the physical Hilbert space as 



H phys := Cyl/ J for J := {s G Cyl s.t. (s,s) p = 0}, (232) 

where the bar denotes the standard Cauchy completion of the quotient space in the physical norm. 

One can make (|229p a rigorous definition if one introduces a regularization. A regularization is 
necessary to avoid the naive UV divergences that appear in QFT when one quantizes non-linear 
expressions of the canonical fields such as F{A) in this case (or those representing interactions in 
standard particle physics). A rigorous quantization is achieved if the regulator can be removed 
without the appearance of infinities, and if the number of ambiguities appearing in this process is 
under control. We shall see that all this can be done in the simple toy example of this section. 




Figure 9: Cellular decomposition of the space manifold E (a square lattice of size e in this example), 
and the infinitesimal plaquette holonomy W P [A]. 

We now introduce the regularization. Given a partition of E in terms of 2-dimensional plaque- 
ttes of coordinate area e 2 (Figure [9]) one can write the integral 

F[N] := / Tt[NF(A)} = lim ^VTr[7V p F p ] (233) 

as a limit of a Riemann sum, where N p and F p are values of the smearing field N and the curvature 
e ab Fl b [A] at some interior point of the plaquette p and e ab is the Levi-Civita tensor. Similarly the 
holonomy W p [A] around the boundary of the plaquette p (see Figure [9]) is given by 

W P [A] = 1 + e 2 F p (A) +0(e 2 ). (234) 
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The previous two equations imply that F[N] — lim e ->o X) P TM-^pWpL an d l ea d to the following 
definition: given s, s' € Cyl (think of spin network states) the physical inner product (|23ip is given 

(235) 



s'P,s) :=lim ( s Il/ diV P exp(iTr[7V p iy p ]), s). 



The partition is chosen so that the links of the underlying spin network graphs border the plaque- 
ttes. One can easily perform the integration over the N p using the identity (Peter-Weyl theorem) 



J dN exp(iTr[NW}) = ^ (2j + 1) Tr[n(W)], 



(236) 



where n (W) is the spin j unitary irreducible representation of SU (2). Using the previous equation 



(s'P,s) 



lim 

e-s-0 



n p (e) 

n E(% 



i) ( s 'Tr[n(wg]), s >, 



(237) 



where the spin j p is associated to the p-th plaquette, and n p (e) is the number of plaquettes. 
Since the elements of the set of Wilson loop operators {W p } commute, the ordering of plaquette- 
operators in the previous product does not matter. The limit e — > exists and one can give a 
closed expression for the physical inner product. That the regulator can be removed follows from 
the orthonormality of SU (2) irreducible representations which implies that the two spin sums 
associated to the action of two neighboring plaquettes collapses into a single sum over the action 
of the fusion of the corresponding plaquettes (see Figure ITOl) . One can also show that it is finit J^l. 
and satisfies all the properties of an inner product [101J . 



£(2j + l)(2fc + l) 

jk 



= £(2fc+i) 



Figure 10: In two dimensions the action of two neighboring plaquette-sums on the vacuum is 
equivalent to the action of a single larger plaquette action obtained from the fusion of the original 
ones. This implies the trivial scaling of the physical inner product under refinement of the regulator 
and the existence of a well defined limit e — >• 0. 



14.4 The spin foam representation 

j P 

Each Tr[n {W p )\ in (|23T[) acts in T-Lkin by creating a closed loop in the j p representation at the 
boundary of the corresponding plaquette (Figure 1111 and I13|) . 

Now, in order to obtain the spin foam representation we introduce a non-physical (coordinate 
time) as follows: Instead of working with one copy of the space manifold S we consider n p (e) 

26 The physical inner product between spin network states satisfies the following inequality 

\(s,s>) p \<Cj2W + V 2 - 23 , 
j 

for some positive constant C. The convergence of the sum for genus g > 2 follows directly. The case of the sphere 
g = and the torus 9 = 1 can be treated individually [TOT]. 
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Tv[u(W p )}> j j [_*J j =£AW I J 



Figure 11: Graphical notation representing the action of one plaquette holonomy on a spin network 
state. On the right is the result written in terms of the spin network basis. The amplitude Nj tTn _k 
can be expressed in terms of Clebsch-Gordan coefficients. 

copies as a discrete folliation {Ep}™!^. Next we represent each of the Tr[n (W p )] in ()237j) on 
the corresponding S p . If one inserts the partition of unity between in Hum between the slices, 
graphycally 




(238) 



where the sum is over the complete basis of spin network states {I7, {j})} - based on all graphs 
7 C E and with all possible spin labelling - one arrives at a sum over spin-network histories 
representation of (s,s') p . More precisely, (s',s) p can be expressed as a sum over amplitudes 
corresponding to a series of transitions that can be viewed as the 'time evolution' between the 
'initial' spin network s' and the 'final' spin network s. This is illustrated in the two simple examples 
of Figure [T2l and [T4]) : on the r.h.s. we illustrate the continuum spin foam picture obtained when 
the regulator is removed in the limit e — > 0. 




Figure 12: A set of discrete transitions in the loop-to-loop physical inner product obtained by a 
series of transitions as in Figure [TT] On the right, the continuous spin foam representation in the 
limit e — > 0. 

Spin network nodes evolve into edges while spin network links evolve into 2-dimensional faces. 
Edges inherit the intertwiners associated to the nodes and faces inherit the spins associated to 
links. Therefore, the series of transitions can be represented by a 2-complex whose 1-cells are 
labelled by intertwiners and whose 2-cells are labelled by spins. The places where the action of the 
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Figure 13: Graphical notation representing the action of one plaquette holonomy on a spin network 
vertex. The object in brackets ({}) is a 6j-symbol and Aj := 2j + 1. 



plaquette loop operators create new links (Figure IT3l and IT4]) define 0-cells or vertices. These foam- 
like structures are the so-called spin foams. The spin foam amplitudes are purely combinatorial 
and can be explicitly computed from the simple action of the loop operator in Hkin- 

The physical inner product takes the standard Ponzano-Regge form when the spin network 
states s and s' have only 3-valent nodes. Explicitly, 



Mp= e n (2?/+ i ) ¥ n 



/CF„ 



uCF,. 




(239) 



where the sum is over all the spin foams interpolating between s and s (denoted F s ^ s / , see 
Figure US]), / C F s ^ s r denotes the faces of the spin foam (labeled by the spins jf), v C F s ^ s > 
denotes vertices, and i/f=0if/ns^0A/ns'^0, z// = lif/ns^OV/ns'^0, and v s = 2 if 
/ns = QA/ns' = 0. The tetrahedral diagram denotes a 6j-symbol: the amplitude obtained by 
means of the natural contraction of the four intertwiners corresponding to the 1-cells converging 
at a vertex. More generally, for arbitrary spin networks, the vertex amplitude corresponds to 
3nj-symbols, and (s, s') p takes the same general form. 






Figure 14: A set of discrete transitions representing one of the contributing histories at a fixed 
value of the regulator. On the right, the continuous spin foam representation when the regulator 
is removed. 

Even though the ordering of the plaquette actions does not affect the amplitudes, the spin 
foam representation of the terms in the sum (|239l) is highly dependent on that ordering. This is 
represented in Figur Jl~5l where a spin foam equivalent to that of Figur dT2l is obtained by choosing 
an ordering of plaquettes where those of the central region act first. One can see this freedom of 
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representation as an analogy of the gauge freedom in the spacetime representation in the classical 
theory. 




Figure 15: A different representation of the transition of Figure [T"2l This spin-foam is obtained by 
a different ordering choice in (|237|) . 

One can in fact explicitly construct a basis of T-L p hys by choosing an linearly independent set 
of representatives of the equivalence classes defined in (|232[) . One of such basis is illustrated 
in Figure 1161 The number of quantum numbers necessary to label the basis element is 6<? — 6 
corresponding to the dimension of the moduli space of SU(2) flat connections on a Riemann 
surface of genus g. This is the number of degrees of freedom of the classical theory. In this way we 
arrive at a fully combinatorial definition of the standard H p h y s by reducing the infinite degrees of 
freedom of the kinematical phase space to finitely many by the action of the generalized projection 
operator P. 




Figure 16: A spin-network basis of physical states for an arbitrary genus g Riemann surface. There 
are 6g — 6 spins labels (recall that 4-valent nodes carry an intertwiner quantum number). 



14.5 Spin foam quantization of 3d gravity 

Here we apply the general framework of Sections IBIand loTTTSI This has been first studied by Iwasaki 
(in the spin foam framework) in [2231 1222] . The partition function, Z, is formally given b}Q 

Z = j V[e]V[A] e l ^ Tv[ehF{A)] , (240) 



' We are dealing with Riemannian 3-dimensional gravity. This should not be confused with the approach of 
Euclidean quantum gravity formally obtained by a Wick rotation of Lorentzian gravity. Notice the imaginary unit 
in front of the action. The theory of Riemannian quantum gravity should be regarded as a toy model with no 
obvious connection to the Lorentzian sector. 
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where for the moment we assume M. to be a compact and orientable. Integrating over the e field 
in (|240j) we obtain 

Z = J V[A] 5{F(A)). (241) 

The partition function Z corresponds to the 'volume' of the space of flat connections on M.. 

In order to give a meaning to the formal expressions above, we replace the 3-dimensional 
manifold M. with an arbitrary cellular decomposition A. We also need the notion of the associated 
dual 2-complex of A denoted by A*. The dual 2-complex A* is a combinatorial object defined 
by a set of vertices v G A* (dual to 3-cells in A) edges e £ A* (dual to 2-cells in A) and faces 
/ £ A* (dual to 1-cells in A). The fields e and A have support on these discrete structures. The 
sit(2)-valued 1-form field e is represented by the assignment of an e € su(2) to each 1-cell in A. 
The connection field A is represented by the assignment of group elements g e G SU (2) to each 
edge in A*. 

The partition function is defined by 

Z(A) = f [] de f J] dg e e**Wf\ (242) 

/GA* eSA* 

where det is the regular Lebesgue measure on R 3 , dg e is the Haar measure on SU(2), and Uf 
denotes the holonomy around faces, i.e., Uf = g\ . . . <?f for N being the number of edges bounding 
the corresponding face. Since Uf £ SU(2) we can write it as Uf = 1 + Ff where G C and 
Ff £ su(2). Ff is interpreted as the discrete curvature around the face /. Clearly Tr[e/t//] = 
Tr[e/F/]. An arbitrary orientation is assigned to faces when computing Uf. We use the fact that 
faces in A* are in one-to-one correspondence with 1-cells in A and label e / with a face subindex. 
Integrating over e/, we obtain 

Z(A) = f J] dg e J] %*...<7f), (243) 

eSA* /SA* 

where S corresponds to the delta distribution defined on C 2 (SU(2)). Notice that the previous 
equation corresponds to the discrete version of equation (|241jl . 

The integration over the discrete connection (J\ e dg e ) can be performed expanding first the 
delta function in the previous equation using the Peter- Weyl decomposition 

6(g) = ]T Aj Tr \j(g)] , (244) 

jgirrep(S!7(2)) 

where Aj = 2j + 1 denotes the dimension of the unitary representation j, and j(g) is the corre- 
sponding representation matrix. Using equation (|244p . the partition function (|243p becomes 

Z ( A )= E / II d 9e II ^ f Tr[j f (gl... 9 :)], (245) 

where the sum is over coloring of faces in the notation of (??). 

Going from equation (|242p to (|245[) we have replaced the continuous integration over the e's 
by the sum over representations of SU(2). Roughly speaking, the degrees of freedom of e are now 
encoded in the representation being summed over in (|245[) . 

Now it remains to integrate over the lattice connection {g e }- If an edge e £ A* bounds n 
faces there are n traces of the form Tt\jf (• • • g e • • • )] in (1245|) containing g e in the argument. The 
relevant formula is 

Pin, ■■= J dg ji(g) ® 32(g) ® •• • ® Ug) = J2 c ^- ^ C Z,-- ^ ( 246 ) 
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where P^ nv is the projector onto Inv[ji ® ji ® • ■ ■ <S>j n ]- On the RHS we have chosen an orthonormal 
basis of invariant vectors (intertwiners) to express the projector. Notice that the assignment 
of intertwiners to edges is a consequence of the integration over the connection. This is not a 
particularity of this example but rather a general property of local spin foams. Finally, (|243[) can 
be written as a sum over spin foam amplitudes 



e e n a 

e:{j}-K/} C:{i}^{e} /SA 



(247) 



where A v (l v , j v ) is given by the appropriate trace of the intertwiners l v corresponding to the edges 
bounded by the vertex and j v are the corresponding representations. This amplitude is given in 
general by an SU(2) 3iVj-symbol corresponding to the flat evaluation of the spin network defined 
by the intersection of the corresponding vertex with a 2-sphere. When A is a simplicial complex all 
the edges in are 3-valent and vertices are 4-valent (one such vertex is emphasized in FigureEJ the 
intersection with the surrounding S 2 is shown in dotted lines) . Consequently, the vertex amplitude 
is given by the contraction of the corresponding four 3-valent intertwiners, i.e., a 6j-symbol. In 
that case the partition function takes the familiar Ponzano-Rcggc 303 form 



z(A) = e n a* n 

C:{i}-K/} /6 A * *>e A * 
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(248) 



were the sum over intertwiners disappears since dim(Inv[j! <g> j 2 ® J3]) = 1 for SU(2) and there is 
only one term in (I246[) . Ponzano and Regge originally defined the amplitude (|248|) from the study 
of the asymptotic properties of the 6j-symbol. 



14.5.1 Discretization independence 

A crucial property of the partition function (and transition amplitudes in general) is that it does 
not depend on the discretization A. Given two different cellular decompositions A and A' (not 
necessarily simplicial) 

T- no Z(A) = t-<Z(A'), (249) 

where n is the number of 0-simplexes in A (hence the number of bubbles in J7a), and r = 
is clearly divergent which makes discretization independence a formal statement without 
a suitable regularization. 

The sum over spins in (|248p is typically divergent, as indicated by the previous equation. 
Divergences occur due to infinite volume factors corresponding to the topological gauge freedom 
((2261) fsee pO] H The factor r in (|249|) represents such volume factor. It can also be interpreted 



28 For simplicity we concentrate on the Abelian case G = U(X). The analysis can be extended to the non-Abelian 
case. Writing g £ U(l) as g = e 1 ® the analog of the gravity simplicial action is 

S(A,{e f },{e e })= e/*>({0e», (250) 

/SA* 

where Ff({8 e }) = E e g/ ^ Gauge transformations corresponding to H225H act at the end points of edges e G JTa 
by the action of group elements {/3} in the following way 

Bf -* B f , 

e e ^e e + /3 s -i3t, (251) 
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as a (5(0) coming from the existence of a redundant delta function in (|243[) . One can partially 
gauge fix this freedom at the level of the discretization. This has the effect of eliminating bubbles 
from the 2-complex. 

In the case of simply connected E the gauge fixing is complete. One can eliminate bubbles 
and compute finite transition amplitudes. The result is equivalent to the physical scalar product 
defined in the canonical picture in terms of the delta measur d^l. 

In the case of gravity with cosmological constant the state-sum generalizes to the Turaev- 
Viro model [358] defined in terms of SU q (2) with q n — 1 where the representations arc finitely 
many. Heuristically, the presence of the cosmological constant introduces a physical infrared cutoff. 
Equation (|249p has been proved in this case for the case of simplicial decompositions in 358], see 
also [357, 232 . The generalization for arbitrary cellular decomposition was obtained in [200 . 

14.5.2 Transition amplitudes 

Transition amplitudes can be defined along similar lines using a manifold with boundaries. Given 
A, Ji\ then defines graphs on the boundaries. Consequently, spin foams induce spin networks on 
the boundaries. The amplitudes have to be modified concerning the boundaries to have the correct 
composition property (|74|) . This is achieved by changing the face amplitude from (A^) to (A^) 1 / 2 
on external faces. 

The crucial property of this spin foam model is that the amplitudes are independent of the 
chosen cellular decomposition [3571 1200] . This allows for computing transition amplitudes between 
any spin network states s = (7, {j}, {t}) and s' = (7, {j'}, {V}) according to the following ruleJ^l: 

• Given M. = E X [0, 1] (piecewise linear) and spin network states s — (7, {j}, {(,}) and s' — 
(7j i L '}) on the boundaries - for 7 and 7' piecewise linear graphs in E - choose any 
cellular decomposition A such that the dual 2-complex J7a is bordered by the corresponding 
graphs 7 and 7' respectively (existence can be shown easily). 

• Compute the transition amplitude between s and s' by summing over all spin foam amplitudes 
(rescaled as in (|249jl ) for the spin foams F : s —> s' defined on the 2-complex 

14.5.3 The generalized projector 

We can compute the transition amplitudes between any element of the kinematical Hilbert space 
Transition amplitudes define the physical scalar product by reproducing the skein relations 

where the sub-index s (respectively t) labels the source vertex (respectively target vertex) according to the orientation 
of the edge. The gauge invariance of the simplicial action is manifest. The gauge transformation corresponding to 
(1226ft acts on vertices of the triangulation A and is given by 

B f -> B f + r ls - Vt, 

e e -> e e . (252) 

According to the discrete analog of Stokes theorem 

£ *>({»«}) = o, 

/GBubble 

which implies the invariance of the action under the transformation above. The divergence of the corresponding 
spin foam amplitudes is due to this last freedom. Alternatively, one can understand it from the fact that Stokes 
theorem implies a redundant delta function in J 2 4 3 1) per bubble in j7a- 

29 If M = S 2 X [0, 1] one can construct a cellular decomposition interpolating any two graphs on the boundaries 
without having internal bubbles and hence no divergences. 

30 Here we are ignoring various technical issues in order to emphasize the relevant ideas. The most delicate is 
that of the divergences due to gauge factors mentioned above. For a more rigorous treatment see 13651 1365] . 

31 The sense in which this is achieved should be apparent from our previous definition of transition amplitudes. 
For a rigorous statement see [36511366] . 
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of the canonical analysis. We can construct the physical Hilbert space by considering equivalence 
classes under states with zero transition amplitude with all the elements of H, i.e., null states. 






Figure 17: Elementary spin foams used to prove skein relations. 

Here we explicitly construct a few examples of null states. For any contractible Wilson loop in 
the j representation the state 



1> = (2j + l)a- [J 




s = 0, 

phys 



(253) 



for any spin network state s, has vanishing transition amplitude with any element of %. This 
can be easily checked by using the rules stated above and the portion of spin foam illustrated in 
Figure [T7] to show that the two terms in the previous equation have the same transition amplitude 
(with opposite sign) for any spin- network state in T-L. Using the second elementary spin foam in 
Figure fT71 one can similarly show that 



phys 



(254) 



or the re-coupling identity using the elementary spin foam on the right of Figure [P71 



V2s + l-y/2e + 1 



i j s 
k I e 



(255) 



where the quantity in brackets represents an SU(2) 6j-symbol. All skein relations can be found 
in this way. The transition amplitudes imply the skein relations that define the physical Hilbert 
space! The spin foam quantization is equivalent to the canonical one. 



14.5.4 The continuum limit 

Recently, Zapata 365 , 366 formalized the idea of a continuum spin foam description of 3-dimensional 
gravity using projective techniques inspired by those utilized in the canonical picture |28) . The 
heuristic idea is that due to the discretization invariance one can define the model in an 'infinitely' 
refined cellular decomposition that contains any possible spin network state on the boundary (this 
intuition is implicit in our rules for computing transition amplitudes above). Zapata concentrates 
on the case with non- vanishing cosmological constant and constructs the continuum extension of 
the Turaev-Viro model. 
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14.6 Conclusion 



We have illustrated the general notion of the spin foam quantization in the simple case of 3 
dimensional Riemannian gravity (for the generalization to the Lorentzian case see |172j ). The 
main goal of the approach is to provide a definition of the physical Hilbert space. The example 
of this section sets the guiding principles of what one would like to realize in four dimensions. 
However, as should be expected, there are various new issues that make the task by far more 
involved. 

14.7 Further results in 3d quantum gravity 

In this part of the article we have reviewed SU (2) BF theory as from the pespective of classical and 
quantum gravity in three dimensions (for a classic reference see [106] ). The state sum as presented 
above matches the quantum amplitudes first proposed by Ponzano and Regge in the 1960s based 
on their discovery of the asymptotic expressions of the 6j-symbols 303] and is often referred to as 
the Ponzano-Regge model. Divergences in the above formal expression require regularization. We 
have seen in Section 114.31 that transition amplitudes are indeed finite in the canonical framework 
where M = Sx R. Natural regularizations are available in more general cases [611 12751 1180] , For 
a detailed study of the divergence structure of the model sec [94, 95 , 96 . The quantum deformed 
version of the above amplitudes lead to the so called Turaev-Viro model [358 which is expected 
to correspond to the quantization of three dimensional Riemannian gravity in the presence of a 
non vanishing positive cosmological constant. 

The topological character of BF theory can be preserved by the coupling of the theory with 
topological defects playing the role of point particles. In the spin foam literature this has been 
considered form the canonical perspective in )273| I274J and from the covariant perspective exten- 
sively by Freidel and Louaprc [182 . These theories have been shown by Freidel and Livine to be 
dual, in a suitable sense, to certain non-commutative fields theories in three dimensions [17911178] . 

Concerning coupling BF theory with non topological matter see [1661 1149] for the case of 
fermionic matter, and [341J for gauge fields. A more radical perspective for the definition of 
matter in 3d gravity is taken in |167j . For three dimensional supersymmetric BF theory models 
see [2431132]. 

Recursion relations for the 6j vertex amplitudes have been investigated in 1151] . They 
provide a tool for studying dynamics in spin foams of 3d gravity and might be useful in higher 
dimensions [93] . 
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Part V 

Conceptual issues and open problems 

In this last part of this review we discuss some important conceptual issues that remain open 
questions to a large degree in the formulation. The description of the completely solvable model 
of the previous section will serve as the main example to illustrate some of these issues. 
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14.8 Quantum spacetime and gauge-histories 

What is the geometric meaning of the spin foam configurations? Can we identify the spin foams 
with "quantum spacetime configurations"? The answer to the above questions is, strictly speaking, 
in the negative in agreement with our discussion at the end of Section 12.21 Physical degrees of 
freedom are extracted from the huge set of kinematical ones by the sum over gauge-histories that 
is realized in the spin foam representation. 

This conclusion is illustrated in an exact way the simple example in 2+1 gravity where M = 
S 2 xR (g = 0) described in detail in previous sections. In this case the spin foam configurations 
appearing in the transition amplitudes look locally the same to those appearing in the representa- 
tion of P for any other topology. However, a close look at the physical inner product defined by P 
permits to conclude that the physical Hilbert space is one dimensional - the classical theory has 
zero degree of freedom and so there is no non-trivial Dirac observable in the quantum theory. This 
means that the sum over spin foams in Eq. (|239l) is nothing else but a sum over pure gauge degrees 
of freedom and hence no physical interpretation can be associated to it. The spins labelling the 
intermediate spin foams do not correspond to any measurable quantity. For any other topology 
this still holds true, the true degrees of freedom being of a global topological character. This 
means that in general (even when local excitations are present as in 4d) the spacetime geometric 
interpretation of the spin foam configurations is subtle. This is an important point that is often 
overlooked in the literature: one cannot interpret the spin foam sum of Eq. (|239[) as a sum over 
geometries in any obvious way. Its true meaning instead comes from the averaging over the gauge 
orbits generated by the quantum constraints that defines P - recall the classical picture Figure [2j 
the discussion around Eq. (|8|), and the concrete implementation in 2+1 where U(N) in (|230j) is the 
unitary transformation representing the orbits generated by F. Spin foams represent a gauge his- 
tory of a kinematical state. A sum over gauge histories is what defines P as a means for extracting 
the true degrees of freedom from those which are encoded in the kinematical boundary states. 

Here we are discussing the interpretation of the spin foam representation in the precise context 
of our solvable 3d example; however, the validity of the conclusion is of general character and 
holds true in the case of physical interest: four dimensional LQG. Although, the quantum numbers 
labelling the spin foam configurations correspond to eigenvalues of kinematical geometric quantities 
such as length (in 2+1) or area (in 3+1) LQG, their physical meaning and measurability depend on 
dynamical considerations (for instance the naive interpretation of the spins in 2+1 gravity as quanta 
of physical length is shown here to be of no physical relevance). Quantitative notions such as time, 
or distance as well as qualitative statements about causal structure or time ordering are misleading 
(at best) if they are naively constructed in terms of notions arising from an interpretation of spin 
foams as quantum spacetime configurations. 

Now this does not mean that kinematical features (that usually admit simpler geometric inter- 
pretation) is completely meaningless. It may be that such kinematical quantum geometric consid- 
erations might reflect true physical features when Dirac observables are considered. For instance, 
the simplicial geometry interpretation of spin foam histories (e.g., the total area of a polyhedron 
associated to a spin network vertex) is the analog of the computation of some classical geometric 
quantity in general relativity (e.g. the area of some closed surface embedded in spacetime). Both 
the former and the latter have no intrinsic physical meaning: in the first case as argued above in the 
second case because of the absent diffeomorphism invariant characterization of the surface. How- 
ever, in classical general relativity, the ability of computing area of surfaces is useful to construct 
a physical meaningful quantity such as the area of the event horizon in Schwarzschild spacetime 
which is a fully diffeomorphism invariant property of the BH spacetime. Similarly, all the machin- 
ery of quantum geometry (and its simplicial geometric interpretation in suitable cases) might play 
an important role in the construction of the physically meaningful objects (Dirac observables) of 
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Figure 18: A spin foam as the 'colored' 2-complex representing the transition between three dif- 
ferent spin network states: it represents a pure gauge-history interpolating between kinematical 
boundary quantum geometry states. The sum over such histories in the spin foam representation 
of the path integral is meant to project out gauge degrees of freedom in order to extract the true 
physical information out of the kinematical Hilbert space. 

quantum gravity [3- 

It is well known (and can be illustrated by many simple examples) that gauge invariant proper- 
ties can often be obtained by the suitable combination of gauge covariant quantities. For example, 
in Yang-Mills theory the holonomy h(jAB) € G from a point A to a point B along a path jab 
connecting the two points is not physically meaningful. Nor is the value of some field mutlplets 
*&(A) and 'i'(B) at points A and B respectively. A physically meaningful quantity arises from the 
gauge invariant combination of the above quantities, namely (^(A)h(jAB)^(B)) . The pure gauge 
character of the quantum geometric interpretation pointed out above in the simple 3d context 
corresponds to the holonomy /i(7ab) € G in the above analogy. As soon as we couple 2+1 gravity 
with matter fields (corresponding to the \l/ in the above analogy) non-trivial Dirac observables 
become available (in fact such theory would have local degrees of freedom). These will have to 
encode the relationship between the gauge dependent geometry and the gauge dependent field 
content in a gauge invariant fashion. 

In four dimensions we do not need to add anything as gravity has already its local degrees 
of freedom yet physically meaningful questions are still relational (for more discussion see [] and 
). Unfortunately, no background independent quantum theory with local degrees of freedom is 
sufficiently simple to illustrate this with an example. There is no lucky example as what are 
free quantum field theories for standard QFT (non linearity and background independence may go 
hand- in- hand). In four dimensions, and in particular, in the case of the spin foam models described 
in this review, the difficult question of extracting physics from the models is very important and 

32 An interesting example where this is strictly realized is the midi-superspace formulation of black holes in the 
context of loop quantum gravity known as the isolated horizon formulation (see |27j for a review). In this framework 
the area of what represents the black hole horizon is a Dirac observable and has a discrete spectrum inherited from 
the area spectrum of the kinematical area in the full theory 1221 . Such discreteness id crucial in recovering the 
Bckcnstcin-Hawking area law for the black hole entropy in LQG [1641 121 j . 
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presents a major challenge for the future. The references discussed in Section [13.31 are important 
encouraging steps into this direction. 

14.9 Anomalies and gauge fixing 

As we mentioned before and illustrated with the example of three dimensional gravity (see discus- 
sion in Section I14.8[) the spin foam path integral is meant to provide a definition of the physical 
Hilbert space. Spin foam transition amplitudes are not interpreted as defining propagation in time 
but rather as defining the physical scalar roduct. This interpretation of spin foam models is the 
one consistent with general covariance. However, in the path integral formulation, this property 
relies on the gauge invariance of the path integral measure. If the measure meets this property we 
say it is anomaly free. It is well known that in addition to the invariance of the measure, one must 
provide appropriate gauge fixing conditions for the amplitudes to be well defined. In this section 
we analyze these issues in the context of the spin foam approach. 

Since we are interested in gravity in the first order formalism, in addition to diffeomorphism 
invariance one has to deal with the gauge transformations in the internal space. Let us first describe 
the situation for the latter. If this gauge group is compact then anomaly free measures are defined 
using appropriate variables and invariant measures. In this case gauge fixing is not necessary for 
the amplitudes to be well defined. Examples where this happens are: the models of Riemannian 
gravity considered in this paper (the internal gauge group being SO (A) or SU(2)), and standard 
lattice gauge theory. In these cases, one represents the connection in terms of group elements 
(holonomies) and uses the (normalized) Haar measure in the integration. In the Lorentzian sector 
(internal gauge group SL{2, C)) the internal gauge orbits have infinite volume and the lattice path 
integral would diverge without an appropriate gauge fixing condition. These conditions generally 
exist in spin foam models and are used to regularize the vertex amplitudes of the Lorentzian models 
(we have illustrated this at the end of Section 15731 for a general treatment see [180] ). 

The remaining gauge freedom is diffeomorphism invariance. To illustrate this issue we con- 
centrate on the case of a model defined on a fixed discretization A as it is the usual setting for 
calculations dealing with the models in four dimensions. 

Let us start by considering the spin network states, at <9A*: boundary of A*, for which we 
want to define the transition amplitudes. According to what we have learned from the canonical 
approach, 3-diffeomorphism invariance is implemented by considering (diffeomorphism) equivalence 
classes of spin-network states. In the context of spin foams, the underlying discretization A restricts 
the graphs on the boundary to be contained on the dual 1-skeleton of the boundary complex 
dA*. These states are regarded as representative elements of the corresponding 3-diffeomorphism 
equivalence class. The discretization can be interpreted, in this way, as a gauge fixing of 3- 
diffeomorphisms on the boundary. This gauge fixing is partial in the sense that, generically, there 
will remain a discrete symmetry remnant given by the discrete symmetries of the spin network. 
This remaining symmetry has to be factored out when computing transition amplitudes (in fact 
this also plays a role in the definition of the kincmatical Hilbert space of LQG). 

A natural view point (consistent with LQG and quantum geometry) is that this should naturally 
generalize to 4-diffcomorphisms for spin foams. The underlying 2-complex J7a on which spin foams 
are defined represents a partial gauge fixing for the configurations (spin foams) entering in the path 
integral. The remaining symmetry, to be factored out in the computation of transition amplitudes, 
corresponds simply to the finite group of discrete symmetries of the corresponding spin foamJ^I 
This factorization is well defined since the number of equivalent spin foams can be characterized 
in a fully combinatorial manner, and is finite for any spin foam defined on a finite discretization. 
In addition, a spin foam model is anomaly free if the amplitudes arc invariant under this discrete 

33 Baez 34 points out this equivalence relation between spin foams as a necessary condition for the definition of 
the category of spin foams. 
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symmetry. We have seen that this requirement (advocated in [84]) can be met by suitable definitions 
of the transition amplitudes {4"51 1229] (see Section [73]) . 

However, it is expected - from the experience in lower dimensional background independent 
models (recall the discussion of the previous section) - that there will be remnants of the gauge 
symmetries acting non trivially on spin foams histories by doing more than simply changing the 
embedding in the discrete regulating structure A as described above. After all, in gravity the 
Hamiltonian constraint generates gauge transformations that hide in themselves the non trivial 
dynamics of the theory. This is illustrated precisely in the simplest scenario of 3d gravity where 
it is well known that bubble divergencies are directly linked to the infinite volume of the gauge 
orbits generated by the curvature constraint |180j . The action of such gauge symmetry related 
spin foams that do not differ only by their embedding (see Section fT4|) . 

It is also important to point out that the view that spin foams represent diffcomorphism equiv- 
alence classes of geometries, which sometimes can be loosely stated in the community, is incom- 
patible with the idea that spin foams should act as a projector onto the solutions of constraints. 
A path integral would only act as a projector if it includes gauge symmetries in the sense of group 
averaging. This related to the fact that physical states (those left invariant by the projector) are 
outside the kinematical Hilbert space, i.e. non-normalizable. 

The discretization of the manifold A is seen as a regulator introduced to define the spin foam 
model. Even when the regulator (or the discretization dependence) eventually has to be removed 
(see next subsection), the theory is presumed to remain discrete at the fundamental level. The 
smooth manifold diffcomorphism invariant description is expected to emerge in the (low energy) 
continuum limit. Fundamental excitations are intrinsically discrete. From this viewpoint, the 
precise meaning of the gauge symmetries of such a theory would have to be formulated directly 
at the discrete level. We have seen that this can be achieved in the case of 3-dimensional gravity 
(recall Section EXT) . 

The previous paragraph raises the question of whether the discretization procedures used in the 
derivation of the spin foams are compatible with the expectation that one is approximating a diffeo- 
morphism invariant fundamental theory. More precisely, can one tell whether the diffeomorphism 
invariance is broken or not by our regularization procedure. This question is a quite important one 
and has been one of the central concerns of the work led by B. Dittrich and collaborators for the 
last few years [1451 11411 11441 11381 [39] . There exist indeed possible discretizations of a field theory 
which maintain the full symmetry content of their continuum relatives. These regularizations are 
called perfect actions and even when they are difficult to construct explicitly in general cases there 
are in principle methods (based on renormalization group ideas) to approach them [1381 140] (see 
also [H] and [12] for more recent discussion and explicit calculations in the perturbative context). 

There are however indications that the discretization procedures used generically do break the 
general covariance of gravity. In the simple case of 3d gravity with cosmological constant this can 
be shown explicitly at the classical level [UJ as well as at the quantum level 298 . In the first 
reference it is shown nonetheless that there exist a well defined perfect action in this case. At 
the quantum level the results of [276] indicate a possible resolution of the anomaly issue. In four 
dimensions, this was first validated by the results obtained by Gambini and Pullin et al. in the 
so-called consistent discretization approach [184] [192] [193] HUJ H90] EDS] [132] [189] OH] [THE] . They 
study the canonical formulation of theories defined on a lattice from the onset. The consistent 
discretization approach provides a way to analyzing the meaning of gauge symmetries directly a la 
Dirac. Their results indicate that diff-invariance is indeed broken by the discretization in the sense 
that there is no infinitesimal generator of diffeomorphism. This is consistent with the covariant 
picture of discrete symmetries above. In their formulation the canonical equations of motion fix 
the value of what were Lagrange multipliers in the continuum (e.g., lapse and shift). This is 
interpreted as a breaking of diffeomorphism invariance; however, the solutions of the multiplier 
equations are highly non unique. The ambiguity in selecting a particular solution corresponds 
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to the remnant diffeomorphism invariance of the discrete theory. More recently, the issue of the 
breaking of the diffeomorphism symmetry by the regularizations used in spin foams has been 
studied by the Dittrich group (see references above) . A possibility of using very simple models to 
study these questions has been open in [43] . 

An important question is whether the possible breaking of the diffeomorphims gauge symmetry 
by the regularization can be made to disappear in the continuum limit. The breaking of a gauge 
symmetry implies the existence of a new host of degrees of freedom that might be in conflict 
with the low energy interaction one is trying to recover in that limiiF^I. In addition the breaking 
of non-compact symmetries by a regulating structure poses serious difficulties for the low energy 
regime even in the case of global symmetries |109j (see |194| ?] and [302] for a recent discussion). 
We mention this point here because it is in our opinion of a major importance. Even though the 
question seems quite difficult to address in the context of the present models at this stage it must 
be kept in mind in search of opportunities of further insights. 

However, the notion of anomaly freeness evoked at the beginning of this section should be 
strengthened. In fact according to our tentative definition, an anomaly free measure can be multi- 
plied by any gauge invariant function and yield a new anomaly free measure. This kind of ambiguity 
is not wanted; however, it is in fact present in most of the spin foam models defined so far: In stan- 
dard QFT theory, the formal (phase space) path integral measure in the continuum has a unique 
meaning (up to a constant normalization) emerging from the canonical formulation. Provided an 
appropriate gauge fixing, the corresponding Dirac bracket determines the formal measure on the 
gauge fixed constraint surface. In the discussion in Section [12. II we have provided references where 
this issue is analyzed. 

14.10 Discretization dependence 

The spin foam models we have introduced so far are defined on a fixed cellular decomposition of 
Ai. This is to be interpreted as an intermediate step toward the definition of the theory. The 
discretization reduces the infinite dimensional functional integral to a multiple integration over a 
finite number of variables. This cutoff is reflected by the fact that only a restrictive set of spin 
foams (spin network histories) is allowed in the path integral: those that can be obtained by 
all possible coloring of the underlying 2-complex. In addition it restricts the number of possible 
3-geometry states (spin network states) on the boundary by fixing a finite underlying boundary 
graph. This represents a truncation in the allowed fluctuations and the set of states of the theory 
playing the role of a regulator. However, the nature of this regulator is fundamentally different 
from the standard concept in the background independent framework: since geometry is encoded 
in the coloringfthat can take any spin values) the configurations involve fluctuations all the way 
to Plank scala 35 !. This scenario is different in lattice gauge theories where the lattice introduces an 
effective UV cutoff given by the lattice spacing. Transition amplitudes are however discretization 
dependent now. A consistent definition of the path integral using spin foams should include a 
prescription to eliminate this discretization dependence. 

A special case is that of topological theories such as gravity in 3 dimensions. In this case, 
one can define the sum over spin foams with the aid of a fixed cellular decomposition A of the 
manifold. Since the theory has no local excitations (no gravitons), the result is independent of the 
chosen cellular decomposition. A single discretization suffices to capture the degrees of freedom of 
the topological theory. 

34 Breaking gauge symmetries comes hand-in-hand with introducing (unwanted) new degrees of freedom. For 
example, longitudinal photon modes would appear if U(l) gauge symmetry is violated in electromagnetism; similar 
spurious modes come to life if diffeormorphism invariance is broken in gravity theories. 

35 Changing the label of a face from j to j + 1 amounts to changing an area eigenvalue by an amount of the order 
of Planck length squared according to 171> . 
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In lattice gauge theory the solution to the problem is implemented through the so-called con- 
tinuum limit. In this case the existence of a background geometry is crucial, since it allows one to 
define the limit when the lattice constant (length of links) goes to zero. In addition the possibility 
of working in the Euclidean regime allows the implementation of statistical mechanical methods. 

None of these structures are available in the background independent context. The lattice (tri- 
angulation) contains only topological information and there is no geometrical meaning associated 
to its components. As we mentioned above this has the novel consequence that the truncation 
cannot be regarded as an UV cutoff as in the background dependent context. This in turn repre- 
sents a conceptual obstacle to the implementation of standard techniques. Moreover, no Euclidean 
formulation seems meaningful in a background independent scenario. New means to eliminate 
the truncation introduced by the lattice need to be developed. For a recent analysis where the 
difference between lattice regularization of background dependent versus background independent 
theories is carefully considered see [320] . For a more general discussion see |318j . 

This is a major issue where concrete results have not been obtained so far beyond the topological 
case. Here we explain the two main approaches to recover general covariance corresponding to the 
realization of the notion of 'summing over discretizations' of |314j . 

• Refinement of the discretization: 

According to this idea topology is fixed by the simplicial decomposition. The truncation 
in the number of degrees of freedom should be removed by considering triangulations of 
increasing number of simplexes for that fixed topology. The flow in the space of possible 
triangulations is controlled by the Pachner moves. The idea is to take a limit in which the 
number of four simplexes goes to infinity together with the number of tetrahedra on the 
boundary. Given a 2-complex J2 which is a refinement of a 2-complex J\ then the set of 
all possible spin foams defined on J\ is naturally contained in those defined on J2 (taking 
into account the equivalence relations for spin foams mentioned in the previous section) . The 
refinement process should also enlarge the space of possible 3-geometry states (spin networks) 
on the boundary recovering the full kinematical sector in the limit of infinite refinements. An 
example where this procedure is well defined is Zapata's treatment of the Turaev-Viro model 
[365] . The key point in this case is that amplitudes are independent of the discretization 
(due to the topological character of the theory) so that the refinement limit is trivial. In 
the general case the definition of the refinement limit has been recently been studied and 
formalized in [32 1J . There is no evidence as to whether the amplitudes of any of the present 
models for 4d converges in such limit. 

It has been often emphasized that such refinement limit may be studied from the Wilsonian 
rcnormalization view point. In the past a renormalization approach for spin foams have been 
proposed by Markopoulou [260, 259 . Also Oeckl |279j has studied the issue of renormal- 
ization in the context of spin foam models containing a coupling parameter. These models 
include generalized covariant gauge theories |278[ 13001 1285] . the Reisenberger model, and the 
so called interpolating model (defined by Oeckl). The latter is given by a one-parameter 
family of models that interpolate between the trivial BF topological model and the Barrett- 
Crane model according to the value of a 'coupling constant'. Qualitative aspects of the 
renormalization groupoid flow of the couplings are studied in the various models. 

• Spin foams as Feynman diagrams: 

This idea has been motivated by the generalized matrix models of Boulatov and Ooguri 
[99| 1286] . The fundamental observation is that spin foams admit a dual formulation in terms 
of a field theory over a group manifold |127l 13091 1308] . The duality holds in the sense that spin 
foam amplitudes correspond to Feynman diagram amplitudes of the GFT. The perturbative 
Feynman expansion of the GFT (expansion in a fiducial coupling constant A) provides a 
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definition of sum over discretizations which is fully combinatorial and hence independent of 
any manifold structure!^!. The latter is most appealing feature of this approach. For a recent 
reference proposing a GFT for the EPRL type of models see [H] . 

However, the convergence issues clearly become more involved. The perturbative series are 
generically divergent. This is not necessarily a definite obstruction as divergent series can 
often be given an asymptotic meaning and provide physical information. Moreover, there 
are standard techniques that can allow to 're-sum' a divergent series in order to obtain non 
perturbative information [3 1 3J . Freidel and Louapre 181 have shown that this is indeed 
possible for certain GFT's in three dimensions. Other possibilities have been proposed in 
[309]. 

Diffeomorphism equivalent configurations (in the discrete sense described above) appear at 
all orders in the perturbation seriea^. From this perspective (and leaving aside the issue of 
convergence) the sum of different order amplitudes corresponding to equivalent spin foams 
should be interpreted as the definition of the physical amplitude of that particular spin foam. 
The discussion of the previous section does not apply in the GFT formulation, i.e., there is 
no need for gauge fixing. 

The GFT formulation could resolve by definition the two fundamental conceptual problems 
of the spin foam approach: diffeomorphism gauge symmetry and discretization dependence. 
The difficulties are shifted to the question of the physical role of A and the convergence of 
the corresponding perturbative series. A prescription that avoids this last issue is the radical 
proposal of Freidel [17.3] where the tree-level GFT amplitudes are use to define the physical 
inner product of quantum gravity. 

In three dimensions this idea has been studied in more detail. In [256] scaling properties of 
the modification of the Boulatov group field theory introduced in |181j was studied in detail. 
In a further modification of the previous model (known as coloured tensor models [202 ) 
new techniques based on a suitable 1/N expansion imply that amplitudes are dominated 
by spherical topology [201] ; moreover, it seem possible that the continuum limit might be 
critical as in certain matrix models [204 , 89, 203. 205 , 332 . However, it is not yet clear if 
there is a sense in which these models correspond to a physical theory. The appealing possible 
interpretation of the models is that they correspond to a formulation of 3d quantum gravity 
including a dynamical topology. No much is known about gravity models in four dimensions 
at this stage. 

Finally, discretization independence is not only a problem of spin foams but also of Regge 
gravity. Even linearized classical Regge calculus is not discretization independent in 4d as shown in 
|146] where a general discussion of the issues raised in this section are considered. It is also natural 
to expect the issue of discretization independence to be connected to the issue of diffeomorphism 
symmetry discussed in the previous section. There is evidence [HJ 11401 1139] that implementation 
of diffeomorphism symmetry into the models will lead to discretization independence. 

14.11 Lorentz invariance in the effective low energy regime 

Finally a very important question for loop quantum gravity and spin foams (closely related to the 
consistency of their low energy limit) is that of the fate of Lorentz symmetry in the low energy 
regime. This question has received lots of attention in recent years as it seems to provide an open 

36 This is more than a 'sum over topologies' as many of the 2-complex appearing in the perturbative expansion 
cannot be associated to any manifold 11281 . 

37 The GFT formulation is clearly non trivial already in the case of topological theories. There has been attempts 
to make sense of the GFT formulation dual to BF theories in lower dimensions |251j . 
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"window of opportunity for observations of quantum gravity effects (see |264[ 1227] and references 
therein). In the context of loop quantum gravity a systematic first-principles analysis of such 
properties of the low energy regime remains open. This is mainly due to the presently unsolved 
difficulties associated with a quantitative dynamical description of physical states. The spin foam 
approach discussed in this review is meant to address the dynamical question and the definition of 
the new models may provide the necessary insights to tackle the relevant aspects for the question 
of Lorentz invariance at low energies. This question is perhaps within reach at the present stage 
of development of the theory, so it deserves all the attention in future investigations. 

In the absence of detailed formulations, certain early pioneering model calculations |187j (see 
also general treatment [333, 334 ) indicated the possibility for certain Lorentz invariance violating 
(LIV) effects associated to a quantum gravity granularity of spacetime measured at rest by a 
preferred frame of observers. The effects take the form of corrections to the effective Lagrangian 
with LIV terms of dimension 5 or higher that are suppressed by negative powers of the Planck 
mass. This makes them in principle negligible at low energies so that the potential quantum gravity 
effects would be observed as corrections to the standard QFT well tested physics. However, it was 
soon realised [109] that the framework of effective field theories severely restricts the possibility 
that LIV remains small when interactions in QFT are taken into account. In fact dimension 5 
or higher LIV terms will generically (in the absence of some fundamental protecting mechanism) 
generate dimension 4 an lower LIV terms (with factors in front that are at best quadratic in 
standard model coupling constants). Such unsuppressed contributions are in flagrant conflict with 
the observed Lorentz invariance of particle physics at low energies. 

Therefore, the early calculations in LQG cited above need to be revised under the light of the 
strong constraints on the way LIV can arise in quantum gravity. A logical possibility is that these 
LIV terms are there thanks to the protecting effect of some custodial mechanism restricting the 
LIV to dimension 5 and higher operators 0. For example such possibility has been considered in 
[194] where it is shown that the large LIV effects can be avoided in the context of Euclidean field 
theory (the protecting mechanism being here the trading of the non-compact Poincare group by 
the Euclidean group). This example is interesting but unfortunately does not provide an answer 
to the question for the physical Lorentzian theory (see discussion in [108] and more recently [302] ). 

There is thus a challenge for LQG and spin foams consisting of showing that the low energy 
limit of the theory is compatible with the observed low energy LI. This can be put in the form of 
two conservative possibilities. The description of the low energy physics in terms of 

1. an effective action with no LIV terms, or 

2. an effective action where LIV terms appear starting from dimension 5 operators on due 
to some fundamental protecting mechanism in LQG for such structure to be stable under 
radiative corrections. 

My view, based on the ignorance of a suitable protecting mechanism in LQG, is that only the 
first possibility is sensible. This would mean that the way the fundamental Planckian discreteness 
of LQG must enter physics is through reference frame independent phenomena (as mass enters 
a Klein-Gordon theory). In my opinion the main ingredient that must be taken into account in 
considering this issue is the construction of the low energy limit of quantum gravity in terms of 
physical or gauge invariant notions (see discussion in Section fl4.8p . The LIV effects obtained in 
LQG so far promote by assumption the kinematical discreteness of quantum geometry to discrete- 
ness of the physical background in a literal fashion (on the one hand discreteness of geometric can 
coexist with Lorentz invariance [324] , on the other hand, at the level of gauge invariant observables 

38 For an example of protecting symmetry see |107j : however, notice that the non-compact nature of the protecting 
symmetry in this case makes it seemingly unsuitable for a sensible possibility in LQG. 
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some discrete aspect could disappear [148] ). The limitations of such naive view point on kinemat- 
ical discreteness versus physical properties in the theory is, in my view, best manifested in the the 
setting of pure 3d gravity as discussed above. 

For completeness we add that in addition to the above (conservative) possibilities, some less 
standard scenario could be realized. This has been the view point of some part of the quantum 
gravity community. For a recent proposal of a completely new theoretical framework see [18] and 
references therein. 
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